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The probabilities of auto-ionization of the 2s*'S, 2s2p*P, 2p*'1D, 3d*'G states of He are 
calculated with continuous wave functions which are obtained by numerical integration of the 
wave equation including the effect of exchange. The probabilities are 4 10", 5 104, 1 10", 
4.9X 10" per second, respectively, for these states. On the basis of these calculations and the 
theoretical energies of the doubly excited states, the identification of the two lines 4320.4 and 
4357.5 observed by Compton and Boyce and by Kruger is discussed. The former has been 
reasonably definitely ascribed to the transition 1s2p*P—2p**P, the state 2p**P not being 
subject to auto-ionization. The latter line cannot be satisfactorily accounted for with respect to 


both the position and the width of the line. 





HERE have been many investigations on 
different aspects of the problem of doubly 
excited helium. Compton and Boyce,'™ and later 
Kruger,!® observed in the extreme ultraviolet 
region two lines \320.4 (v=312,117 cm™@) and 
\357.5 (v= 279,715 cm~) which they ascribed to 
transitions from the doubly excited to the ordi- 
nary singly excited states of helium. Considerable 
interest on the energy states of doubly excited 
helium has been aroused by the suggestion, now 
abandoned, that the corona lines may be due to 
transitions among the doubly excited states of 
helium.? This has led to the theoretical calculation 
of these energy states by many authors.* Many 
attempts‘ have been made to observe the spec- 
1 (a) K. T. Compton, and Boyce, J. Frank. Inst. 205, 497 
(1928); (b) P. Gerald Kruger, Phys. Rev. 36, 855 (1930). 

* Rosenthal, Zeits. f. Astrophys. 1, 115 (1932); Goudsmit 
and Wu, Astrophys. J. 80, 154 (1934). 

*F. G. Fender and J. P. Vinti, Phys. Rev. 46, 77 (1934); 
T. Y. Wu, ibid. 46, 239 (1934); T. Y. Wu and S. T. Ma, 
ibid. 48, 917 (1935); Chin, J. Chem. Soc. 4, 344 (1936); W. 
Wilson, Phys. Rev. 48, 536 (1935). 

‘Kiang, Ma, and Wu, Phys. Rev. 50, 673 (1936); Chin. 


a 2, 117, (1936); F. Bundy, Phys. Rev. 52, 452 


trum of doubly excited helium under experi- 
mental conditions that would have been most 
favorable for their excitation according to the 
theoretical excitation probabilities of these states 
by electron impact.’ All these attempts have 
failed. Approximate calculations of the natural 
lifetime of the 2s3s*S state as the result of auto- 
ionization have been carried out by Kreisler* and 
by Kiang, Ma, and Wu.‘ They found a natural 
lifetime of the order 10-'*-10-"* second. 

It was pointed out by Kiang, Ma, and Wu that 
the failure to observe the spectrum of doubly 
excited helium may be understood on the basis 
of the theoretical great natural width of these 
states, since the spreading of a “‘line’’ over a 
spectral region of the order 100-1000 angstroms 
renders its observation exceedingly difficult. It 
was pointed out, however, by these authors that 
on accepting this theoretical probability of auto- 
ionization and on assuming the probabilities for 


5 Massey and Mohr, Proc. Camb. Phil. Soc. 31, 604 


(1935). 
® Kreisler, Acta Physica Polonica 4, 15 (1935). 
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TABLE I. Probabilities of auto-ionization. 











State i State f Atomic wave functions Ei in Rh P in sec." 
2s?1§ Isks 4S v2e~(1—Br)e*", B=0.8469 — 1.4440 4X 10" 
2s2p*P iskp*P voe~(1—Br)e*", Yop~re-7, B=0.874, y = 0.862 — 1.5043 5X 10% 
2p?'!D iskd }D Vop~re—7", y=0.815 — 1.328 1X10" 
3€7'G iskg 1G vsa~r*e*", 5 = 0.659 — 0.868 4.2 10" 
3d? 'G 2skg 1G vsa~re*", 5 =0.659 + 0.868 4.9X 10" 








the other doubly excited states to be of the same 
order of magnitude, one would expect a width of 
the order 1 angstrom for the lines \320.4 and 
\357.5 whereas the observed lines appear to be as 
sharp as the ordinary lines of He II. The situation 
is hence an unsatisfactory one. Either the calcu- 
lated probability of auto-ionization for 2s3s*S is 
grossly in error and the assumption that all the 
doubly excited states have about the same life- 
time is not justified, or an explanation of the 
\320.4 and A357.5 lines must be sought elsewhere. 
The purpose of the present study is to make more 
calculations for a few doubly excited states of 
helium and to discuss the possible identification 
of the \320.4 and A357.5 lines. 


PROBABILITY OF AUTO-IONIZATION. 


The general theory of the process of auto-ioni- 
zation is a well-known one.’ Let V;, V; be the wave 
functions of the initial (a doubly excited atom, 
for example) and the final (an ion plus an ejected 
electron) state. The probability of auto-ioniza- 
tion is then 


4x? 4r? * é, 
a rapatt| foils) 
h h kj 


where the integration is taken over the coordi- 
nates of all the electrons. For a triplet or a singlet 
state, we have, respectively, 


2 
’ 


(1) 





1 


1 
Wy = Wel ve(2) Fve(Ive(2)), (3) 


where ya, ¥, ¥- are the discrete atomic wave 
functions and y; is the continuous wave function 
of the ejected electron. Kreisler employed for yx 


= and Shortley, Theory of Atomic Spectra (1935), 
p. 369. 





the continuous wave function of hydrogen, 
thereby approximating the atomic field due to 
the nucleus and an electron by that of a bare 
nucleus of charge 1; whereas Kiang, Ma, and Wu 
employed, for 2s3s*S, 


1/8xm\' 1 . 
n=-(—") — cos ( r—a). (4) 


r 


This amounts to employing for y its asymptotic 
form for large r which does not hold for small r, 
and with a phase a@ so chosen that (4) is not valid 
for large r. Both calculations neglected the effect 
of exchange on the solution y of the wave equa- 
tion. As the value of the matrix element Vj; in 
(1) depends rather strongly on the overlap of the 
various wave functions, it is thought that the 
insufficiently accurate y of the earlier calcula- 
tions may cause considerable error in the proba- 
bility of auto-ionization. In the present calcula- 
tion we shall seek a solution y, taking into 
account the effect of exchange and the actual 
atomic field. 

The function W; is the solution of the 
Schrédinger equation 


2Z 2Z 2 
| atast e+ —+— —=}y,=0, (5) 


és ee 


and y, is the solution of 


r 


2Z 
[a+Et+—p.=0, (6) 


where energy is in units of the ionization energy 
of hydrogen and distance in units of the first 
Bohr radius of hydrogen. Substitution of (3) into 
(5) and making use of (6) give the integro- 
differential equation® 


triplet 
[A+k?— V(r) ](ryir) = F U(r) (re) (7) 


singlet 


8 Cf. P. M. Morse and W. P. Allis, Phys. Rev. 44, 269 


(1933). 
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in which 


1 Z 
Vrs) =2 f (—-—)veraar (8) 
T12 =e 


and 


1 
-Ulrs) =2 f blr adr dee: 


Tie 


~(k*—E,) f VerWalr)dr. (9) 


k? is the kinetic energy of the ejected electron so 
that E=E.+k?*. V is the energy of the ejected 
electron in the atomic field and U(r) represents 
the effect of exchange. The subscript 7 is the 
azimuthal quantum number of the ejected elec- 
tron. For a given initial state V;(mu, la; mp, 1) and 
a given ¥.-(m.,/.), the value of / is determined by 
the following rule for the non-vanishment of the 
matrix element V;; and hence for the presence of 
auto-ionization : The states 7 and f must (7) be of 
the same parity, (77) have the same value of J, 
and in the case of L, S coupling, (iz) have the 
same values of L and S. To obtain the solution 
of (7), since the effect of exchange is expected to 
be small, we shall, as an initial approximation, 
obtain the solution of the homogeneous equation 


[A+k?— V(r) ](ryer) =0 (10) 


; . ; , 0 

by numerical integration. The solution y¥,: so 
obtained is substituted into (9) and Eq. (7) is 
solved again by numerical integration. If the 


° (1) ° . ‘ 
solution ¥,; so obtained differs considerably from 
ree ae , ; 
the initial ¥x:, Eq. (7) is solved again by using 
a. : , . 
vx: in U(r). The process is repeated until the ith 


solution dar of (7) gives a value for V;, differing 
by less than 20 percent from the one obtained 


with the previous solution -_—: It is found that 
the successive approximations oscillate about 
and converge upon the final solution very 
rapidly. 

For the discrete wave functions, the simple 
variational wave functions of the earlier work? 
are employed. The continuous wave function is 
normalized according to (4) where cos (kr—a) is 
the asymptotic form of rpg. 

As an example, the equation for the radial part 
of the continuous wave function for the auto- 


AUTO-IONIZATION 





IN HE 


ionization of 2s2p*P is 


d? 2 2 2 
— (rer) + |2.49574+-+ (4 += Jew —_ = |orray 
dr? r r r? 
64 1 r 
= - et — f re~?" (rv) rdr 
3 r? 0 


» 5 
+rf <*(rburir| 
r r? 


The probabilities of auto-ionization of a few 
states of doubly excited helium calculated in this 
manner are given in Table I. 

Concerning the magnitude of the probability 
of auto-ionization in general, very little can be 
said on the basis of these calculations. It appears, 
however, from the values of P for the auto- 
ionization of 3d? 'G into 1s Xg'G and 2s Xg'G that 
P decreases very rapidly as the overlap of the 
atomic wave functions in the initial and the final 
states decreases. 


THE 2320.4 AND 4357.5 LINES OF HELIUM 


As the above calculations seem to indicate 
a probability of auto-ionization of the order 
10'*-10'* per second for the doubly excited states 
of helium, and as the observed lines \320.4 and 
\357.5 do not seem to show any diffuseness, we 
would first try to identify them with transitions 
from such states as not being subject to auto- 
ionization. Of the configurations 2sns, 2snd, 2pnp, 
3dnd, 2snp, 2pnd, only the states 2pnp'P, *P; 
2pnd'D, *D; 3dnd'D, *D do not undergo auto- 
ionization. The energy of 2p? *P has been calcu- 
lated by Wilson by the method of self-consistent 
field and is E(2p?*P)=—1.4018 Rh measured 
from the state of the naked nucleus.* This value 
is lower than that obtained by the use of simple 
variational wave functions, but a still better 
approximation to the true energy value of 2p? *P 
can be made by means of the method of Steven- 
son.® This method consists in using for the varia- 
tional wave function the function (1+¢ cos 612) ¥ 
where W is the correct combination of one- 
electron wave functions determined by the 
method of self-consistent field, and 62 is the angle 
between 7; and rz of the two electrons. Applica- 


* A. F. Stevenson, Phys. Rev. 56, 586 (1939). 























tion of the variational principle to the determi- 
nation of the parameter c leads to the secular 
equation for the energy of the system. For our 
purpose, we shall evaluate the integrals involved 
in the theory by means of the variational wave 
functions of Wu and Ma instead of Wilson’s self- 
consistent field solutions. Since the energy cor- 
rection given by Stevenson’s method is small in 
any case, this procedure is completely justified. 

Using the notation of Stevenson, and with 
Yop~re—” where y = 0.836, we find 


Hu=77/192, Hiu=4y2/5, ME=—0.0091 Rh. 


Thus the corrected energy of 2p? *P is —1.4018 
— 0.0091 = —1.4109 Rh. With this value, the 
calculated frequency of the line 1s2p*P—2p? *P 
is 313,348 cm which differs by 0.4 percent from 
the observed 312,117 cm— for the line 4320.4. As 
the true value of E(2p? *P) lies probably some- 
what lower than —1.4109 Rh, one may feel 
reasonably sure of the identification of (320.4 
with the transition 1s2p*P —2p? *P. 

With the other line \357.5 (v=179,715 cm~), 
the situation is far less satisfactory. From the 
calculated energy states of doubly excited helium, 
it is found that only the lowest one, namely, 
2s2p*P, lies within 280,000 cm from the ioniza- 
tion limit of He I, the other states all being too 
high. On assuming the line \357.5 to originate 
from 2s2s*P, it may then be compared with the 
following transitions 


154s*S—2s2p*P = 282,120 cm; 
1s4d*D — 2s2p*P = 280,968 cm. 


But then the auto-ionization of 2s2p*P would 
broaden the line to about 0.3A which is ap- 
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parently too great compared with the observa- 
tion. Also it is not clear why the transition 
1s2s*S—2s2p*P, which is expected to be stronger 
than the two above, did not appear on the same 
spectrogram. Hence it appears that, unless the 
calculated energies of the doubly excited states 


and the probabilities of auto-ionization are seri- : 


ously in error, the line \357.5 cannot be satis- 
factorily explained. 


SATELLITE LINES IN CV AND CVI 
OBSERVED BY EDLEN 


Edlén and Tyrén observed groups of weak 
lines on the long wave-length side of the lines 
1s?4\S—1s2p'P, 1s°S—2p?P in C V and C VI, re- 
spectively.’° These satellite lines they ascribed to 
1s?2s°S—1s2s2p?P, 1s°2p°P —1s2p*P, *D, *S, etc., 
in C IV, and 1s2s**S—2s2p'*P, 1s2p+%pP 
—2p2p'P,'D,'S in C V, respectively. An attempt 
has been made by the writer to correlate the 
individual satellite lines with such transitions on 
the basis of energies of such doubly excited states 
calculated by the variational method." Most of 
these doubly excited states are capable of auto- 
ionization and accordingly the satellite lines 
should show some broadening. If, however, the 
probabilities of auto-ionization of such states as 
1s2p? 2D are of the same general order of magni- 
tude as in He, namely, 10'*-10"* per second, the 
width of these lines in the region of 40A would be 
0.001-0.01A, which is not hopelessly inconsistent 
with observation, in distinction to the case of the 
4357.5 line of He. 


1” Edlén and Tyrén, Nature 143, 940 (1939). 
u T. Y. Wu, Phys. Rev. 58, 1114 (1940). 
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Further Studies on the Origin of Cosmic Rays 
Helium Annihilation Rays and the Cause of Their Variability with Time 


RosBert A. MILurKan, H. Victor NEHER, AND WILLIAM H. PICKERING 
California Institute of Technology, Pasadena, California 


(Received October 16, 1944) 


We find definite evidence that a new band of rays which we interpret as helium annihilation 
rays does come in vertically at about the predicted latitude. We present a discussion of the 
possible composite character of the so-called silicon annihilation band and of the so-called oxygen 
annihilation band. We then bring forward an explanation of the cause of our large and already 
reported variability in the cosmic-ray intensities found in high altitude electroscope flights at 
Bismarck, Omaha, and Oklahoma City. We also make a new and more accurate determination 
of the value of the field sensitive and the non-field sensitive components of the incoming 


cosmic rays. 





I. HISTORICAL SUMMARY OF THE DISCOVERY OF 
THE SILICON COSMIC-RAY BAND 


N our expedition to India in late 1939 and early 
1940 we had found unambiguous evidence 
that in going from the magnetic equator near 
Bangalore north to Peshawar (Mag. Lat. 25° N) 
there was no change at all in the intensity of the 
vertically incoming cosmic rays between Banga- 
lore and Agra (Mag. Lat. 17.3), but in going the 
next 7.7°, from Agra to Peshawar, there was a 
marked increase in that intensity, amounting 
to more than 20 percent. 

We had gone to India primarily to look for 
just this cosmic-ray effect, for we had predicted 
on the basis of the hypothesis that the birth 
of the cosmic rays in interstellar space consists 
merely in the complete transformation into a 
“charged particle pair’’ (which we treat as an 
“electron pair’’') of the total rest-mass energy of 
the most abundant of the atoms found in inter- 
stellar space, namely, the atoms of silicon, 
oxygen, nitrogen, carbon, and helium. 

If this were so we had computed from the 
Einstein equation E=mc* that each electron of 
the ‘electron pair’’ arising from such a trans- 
formation of the silicon atom (we have called 

1A mesotron pair is barred out by its short lifetime. The 
only choice is between an electron pair and a proton pair, 
but the difference between the latitude of entrance of 
electrons and protons entering the earth’s magnetic field 
from this mode of origin is in no case, not even in the case of 
the He annihilation rays, large enough to be detected with 
the resolving power of the experimental techniques we have 
so far used. However, we hope to obtain a resolution which 
will differentiate between these alternatives. See footnote 4, 
ao Rev. 61, 398 (1942). We miscalculated the difference 
mn 


hys. Rev. 63, 245 (1943), so that the conclusion there 
reached is not yet{a valid one. 


these ‘‘silicon annihilation rays’’) would have an 
energy of 13.2 billion electron volts (Bev), not 
enough to break vertically through the earth’s 
magnetic field at the magnetic equator in India, 
where we had computed that the requisite in- 
coming charged particle energy was 17 billion 
electron volts, but it would have just about 
enough incoming energy to break through at 
somewhere near 20° N Magnetic Latitude, i.e., 
somewhere between Agra and Peshawar. 

Having then found experimentally by balloon 
flights to very high altitudes a plateau of con- 
stant vertically incoming cosmic-ray energy be- 
tween Bangalore and Agra and a new band of 
rays coming in between Agra and Peshawar, we 
provisionally identified this new band as the 
cosmic-ray silicon annihilation band. 

We found additional evidence for the correct- 
ness of the foregoing interpretation in the fact 
that when in our high balloon flights we again 
sent up to the very top of the atmosphere in 
place of the pairs of counters (one parallel to 
the other and some six inches above it) either a 
recording electroscope or a single counter, both of 
these two latter receivers being of the kind that 
respond to rays coming in from all directions 
rather than merely to rays coming in vertically, 
then, instead of finding this plateau of constant 
cosmic-ray intensity between Bangalore and 
Agra succeeded by a sudden rise in intensity 
between Agra and Peshawar, we found both the 
electroscope and the single-counter readings 
rising continuously both between Bangalore and 
Agra and also between Agra and Peshawar. 


295 














296 MILLIKAN, NEHER, 


This last behavior is just what is to be expected 
when the electroscope or the single counter is 
the measuring instrument. The reason for this 
difference can be seen from a glance at our dis- 
cussion of the Lemaitre-Vallarta curves repro- 
duced in one of our preceding papers.? 


Il. THE MORE PRECISE MEANING OF 
THE SILICON BAND 


It is desirable to call attention to one further 
characteristic of this band of incoming cosmic 
rays that we have above and, in previous papers, 
designated simply as the silicon-band, namely, 
its possible composite character. Although, ac- 
cording to Bowen’s measurements® Si (at. wt. 
28.3) is probably ten times more abundant in 
interstellar space than are any of the three 
closely adjacent elements Al (at. wt. 27.1), 
P (at. wt. 31), and S (at. wt. 32) these are all 
so close to Si in weight that their annihilation 
rays would be quite inseparable by the type of 
tests we have so far made. 

Hence, for simplicity we designated, and shall 
continue to designate, this observed band of 
incoming rays the ‘‘silicon-band”’ without wishing 
to deny to it the possibility of a composite 
structure. Indeed, P. Swings, who, with his 
collaborators, has recently made elaborate studies 
on the abundance of both atoms and molecules 
in interstellar space,‘ says, “It may safely be 
assumed that the atomic populations in inter- 
stellar space of H, O, N, and C per cm? are 
between 1 and 10-*.”’ These estimates of the 
actual interstellar abundance of these atoms are 
not out of harmony with Bowen’s measurements 
of the relative-interstellar abundance of atoms. 
From these he finds? the He atom 1/10 as 
abundant as the H atom, and the atoms of C, 
N, O, and Si each of about one-tenth the 
abundance of the atom of He. 

Swings also writes us ‘“‘that the interstellar 
abundance of the carbon compound molecules 
CH and CH*+ must be of the order 10-* per 
cm?.”’ In other words, the most abundant inter- 
stellar molecules would have only 1/10,000 of the 

*See Fig. 2, R. A. Millikan, H. V. Neher, and W. H. 
Pickering, Phys. Rev. 61, 401 (1942). 


See R. A. Millikan, H. V. Neher; and W. H. Pickering, 
Phys. Rev. 61, 399 (1942), footnote 5. 

4 P. Swings, Astrophys. J. 95, 270 (1942). See also ibid. 
92, 292 (1942), and Astronom. Soc. Pac. 54, 3 (1942), and 
J. Roy. Astronom. Soc, Canada 35, 75 (1941), 
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abundance of the foregoing abundant atoms C, 
O, N, and Si, so that so far as the origin of the 
“field sensitive cosmic rays’’® is concerned, no 
molecules and no atoms heavier than the atoms 
causing the Si band need to be considered at all. 
This introduces a very great simplification into 
the whole problem of origin. 


Ill. THE 1940 PREDICTIONS OF COSMIC-RAY 
BANDS AND OF PLATEAUS OF CONSTANT 
COSMIC-RAY INTENSITIES 


If, then, our hypothesis as to origin is correct, 
there can be four, and only four, cosmic-ray 
bands in addition to the Si band to look for, 
namely, the O, N, C, and He bands, which 
should have the energies required to enable them 
to begin to break vertically through the resist- 
ance of the earth’s magnetic field at a series of 
four different latitudes as the observer goes north 
from the latitude of first entrance of the silicon 
band. 

However, on account of the dissymmetry of 
the earth’s magnetic field we had computed as 
early as 1939 that the silicon band (energy 
13.2 Bev), though it could not get vertically 
through at the equator in India, where the 
required energy was 17 Bev) could get through 
at the equator off the coast of Peru where we 
computed’ that the requisite energy was but 
13 Bev. 

As soon, then, as we had found the silicon 
band in India we were able to predict, and in 
fact published the prediction, that there should 
be a very long plateau of constant vertically 
incoming cosmic-ray intensity extending from 
Mollendo, near the magnetic equator in Peru, 
clear up to the middle of Mexico (Mag. Lat. 32° 
to 34°) at which latitude the oxygen annihila- 
tion band should begin to come in and of course 
remain in at full strength for all more northerly 
latitudes. 

But the oxygen band like the silicon band 
might have a bit of complexity about it for, 
according to Bowen’s table,’ the Ne atom has 
one-tenth the interstellar abundance of O. Be- 

5 We define non-field sensitive cosmic rays as those found 
at the magnetic equator in India when the measuring 
instrument that is carried to the top of the atmosphere is an 
electroscope. Field sensitive rays are the additional rays 


that come in when the electroscope measurements are 
carried out at successive latitudes between the equator and 


the pole, 
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cause of its weakness and closeness to O (atomic 
numbers 10 and 8, respectively) we had no hope 
of separating the two effects and decided to 
designate the whole effect as the “O annihilation 
band,” as we had done with the similarly com- 
plex silicon band. 

But because of the equality in the interstellar 
abundance of O and N at a latitude a very small 
number of degrees north of the first appearance 
of the O band, the N band should add its equal 
intensity to the total vertically incoming energy. 
We of course had little hope of separating these 
O and N bands, but we had confidence that the 
vertically incoming cosmic-ray intensity at Pasd- 
dena (Mag. Lat. 40.7° N) would represent the 
addition of the intensity of the O and N bands 
and would therefore be notably larger than that 








to be found at the first latitude at which the O 
band came in in its full strength. 

Again, we predicted from our hypothesis that 
in going a very few degrees north of Pasadena 
there would be another increase in vertically 
incoming cosmic-ray intensity roughly equal to 
half the increase between Mag. Lat. 33° N and 
40.7° N, which latter increase should be due to 
the joint effect of the O and N bands, for accord- 
ing to Bowen O, N, and C all had equal 
abundance. 

We expected the most crucial test of this whole 
hypothesis to come in the testing of the prediction 
that between say Mag. Lat. 44° N and 54° N there 
should be found another flat plateau of constant 
vertically incoming cosmic-ray energy since there 
are no abundant atoms at all of weights between 
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those of C and He and we computed from Lemaitre 
and Vallarta’s curves that the He atom-annthilation 
ray could not break vertically through the earth's 
magnetic field farther south than at about Mag. 
Lat. 54° N. 


IV. REASONS FOR MAKING THE 1940 SERIES OF 
INCOMING COSMIC-RAY ENERGY MEASURE- 
MENTS WITH ELECTROSCOPES AT A 
SERIES OF DIFFERENT LATITUDES 


All the foregoing predictions had been made 
upon our return from India in the early spring 
of 1940. We knew, however, that we could not 
develop the necessary new equipment for putting 
these predictions to careful experimental test in 
less than a year’s time. But we were already well 
supplied with excellent recording electroscopes 
such as we had used, in part, in India. Since, 
then, the total cosmic-ray energy coming into 
the atmosphere at any latitude, no matter from 
what diréction, is an important measurement 
which could be made most accurately with 
electroscopes, since, further, a series of such 
measurements taken in a series of latitudes in 
such rapid succession as to reduce the likelihood 
of time-changes was a matter of considerable 
importance, and since this was the only accurate 
way of getting the ratio of the field sensitive to 
the non-field sensitive incoming rays, we decided 
to interrupt, in the summer of 1940 for a period 
of two weeks, the preparations for the above- 
mentioned vertical counter tests planned for the 
year 1941-1942, in order to make a series of 
flights with electroscopes in September, 1940. 

The absorption of all three of us in war projects 
has delayed the working up and presentation of 
these data until the present. 

Meantime, the crucial vertical counter tests of 
the above-mentioned series of predictions as to 
cosmic-ray bands and plateaus were carried out 
in December, 1941 and March and April, 1942, 
and have already been published® with results 
which appear to confirm all of the predictions made. 
However, we have not yet located with the aid 
of vertical counters the latitude of entrance of 
the vertically incident He annihilation rays. We 
hope to do this at the close of the war, but the 
present data, taken with electroscopes, throw 


*R. A. Millikan, H. V. Neher, and W. H. Pickering, 
Phys. Rev. 63, 234-245 (1943). 


AND PICKERING 


some light on this point and on other points of 
importance, as indicated below. 

This expedition, then, was undertaken in the 
summer of 1940 because (1) we had the equip- 
ment ready and could carry out these high 
flights with electroscopes to best advantage in 
the more stable weather conditions usually pre- 
vailing in the Middle West in late summer and 
early fall; (2) because we had already found in 
very high electroscopes flights in high magnetic 
latitudes (Omaha, Bismarck, Saskatoon) as yet 
unexplained variabilities with time’ which were 
conspicuously absent in India® and which we 
thought might disappear in part, at least, if a 
series of flights at different latitudes was taken 
all in one season and in as rapid succession as 
possible; (3) because such a series of more accu- 
rate electroscope flights would establish with 
greater certainty the relation between the field 
sensitive and the non-field sensitive portions of 
the incoming cosmic rays such as had been 
thus far defined only in terms of the total 
integrated value of the ionization found at a 
given latitude through ionization measurements 
made with electroscopes or single counters as 
they are carried from the lowest point at which 
cosmic-ray ionization can be detected up to the 
top of the atmosphere ; (4) because on account of 
the shape of the Lemaitre-Vallarta curves in 
high magnetic latitudes there was a possibility 
of throwing some light, even with electroscope 
flights, on the reality, and latitude of entrance, 
of the helium annihilation band. 


V. THE NEW MEASUREMENTS 


At Bismarck, North Dakota three flights were 
made on the same day, August 21, 1940, the 
first pair of balloons used for carrying up the 
recording electroscope being released at 8:45 
A.M., the second at 9:50 a.M., the third at 
1:46 p.m. Four days later, namely, on August 
25, 1940, three more similar flights were made at 
Omaha, the first ascent starting at 9:37 A.M., the 
second at 10:36 a.m., the third at 11:41 A.M. 
This process was repeated at Oklahoma City on 


™” See R. A. Millikan, H. V. Neher, and W. H. Pickering, 


Phys. Rev. 61, 405-406 (1942). F 
8 See R. A. Millikan, H. V. Neher, and W. H. Pickering, 


Phys. Rev. 61, 408 (1942), and especially Fig. 1 on page 
408. 
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Fic. 2. Soft He rays make Bismarck-Omaha 
difference only at high altitudes; penetrating 
carbon rays cause difference, extending down 
to low altitudes, between Oklahoma City and 
Omaha curves. 


August 29, at Fort Worth on September 2 with 
two flights, and at San Antonio on September 3 
with two flights. All of these flights, then, were 
taken within a period of thirteen days, during 
which the atmospheric conditions were particu- 
larly stable. The set of resulting curves, there- 
fore, has a much better chance of being free 
from errors caused by possible changes with time 
in the incoming cosmic-ray intensities than have 
those we have previously published, some of 
which were taken some years apart. Further- 
more, we have never taken a series of flights 
with different electroscopes that showed through- 
out as great a self-consistency in each given 
locality as do these. 

Figure 1 shows just how we treated the data 
at Omaha, for example. It is altogether typical 
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of the data taken in all five of the different 
latitudes. 

In this figure every one of the readings of the 
photographic film for each of the three flights 
made at Omaha is recorded. Our procedure for 
getting the final mean curve for this locality 
from these readings was to draw separately the 
full curve of each flight; then to read off on each 
of these three curves the ionization at a given 
altitude, say the ionization corresponding to 2 
meters of water beneath the top of the atmos- 
pheré; then to average these three readings at 
2 meters; and having obtained this sort of an 
average at a whole series of chosen altitudes up 
to as near the top of the atmosphere as the 
readings extend, to draw through these averaged 
points the final curve. This is the full line shown 
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in Fig. 1. It is reproduced again in Fig. 2 along 
with the mean curves obtained in this way at 
all of the five indicated latitudes. The area under- 
neath each of these curves is the total integral 
of the cosmic-ray energy entering the atmosphere 
at the latitude in question. It can be reduced to 
electron volts by multiplying by 31, the average 
number of electron-volts required to produce an 
ion in air. 
VI. THE BEARING OF THESE CURVES ON THE 


REALITY, AND LATITUDE OF ENTRANCE, 
OF HELIUM ANNIHILATION RAYS 


From as careful a study as we have been able 
to make of the Lemaitre-Vallarta curves repro- 
duced in Fig. 2 of our earlier paper,’ supple- 
mented by a little adjustment, especially in the 
case of the O band because of our own observa- 
tions shown in Fig. 4, page 403 of reference 9(a), 
and .because of the evidence presented in pages 
239 and 242 of reference 6, that off the west 
coast of South America the electronic energy 
necessary to get vertically through the resistance 
of the earth’s magnetic field at the equator is 
13 Bev instead of the 15 Bev assumed by 
Lemaitre and Vallarta (probably a correct as- 
sumption for computing incoming energies north 
of Pasadena), we have estimated in Table I the 
critical magnetic latitudes for the entry of the 
cosmic. annihilation rays as measured by elec- 
troscopes. 

The table brings to light the bad overlapping 
of the various bands and the impossibility of 
unscrambling them, save in the case of the helium 
band. But because of the wide difference in 
energy between helium and carbon, the former 
is able to appropriate to itself entirely all lati- 
tudes above 51°, so that on the basis of the atom- 
annihilation hypothesis, if in going north new 
cosmic rays come in between 51° and 59°, they 
can only be helium annihilation rays. In fact, 
as shown in Fig. 2 there can be no uncertainty 
about the appearance of a new band of rays 
between Omaha and Bismarck and a band so 
little penetrating, as helium rays should be, that 
it cannot throw its influence down even to the 
2-meter level, after which depth the Bismarck 


*(a) Millikan, Neher, and Pickering, Phys. Rev. 61, 397 


(1942). (b) See also Lemaitre and Vallarta’s original papers 
noch te Rev. 49, 719 and 50, 503 (1936); also 43, 87 
1 . 


curve is seen in Fig. 2 to coincide completely 
with the Omaha curve. 

This is not at all so with the difference between 
the Omaha and the Oklahoma City curve be- 
cause this difference is primarily due to the 
carbon band, for carbon rays have an energy 
of 5.6 Bev and this is sufficient to throw their 
influence down nearly to sea level. If, on the 
other hand, one assumed incoming cosmic-ray 
energy lying between the C and He bands, then 
the Oklahoma City curve would certainly be 
expected to run together with the Omaha curve 
at a much higher altitude than Fig. 2 shows to 
be the case. 

We expect as soon as possible after the war to 
locate the latitude of entrance of the vertically 
incoming helium rays by the vertical counter 
technique, but the evidence here as to its exist- 
ence and location is so good that we feel con- 
fident that this prediction will be verified as well 
as this has been done with the whole series of 
predictions discussed above. 


Vil. THE VARIABILITY IN INCOMING 
COSMIC-RAY INTENSITIES 


We have published on several occasions in 
recent years the discovery of large and thus far 
wholly inexplicable changes in incoming cosmic- 
ray intensities as measured by the electroscope 


TaBLE I. Critical magnetic latitudes in North America for 
entry of cosmic annihilation rays. 








Latitude of Latitude of 








first entry full entry Latitude of 
on western on eastern vertical 
Atom Energy horizon horizon entry 
He 1.88 Bev ~51° ~59° ~54° 
> 5.6 Bev ~32° ~49° ~42° 
N 6.6 Bev ~27° ~48° ~39.5° 
O 7.5 Bev ~22° ~46° ~33.5° 








technique in the northern part of the United 
States.!° 

Thus, between August 21 and September 3 of 
1937 we made good flights with electroscopes at 
Omaha to such heights as to get definitely over 
the maximum of ionization in our electroscopes. 
The value of this ionization at the top when 
reduced to standard air density was 338 ions 


See R. A. Millikan and H. V. Neher, Phys. Rev. 56, 
491 (1939). 
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per cm’. In 1938 we went purposely to Omaha 
in winter (December 22 and 23) and repeated 
these experiments and found a maximum ioniza- 
tion at the top of 364 ions, an increase of nearly 
8 percent. Again, the maximum ionization shown 


TaBLE II. Time changes in the earth’s magnetic field. 
H in c.g.s. units. 











1915 Annual changes 1940 
Bismarck 0.161 — 0.00027 0.154 
Omaha .197 — .00049 .185 
Oklahoma City 241 — .00063 225 
St. George .247 — .00037 .238 
Mt. Wilson .265 — .00035 .256 


Magnetic records from Mount Wilson, California; 
Tucson, Arizona; Cheltenham, Maryland 
1937 September 3 Nostorm 
1938 July 10 Nostorm 
July 15 Small storm July 13 to 17 
December 22 H about 0.00030 lower than normal 
1940 August, September No storms. 








at the top in the September 25, 1940 flights 
herewith reported (see Figs. 1 and 2) is 413 ions 
per cm*, an increase over December, 1938 of 
13.4 percent. 

Again, the maximum which we obtained in 
Bismarck in the summer of 1938 (July 5) was 
374 ions/cm*, while the 1940 maximum found at 
Bismarck August 21 and reported in Fig. 2 was 
485 ions/cm’, or an increase over 1938 of 29.7 
percent. 

Still further, our maximum obtained in Okla- 
homa City in mid July (11, 12, 13) of 1938 was 
280 ions/cm*, while the maximum shown in 
summer of 1940, August 29 (see Fig. 2), is 319 
ions/cm’, an increase over 1938 of 14 percent. 

Contrast the foregoing findings with the follow- 
ing facts. We made electroscope flights at San 
Antonio in the summer of 1936 and got a maxi- 
mum of 234 ions/cm* and repeated these flights 
on September 3, 1940 and found (see Fig. 2) 
precisely the same maximum as in 1936, viz., 
234 ions/cm’. Also, not only was our maximum 


at Madras, India in 1936 the same as in Febru- 


ary in 1940, but as shown in the report of our 
India work" the two curves taken four years 
apart are everywhere indistinguishable. 

In a word, while the cosmic rays coming in at 
Madras and at San Antonio seem to remain 


amazingly constant, yet in the northern part 


See H. V. Neher and W. H. Pickering, Fig. 1, Phys. 
Rev. 61, 408 (1942). 


of the United States they show amazing fluc- 
tuations. 

We seek the explanation of this strange be- 
havior in the atom-annihilation hypothesis as 
follows: We assume for simplicity that the cosmic 
processes creating the cosmic rays cause an es- 
sentially constant incidence of all the five cosmic- 
ray bands on the earth. According to Bowen, He 
is ten times more abundant than C, N, O, or Si, 
and even from the energy standpoint, if the 
probability of the transformation of the whole 
rest mass into an electron pair is the same for all 
atoms, the energy in the helium band is some 3.5 
times that in the C, N, and O bands, and some 13 
times the energy of the Si band. 

Bismarck and Omaha both lie close to the 
southern edge of the polar cap of He annihilation 
rays, which we have computed as located for 
vertical-counter measurements at 54° N Mag. 
Lat. and for electroscope measurements at be- 
tween 51° N Mag. Lat. and 59° N Mag. Lat. With 
even slight changes in the earth’s magnetism this 
edge of the He polar cap will move north or south 
and, since the resistance to incoming electrons is 
very small in northern latitudes, the relative 
effect of changes in the earth’s magnetic field may 
be correspondingly large depending upon the kind 
of cause that makes the change. In equatorial 
latitudes, where the resistance to incoming elec- 
trons is say 8 times larger than at Bismarck or 
Omaha, the fluctuations in the earth’s magnetic 
field will have smaller influence and no influence 

TaBLeE III. Tucson values of H corresponding to the 


year and month at which the foregoing cosmic-ray measure- 
ments were made. 





—_ 


Omaha Bismarck 





Oklahoma City 





Sept. 1-30, 1937, 261897 July 1-31, 1938, 261827 
Dec. 1-31, 1938, 261627 Aug. 1-31, 1940, 261417 
Aug. 1-31, 1940, 261407 


July 1-31, 1938, 26182y 
Aug. 1-31, 1940, 261417 








at all unless the point of observation is near the edge 
of one of the five polar caps corresponding to the 
latitude of entrance of one of the cosmic-ray bands. 

In the present experiments according to Table I 
the helium rays from the western horizon are 
already passing through the electroscope in the 
upper levels of the atmosphere and a very little 
change in the earth’s magnetic field that would 
make the band move say 2° southward, according 
to Lemaitre-Vallarta curves, would throw the 
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greater part of the helium band’s energy into the 
top of the Omaha curve. It is presumably actually 
there 1n these observations in which the Omaha 
curve runs higher than we have observed it to do 
in any of our other trials there. 

Oklahoma City is the only place, save Bismarck 
and Omaha, at which these changes have been 
observed, and here they were in smaller amount. 
Oklahoma City is actually quite close to the 
computed latitude of entrance of the carbon 
band. The reason for the high values of the 
Bismarck, Omaha, and Oklahoma City curves in 
1940 is that then the edge of the He and C polar 
caps were farther south than at the time of previ- 
ous readings. Furthermore, as Table I shows, the 
width of the slope of the rise to a new plateau 
increases greatly for the more southerly bands. 
This is another reason why a small change in the 
earth’s magnetic field has little effect in the more 
southerly latitudes, but a large effect in the 
northerly ones. 

According to this hypothesis, an increase in the 
intensity of the cosmic rays at Bismarck and 
Omaha means a southward movement of the edge 
of the helium polar cap and such a movement in 
turn means a weakening in the strength of the 
earth’s magnetic field. Such a weakening is 
actually shown in Table II which has just been 
prepared by Dr. S. B. Nicholson of the Mount 
Wilson Observatory who is looking for correla- 
tions between sun-spot activity and magnetic 
storms, and who also tells us that magnetic 
storms in general correspond to a weakening 
rather than a strengthening in the earth’s field. 
Further evidence of such weakening of the earth’s 
field between May, 1936 and May, 1938 can be 
seen in Fig. 1,!* which we took from the paper of 
S. E. Forbush," of the Carnegie Institution of 
Washington. 


12 See R. A. Millikan and H. V. Neher, Phys. Rev. 56, 


492 (1939). 
13S, E. Forbush, Phys. Rev. 54, 975 (1938). 


In further check on the foregoing apparent 
correspondence between cosmic-ray intensity and 
the strength of the earth’s magnetic field in the 
western part of the United States, Director L. O. 
Colbert of the Coast Geodetic Survey has kindly 
furnished us with the data which we have put 
together in the form shown in Table III. This is 
more definite than Table II in indicating that at 
every one of the above observed times (year and 
month) in whith we observed an increase in 
cosmic-ray intensity over its value at the same 
station at some other year and month there had 
also occurred a corresponding small decrease in 
the horizontal component H of the earth’s mag- 
netic field at the Tucson magnetic observatory. 


VIll. THE FIELD SENSITIVE COMPONENT 
OF THE COSMIC RAYS 


‘In view of the greater accuracy in these 
electroscope measurements than of our preceding 
ones, we have computed from the new curves of 
Fig. 2 combined with the unchanged measure- 
ments at Madras," that the field sensitive com- 
ponent defined through these electroscope readings is 
65 percent of all incoming cosmic-ray energy. All 
this is due to incoming charged particles, cer- 
tainly electrons or protons. There is at present no 
way of finding what fraction of the non-field 
sensitive component is due to charged particles 
and what to photons, though part of it is 
unquestionably due to the former. 

This investigation, like many preceding ones, 
has had support from the Carnegie Corporation 
of New York administered by the Carnegie 
Institution of Washington. For this aid, as well as 
for the aid given us by Dr. Seth B. Nicholson and 
Director Colbert, the authors wish to express 
their keen appreciation. 


4 See Figs. 1, 7, and 8 of Millikan, Neher, and Pickering, 
Phys. Rev. 61, 397 (1942). 
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The paper contains the results of calculations of (1) repulsive exchange forces, (2) attractive 
van der Waals forces between H? and H atoms. The method used is that developed in an earlier 
communication. Results for small distances of separation differ considerably from published 
data on interaction energies deduced from collision cross sections. Reasons for the discrepancy 
are discussed. Figures 1 and 2 represent potential energy curves for the region around the 
van der Waals minimum; they exhibit the asymmetry of the interaction. 





I. THEORY 


HE forces between Hz molecules have been 
calculated with the use of a simple molec- 

ular wave function in a recent publication.’? The 
method and many of the details there developed 
lend themselves readily to the calculation of the 
interaction between an atom and a molecule, a 
matter regarding which there is some experi- 
mental information.* In the present paper we 
present the significant points of the theory and 
the results, making use throughout of the nota- 
tion of reference 1. The simplifications formerly 
made will again be used ; in particular, a Heitler- 
London function with adjustable screening con- 
stant is chosen to represent the molecule. The 
atom will in this case be described by a function 
having the same screening constant as the mole- 
cule, a procedure which reduces the complexity 
of computations greatly without, it seems, fal- 
sifying the results in an essential way. To show 


. the effect of the somewhat arbitrary choice of 


screening constant, the calculations are carried 
out for several values of this parameter. The 
preceding paper! will hereafter be referred to as I. 

Let the nuclei of the three interacting atoms 
be localized at the points a, b, and c. These letters 
will also be used to denote the electronic orbital 


1H. Margenau, Phys. Rev. 63, 131 (1943). 

? Meanwhile, an interesting paper by J. DeBoer, Physica 
9, 363 (1942), dealing with the same problem, has come to 
my attention. His final results, which do not differ greatly 
from ours, are obtained more simply and with much less 
labor by neglecting most of the exchange integrals, to the 
investigation of which reference 1 was chiefly devoted. As 
DeBoer admits and as shown in reference 1, these are 
quite as large as the terms retained in his calculation. It 
appears that because of this difference in treatment, 
DeBoer obtains agreement with observation for Z = 1.166, 
whereas we were forced to assume a greater value. 

+I, Amdur, J. Chem. Phys. 11, 157 (1943). 


functions centered about these points. Thus, in 
the notation of I, three state functions describing 
the three-particle system can be formed for which 
the Z component of the total spin is zero: 


Vi=|aa bB cal, 
W.,=|aB ba cal, (1) 
VW;=|aa ba cB. 


Of these, we can construct by known rules two 
doublet functions: 


Wa =/}(Vi— V2), va=V/3(Vs—V2). 


In general, the correct VY is a combination of 
these two, and the energy E is given by 


Haa—EAsa Hap—EAap 
=(0) 
|\Has—EAazp Happ—EAgpep 


We note, however, that Ya represents a bond 
between atoms a and b, wg a bond between a 
and c. (The criterion for existence of a bond 
between two atoms is symmetry with respect to 
an interchange of the bonded nuclei.) If, then, 
we suppose that a and } form our molecule, we 
may limit our consideration to Wa alone. In this 
case, the determinantal equation reduces to 


Haa 


Aaa 





(2) 


The functions a, b, c used in Eqs. (1) are all 
one-electron orbitals and hence include no ionic 
terms. At close distances of separation ionic 
states certainly make a contribution to the 
energy and, in an exact calculation, have the 
effect of lowering it. Even when the functions 
describing the atom and the molecule at infinite 
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separation are not exact, as is the case in the 
present instance, a similar effect is to be ex- 
pected. Neglect of ionic states is one of the more 
serious approximations made in these cal- 


culations. 
In the notation of I, 


Haa=}(Au+Ae2—2Aw), 
Aaa =}(A1r +422 — 2A12). 


When the Hamiltonian function for the three- 
particle system is written down and use is made 
of the Schrédinger equations satisfied by a, 3, 
and c¢, there results 


H{a(1)b(2)c(3)} = {3D—(1—Z)(ai+62+7s) 
— (agtast+hitBetrvi+72) 
+ pret pist+ post Ew }a(1)b(2)c(3) 


where again D= —}Z*e*/ao, and Ey is the elec- 
trostatic energy of the three nuclei. If, in analogy 
with I, we define 


e=(abc|H|abc), ui:=(bac|H|abc), 
us=(cha|H\abc), us=(ach| H|abc), 
v, = (bca| H| abc) 


and denote the corresponding matrix elements of 
unity by primed symbols, we find 


(3) 


Hy=e—-u, Hy=—m4+01, Ho.=€— U4, 
and hence by (3) 
Haa=e+u1—}3(U2t+u4) —21, 


Aaa = +n)’ —_ 4(u2’ +4’) —v;'. 


(4) 


The exchange integrals now have the following 
form: 


€=3D+3xKA+CatCact Cre, 
u; = A?.~4{3D+Xat Bact Boe 
+«(A+2T)+X'a}, 
where X’.. now has the simpler form 
(acpbc) — (ayb) 
2 m 


Uue=u, with interchange of 5 and c 


, (5) 





u4= 4, with interchange of a and c 
01 = AgAachse{ 3 D+ (x — 1)(Tas+ Tact Tre) + Vase} ; 
xk=Z-—1. 


Henceforth quantities bearing subscripts a4, i.e., 
pertaining to the molecule, will be written without 
indices. When atom and molecule are very far 
apart, the exchange integrals above reduce to 


e=3D+C+3x«A, e =1, 
u,=A?{3D+X+K(A+2T), mu’ =A’, 
Ue = U4=0,=0, Ue’ =u’ =v,’ =0. 
On using these in (4) and (2) there results 
C+2xA +A?(X +2xT) 
1+? . 





E°=3D+«A+ (6) 


This, however, is the energy of a single molecule 
[cf. I, Eq. (15) ] plus D+«A, an increment which 
represents the energy of one hydrogen atom 
computed with the function a, that is, D+«A 
= faHatomadr. Hence, in order to find the ex- 
change energy between the molecule and the 
atom, it is necessary to subtract E® given by (6) 
from E computed in accordance with Eq. (2). 
In proceeding further, prohibitive labor is 
encountered unless use is made of the reduction 
described in I. The approximations called for are 
different when the distance of approach is small 


TABLE I. Exchange energies for small distances of 
separation. Atoms a and } form a molecule. Ra. and R,, 
are the distance of atom c from a and 3, respectively. 
R is the distance of c from center of molecule. 











Z Rae(A) Robe(A) R(A) E —E® (ev) 
a 1.071 0.57 0.17 0.20 104 
b 1.071 0.42 0.42 0.20 72 
c 1.071 0.84 0.10 0.47 117 
d 1.071 0.60 0.60 0.47 26 
e 1.20 0.83 0.09 0.46 142 
f 1:20 0.58 058 0.45 29 
g 1.30 0.82 0.08 0.45 150 
h 1.30 0.45 0.45 0.26 56 








and when it is large. We consider the case of very 
small separations first. 


Il. SMALL DISTANCES OF SEPARATION 


Let atoms b and ¢ come very close together. It 
is then necessary to retain all terms in the 
exchange integrals, but the slowly varying func- 
tions bearing indices 6 and ¢ may be replaced by 
their limiting forms in which the distance 
between atoms 0 and ¢ is put equal to zero. In 
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Fic. 1. Interaction between He and H, calculated for 
Z=1.071. Curve a represents approach of H along bisector 
of molecular axis, curve } along axis. 


that case 
Xa—Xec—Mes—2T ow, 
Xt-—2(abpab) — 2(aBa), 
Uste>X ab +See +See — (aBa) + Mos. 


The only function appearing here which was not 
explicitly given in I is 
¥o(2) 


3 er 
M=2T-- 
41+S+(1 /3)S2 

As before, S is the distance indicated by the 
subscripts, measured in units a/Z. 

Computations of E were made for 3 different 
values of Z and for certain positions of atom c 
relative to the molecule which were selected for 
numerical convenience. The results are collected 
in Table I. It will be noted that cases a, ¢, e, 
and g represent positions of atom c on the axis 
of the molecule, the others on the bisector of the 
axis. In the former group, the atom is close to 
one of the molecular nuclei; hence the inter- 
action energy is very large. For these small 
separations, the distance R between atom and 
center of the molecule is not a significant 
parameter. As will appear later, the best value 
of Z for the present problem is probably near the 
first of those here chosen. The results indicate, 
however, that the choice of Z is not very impor- 
tant at these small distances of separation. 

Comparison with the results of Amdur? yields 
very poor agreement indeed. If we limit our con- 
sideration to cases where the atom is situated on 
the bisector of the molecular axis—this being the 








more probable condition of impact—instances 
d and f should be comparable with the data 
given by Amdur, who finds approximately 
8 and 8.5 ev instead of our 26 and 29 ev. This 
discrepancy by a factor greater than 3 can hardly 
be accounted for in its entirety by inaccuracies 
in the present calculations, or by our omission 
of ionic terms in the state function. It seems 
indeed that the assumption of spherical sym- 
metry of the interaction, made in the deduction 
of the experimental curve from scattering data, 
must be at fault. The discrepancy is materially 
reduced if the collision radius in reference. 3 is 
regarded as the distance of closest approach, not 
to the center of the molecule, but to the nucleus 
of one of the molecular atoms. 


Ill. LARGE DISTANCES OF SEPARATION 


For large distances of separation, a different 
set of approximations must be made in evaluat- 
ing E. Numerical computation shows that the 
quantity X’,, appearing in (5) may be neglected. 
Furthermore, we split Ha and Axa into their 


ANT | 
| 


mai : 
FE 
NY 


T 


a 
4 














Milhvoits 








ie} 
7 














a 



































Fic. 2. Interaction between Hz and H, calculated for 
Z=1.20. Curve a represents approach of H along bisector 
of molecular axis, curve } along axis. 


large and small constituents: 


(0) (2) (0) (2) 
Haa=HaatHaa, Asa=AaatAaa, 


where the large parts H® and A® are those 
which were retained in arriving at Eq. (6). Then 
(0) (2) co), (2) 
rat 44 Has, Has” Haa ——J 


Aaa (0) (0) (0) 2 
Aaa Aaa (Aaa) 





if squares of the small terms are neglected. The 
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TABLE II. Exchange energy E=E® in units Ze*/ao. S is 
the distance between the atom and the center of the 
molecule in units ao/Z. In position a, the atom lies on the 
bisector of the molecular axis; in position 5, on the axis. 








Z=1.071 Z=1.2 Z=13 
Pos. a Pos. b Pos. a Pos. b Pos. @ Pos. b 


0.(2)1607 0.(2)192 0.(2)332 —-0.(2)442 Hest 0.(2)644 
0.(3)331 =—-0.(3)412 0.(3)760 = =—0.(2)1070 =: 0.(2)1119 ~—-0.(2) 1648 
0.(4)638' = 0.(4)748 0.(3)1695  0.(3)238 0.(3)243 .(3)383 
0.(4)1015 = 0.(4)1277 | :0.(4)323 = «(0.(4) 495 0.(4)486 0.(4)822 
0.(5)1970 = 0.(5)2141 0.(5)637 —-0.(5) 947 0.(5)983 0.(4)1596 
0.(6)358  —0.(6)379 0.(5)1261 —0.(5)183 0.(5)2001 0.(5)334 


i) 





_ 
SOMA H 








first of the three fractions on the right of this 
equation is E°®, hence 


0 (0), —1, (2) 0 (2) 
E-E = (A4a) (Has—E Aaa). 
When all terms on the right are inserted and the 
exact expression for U,s- is replaced by its 
limiting form for Ra.—0: 


e 
Uste=—+X(S)+2T(5), 


the result is 
E— E®=(1+A*)™{ Cact+ Coc +A?(Bact+ Boe) 
+4A,.[4Eo’ —B—Xoc—2kT ec] 
+$A3.[4Eo! —B—Xve—2xTre] 
+AAgeAs-[ $ Eo’ —s~!—(«—1)T 
+«(A —Toc—Tre—X(S))]}. (7) 


The notation is explained in I. If all lengths in 
this expression are taken in units a)/Z, the energy 
is in units Ze?/dp. 

The choice of Z calls for comment. It would 
seem at first sight that some mean value between 
that proper for the molecule and that for the 
atom should be used. The minimum of the 
molecular energy occurs for Z = 1.16; the correct 
value for the atom is of course 1. In some of our 
calculations we have taken Z to be 1.071 (which 
makes s=1.5). In our former work (I) on the 
four-electron problem there appeared a peculiar 
difficulty: The value of Z corresponding to the 
minimum of molecular energy produced an at- 
tractive exchange force. On analyzing the situ- 
ation we found reason to accept a larger value 


of Z, the chief argument involving the greater 
repulsive effect between the electrons. A similar 
difficulty does not appear in the three-electron 
problem ; the exchange force is repulsive for all 
Z that were tried. Nevertheless, if the former 
argument is valid, a value in the neighborhood 
of 1.2 should be a better choice. Hence calcula- 
tions were made for that case. The difference 
between the results for Z=1.071 and 1.2 is 
surprisingly small. To complete the picture we 
have also studied Z=1.3. But the minimum re- 
sulting from the composition of exchange and 
van der Waals forces is much deeper in that case 
and may be dismissed from consideration. 

Results for E—E® are given in Table II. In 
position a, the atom is situated on the bisector 
of the molecular axis, in } on the axis. 

There remains now the calculation of the long 
range attractive forces. The method is outlined 
in I. On using for the atom an oscillator strength 
f=1 and a corresponding resonance energy 
0.487 e?/ao (which gives the correct coefficient of 
R-* in the interaction between H atoms), we 
obtain for the dipole-dipole energy 

ao 6 ¢? 
E.:= -8.40(— om, (8) 
R ao 


The dipole-quadrupole energy, calculated from a 
formula previously given,‘ becomes 
ao 8 ¢ 


E24= -148(— am, (9) 
R ao 


When (7), (8), and (9) are finally compounded, 
there results the total interaction which is 
plotted, for 2 values of Z, in Figs. 1 and 2. There 
are at present no empirical data with which they 
may be compared. The minimum comes at about 
the same place as it does for the two similar 
positions in the bimolecular problem, but it is 
slightly less deep. The ‘“‘shape”’ of the molecule 
again manifests itself in the differences between 
curves a and b. 


4H. Margenau, J. Chem. Phys. 6, 896 (1938). 
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‘mer 
100d In previous work on this subject, the second virial coefficient has been calculated with 
ula- numerical success, but on the assumption of intermolecular forces which take inadequate 
ence cognizance of the known structure of the water molecule. In the present paper, all long range 
2 is forces are computed in a semi-empirical way from dipole moment, quadrupole moment in 
+ in accordance with the Bernal-Fowler analysis, and optical dispersion data. Interactions caused 
; by quadrupole moments cannot be neglected. Short range repulsive forces are thus left unde- 
1 re- termined. These are chosen in such a way that the second virial coefficient, in its dependence 
and on temperature, approximates its measured values as closely as possible. It is found that a 
case reasonable fit cannot be obtained with the use of the hard sphere model, nor with any model 
which permits the polar component of the force to be effective at arbitrarily small distances of 
i separation. More satisfactory agreement results if the polar forces are eliminated at all dis- 
| tances smaller than the potential minimum. The best potential satisfying these various require- 
ctor ments is given at the end of the paper [formula (22)]. 
long 
ined - 
ath INTRODUCTION so-called induction effect. (3) The forces arising 
rey HE present paper has as its principal from a mechanism similar to the preceding, but 
t of purpose the determination of the forces through fields caused by the rapidly rotating 
we acting between water molecules. To this end, ™0ments associated with electronic motions 
we rely on the theory of van der Waals forces to within a molecule instead of permanent moments. 
provide an expression for the interaction at large This is known as the dispersion effect. It: is 
(8) distances of separation. The short range repulsive Present regardless of whether the molecule has 
forces are at present beyond purely theoretical @ multipole or not. 
treatment. It might be expected, however, that The short range forces are repulsive; their 
ma when the attractive forces are known, available 0O™8'™ lies partly in the electrostatic repulsion 
empirical data, such as the values of the second between the molecular nuclei, partly in electron 
virial coefficient at different temperatures, would changes. They are not of spherical symmetry, 
(9) fix within reasonable limits the repulsive forces but their angular dependence has not been 
as well. This will prove to be the case. described for the case of water molecules. 

To give a general survey, we list first the We now give a brief survey of the history of 
led, various types offforces which compose the total the present problem. The theory describing 
Be. interaction and” state their physical origins. effect (1) of the preceding : paragraph was de- 
acre Among the long range forces there are: (1) The veloped by Keesom' at the time when all a ow 
hey attractions or repulsions between the permanent Waals forces tar thought 3 ee = 
jout multipoles borne by a water molecule. These ™#"ent multipoles. His treatment of dipole 
ilar include not only the customary dipole forces to interactions is adequate for the problem at hand; 
t is which the attention of most writers has been his work on quadrupoles, which considers only 
cule limited, but also quadrupole and higher types. linear structures, is in need of generalization 
een For it will be seen that H,O carries a sizable >efore application to water molecules can be 

quadrupole moment whose presence may not be made. Furthermore, Keesom considers molecules 
neglected. Since the Boltzmann factor makes having only a dipole, or only a quadrupole 
alignment of the multipoles more probable than renee" ‘The exact nen for the coins 
other orientations, the net effect, called align- tion of both is very eawinay ~ . = . =" 
ment effect, is attractive. (2) Attractions caused a here; an approximation will sulice 
by the induction of additional multipole moments palhctsinitgeacane 

by the fields of permanent moments—the 1 W. H. Keesom, Physik. Zeits. 22, 129 (1921). 
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The induction effect was considered by Debye? 
and Falkenhagen.* Formulas developed by them 
may be used in this work. For the theory of dis- 
persion forces introduced by London, we refer 
the reader to a review article by one of the 
authors.‘ 

An early analysis of the H,O problem® which 
included consideration of all effects here men- 
tioned resulted in the prediction that the 
molecule, in order to produce virial coefficients 
of the known magnitudes, possesses a quadrupole 
moment of about 5X10-** e.s.u. In this work, 
the repulsive forces were treated somewhat 
inadequately, their effect being described by an 
infinite repulsion at a distance of separation 
equal to the kinetic theory diameter. (Hard 
sphere model.) This fiction, frequently employed 
in the past, makes the results extremely sensitive 
to the choice of that diameter, which is empiri- 
cally ill defined and therefore sometimes leads to 
a specious semblance of success. Wohl*® and 
Briegleb,’? using similar methods and approxi- 
mations, discussed the situation from a general 
chemical point of view. 

More recently, the problem has been recon- 
sidered by Stockmayer.* He shows that the 
temperature behavior of the second virial coef- 
ficients of water and ammonia is described with 
numerical correctness by assuming the molecules 
to be point dipoles. His calculation permits much 
latitude in the choice of other parameters of the 
interaction and shows, in fact, that the forces 
producing such agreement are far from unique. 
In a sense, this type of success is satisfying, but 
it teaches very little about the actual properties 
of the molecule and augments our knowledge in 
no fundamental way. Since our prime concern is 
with the actual forces, the assumption of point 
dipoles is inadequate. 

The work of Goff, Anderson, and Gratch® 
contains the significant indication that the hard 
sphere model, when used in conjunction with a 


2 P. Debye, Physik. Zeits. 21, 178 (1920). 

3H. Falkenhagen, Physik. Zeits. Sowjetunion 23, 87 
(1922). For a correction to this paper, see reference 5. 

4H. Margenau, Rev. Mod. Phys. 11, 1 (1939). 

5H. Margenau, Phys. Rev. 36, 1782 (1930). 

* K, Wohl, Zeits. f. physik. Chemie B23, 105 (1933). 

7G. Briegleb, Zeits. f. physik. Chemie B14, 36 (1931). 

8 W. H. Stockmayer, J. Chem. Phys. 9, 398 (1941). 

*J. A. Goff, J. R. Anderson, and S. Gratch, Heating, 
Piping and Air Conditioning 15, 319 (1943). 
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plausible set of force parameters, is not capable 
of yielding the correct temperature dependence 
of the second virial coefficient for water. This is 
true even when the essentially spherical force 
field, employed as an approximation by these 
authors, is replaced by a field of dipole symmetry 
as will later be shown in this paper. 

We consider the various effects of interest in 
the order in which they were listed in the 
second paragraph of this introduction. 


INTERACTION BETWEEN PERMANENT 
MULTIPOLES 


For the sake of completeness and systematic 
development, we include in this section a few 
results that are known. The method is perhaps 
more direct than those found in the earlier 
literature. 

Let a set of discrete point charges e; be situated 
at vector distances r;= (xj, yi, 2;) from an origin 
0. Another set of charges « is placed at distances 
or. = (&, m, pr) from an origin 2; the displacement 
of 2 relative to 0 is R. The mutual static potential 
energy between these sets is then 


C5€ 


7 : 
a |R—r+o| 





If we introduce “relative displacement vectors”’ 
0, =1:i—, V can be expressed in two equivalent 
simple forms, the operator form 


| 1 
V=< ee exp (-va-v)(—), (1) 


which results from a Taylor development, and 
the Legendre form 


€:€ ga\' 
V= s R u P (cos ain) (=) ’ (2) 


obtainable on application of a well-known ex- 
pansion theorem in terms of Legendre poly- 
nomials. In the last expression, aa is the angle 
between ea and R: cos aa=(2:— {/on) provided 
we take the Z axis along R, as will henceforth be 
understood. Equation (2) is the more useful for 
our present purposes. If 


6: = La = 0, 
» 
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as we are supposing, the terms for which ]/=0 
and 1 make no contribution to (2). For suffi- 
ciently compact charge distributions and suffi- 
ciently large R, this series converges rapidly. In 
most molecular problems, terms beyond /=4 
need not be considered. To this approximation, 
insertion of the Legendre polynomials gives 


V=-—- (1/R*)>> €,€\( 22:6, — xika —yim) 
‘ 
r 
+(3/2R*) VeweLr?f,—ipr? 


+ (2x:& + 2yim — 32:fn) (81 — fr) J 
+ (3/4R®) Seve(r i? py? — 522 py? — 572g? 
f 


» 


- 152,°5.?+ 2(42:6: —xit. —yim)? ] 
+(1/2R®) ¥ee[3(r2+ pr?) 


d 


X (42if, —xit — yim) 
+5(2°+ 67) (3xi&+3yim —42ih) J. (3) 


Summations 1 to 4 represent, respectively, 
dipole-dipole, dipole-quadrupole, quadrupole- 
quadrupole, dipole-octupole interactions. 

In physical problems, it is usually necessary 
to calculate 


fee (— V/RT) sin 0,40, sin 02d 60d ¢, 


- f exp(—V/kT)d9, (4) 


6, and @ being the polar angles between fixed 
lines in the molecules and Z, ¢ the relative 
azimuth. If only the first summation of V is 
retained, this integral can be evaluated exactly, 
for V then has the simple form 


— (pip2/R*)(2 cos 0; cos A 
—sin 0, sin @,cos g), (5) 


where the ~’s are dipole moments. But in general, 
a series expansion of exp(—V/kT) must be 
performed before integration. The number of 
terms to be retained in this series depends on T. 
In the application to be made later, one should 
strictly use terms up to V?°, but the labor in- 
volved in such a computation appears pro- 
hibitive. We therefore adopt the following 
approximate procedure. 


Since §V"dQ vanishes for odd n, the first 
term which occurs in the expansion is § V°dQ. 
This will be computed exactly for the whole of 
expression (3). For sufficiently high temperatures 
this is the decisive contribution to (4), though 
that is not true in our application. We then 
evaluate (4) including the dipole term only, but 
augmenting the coefficient p:)2/R* in (5) in the 
proper way so that the mean square of the ap- 
proximate V is nearly equal to { V*dQ. This will 
require the addition of a term in R- to pipo/R, 
a term characteristic of the quadrupole inter- 
action. We proceed now to the calculation of 
S V°dQ=8r( V?)w. 

In squaring V, the cross terms between dif- 
ferent summations of Eq. (3) may at once be 
omitted, for they will be seen to vanish when the 
integration over angles is performed. A great 
number of others may also be omitted; the jus- 
tification for this curtailment will appear later 
on. The significant terms to be retained are given 
here in part: 


V? = (1/R®) dese eneu (42:2 50r0, 
+x 5x jf,E, +99 5m) 
+ (9/4R*) Veen lrer Pit, 
+2;:25pr*p,?+4x 25x 2jh,5,+---) 
+(9/16R")D)---+(1/4R")D +--+. (6) 


To obtain ( V*), it is thus necessary to compute 
averages of the type (x;"xj")w, (xi"Vj")w, etc. This 
procedure is somewhat complicated by the fact 
that the vectors r; and r; are not independent— 
since the whole molecule is being orierited in the 
evaluation of the averages in question—but lie 
along rays always making a fixed angle 6;; at 0. 
However, averages over Greek and Latin letters 
may be carried out independently because they 
refer to different molecules. 

A convenient method for carrying out the 
averaging process with respect to connected 
coordinates appears to be this. We introduce 
Euler angles a, 8, y and let R(a, 6, y) be the 
matrix which, when applied to the vector r, 
produces rotation about the origin through 
a, 8, y. Starting with any special, permitted 


position of the vectors r;’° and rz, such as 
r:° = rik, r.° = ro(sin 612j-+cos 62k), 


10 We now write 1, 2 for the former indices 7, . 
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we obtain their most general position by applying to them the matrix R: 


= Rr,’, 


These vectors have components 


t= Rr,°. 


x= —r, sin B cos y| X2=7e sin 62.(sin a cos 8 cos y+cos a@ sin y) —72 COS 632 sin 8 cos 


v1 =r, sin 8 sin y 


2:=7, cos 8 


Ye=Pe Sin 642(—sin a cos B sin y+Ccos a cos y)+72 COS 612 sin B sin 


Z2=P2 Sin Oy sin a sin B+72 cos 62 cos B. 


The desired average over any quantity Q is then simply 


1 2r 5 2r 
n=— f da f sin ads f dyQ, 
8x? Jo 0 0 


with 6,2 held fixed and Q expressed through the preceding relations. The results needed here are: 


(204% 5) av = (1/3)rir; cos 63; 
(1 j)w = 0 


(x43x5)m = (1/5)r.87; cos 04; ; 
(42x 7) = (1/15)r3?r 7(14+2 cos? 6;;) 
(xs*¥ 3) = 0 


(x 52:52 537) wv == () 


(x4%y 3?) = (1/15) 27 7(1+sin? 6;;) 


(0579 7?) = (1/15) rv? cos 65; 


In addition to this list we note that averages 
over terms of odd power, such as x*y, are always 
zero. Many other relations may be produced by 
noting two operations with respect to which the 
formulas above are invariant: (1) interchange of 
any two coordinates; (2) interchange of the 
indices and j throughout. If averages over 
higher powers of V than the second were needed, 
similar integrations over a greater number of 
connected coordinates would have to be per- 
formed. 

Combination of (6) and (7) yields finally a 
relatively simple result for (V?)w: 


2 
(V?) av = D616 7€r€, a" 3P1(cos 843) prpyPi(cos A,) 
ij 
Au 


1 
te lreriPs (cos 653) Pr pu (cos Ay) 


+17 sP1(cos 0:3) pr? p,2P2(cos Ar,) ] 


(xsViX7y sw = (1/30)r,2r 7(3 cos? 6:;—1) 
(©: 7?) = 0 

(XV i2iZj)m = 0 
(xiy 327?) =0 


(7) 
(Xi (0 525) = O 
08? 85;—C08 85 
105 
2 cos® 6;;+3 cos 6;; 


35 





(2,72 7x :X ju =rr? 





(2°27 )w=rer# 


(2,°X;:2 579 j)m = 0 


(2;2x,2;*)» = 0. 





14 
+——15*1 'Po(cos 6;;) pr? py? P2(cos Ax,) 


+ 
+—l[r®r °P3(cos 0:5) prpyPs(cos %,) 
3R” 


+r 7 ;P1(cos 9:3) pr*p,*Ps(cos ,) }}- (8) 


The P’s are again Legendre polynomials. Here, 
obviously, the @’s with Latin indices refer to one, 
those with Greek indices to the other molecule. 
This answer, the general form of which (though 
hardly the numerical coefficients appearing in it) 


’ could have been predicted from symmetry con- 


siderations, immediately suggests the introduc- 
tion of certain invariants: 


Dipole-invariant : L=Dewriejr;P3(cos 0:5), 
ij 
Quadrupole-invariant : M = Se,r,*ej7r 7? P2(cos 644), 
ij 


Octupole-invariant: M=DLezw,*e7P3(cos 6%). 
ij 





are: 


(7) 


5 Oy) 


Pry) 


(8) 


Here, 
to one, 
lecule. 
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x in it) 
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Fic. 1. Electrical structure of the water molecule assumed 
in calculation of moments. Note: In the figure the minus 
signs should be plus, and the plus sign should be minus. 


Equation (8) then takes the form 


2 1 
(V2)y¥=—L©L@ 4+—(MML@+4L0 Ye) 
3R* R® 


14 4 j 
+o] Mom +-Love+vorey| (9) 
R*L 5 3 


The invariants may also be expressed as func- 
tions of multipole moments. On defining, as 


usual, 
dipole moment = p= >>; e:%;, 


components of quadrupole moment 


=0..=) :ex7, Qy=dDiexyi, etc., 
components of octupole moments 
=O0,.2=) :¢x3, Osey=) 1 ¢x;*;, etCc., 
they read: 
L=p*, 
M=Qee+Qhy+Q:et+(3Qry— Qs On) 
+(3Q:2—QzzOzs) + (3Qh2— Oy Oss): 
N =Oree+Oyyy+Orts +60 rey +60r22-+ 60 ryy 
+60 }y2+6022e+60y22+150 cys 
— 30 z22(Osyy + Oz22) — 30 yyy( Ozzy + Oyee) 
— 30 c22(Osse — Oyys) — 30 cacOyys 
— 30 zy c0y — 30 zyyO ze. 


(10) 


We now return to our original plan for the cal- 
culation of § exp (—V/kT)dQ. This involved 
the replacement of the correct V of Eq. (3) by 
an approximate V’ of the form 


l ‘= { u R 
where 


f=2 cos 6; cos 6.—sin 6; sin 6 cos ¢ 
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in such a way that (V?),=(V"")« as nearly as pos- 
sible. Since (f?)« = % we have 


2rL? 2Lu 
(V"? ny -{—+—“+u| 


On comparing this with (9), in which all super- 
scripts may now be omitted because we are 
dealing with two similar molecules, we find 
u=(3M/2R°) so far as the first two terms on the 
right of (9) are concerned. This choice makes the 
contribution of the 1/R'® term somewhat too 
small for the case of H:,O, but this is not a 
serious matter since that contribution turns out 
to be small anyway. Henceforth, then, we use 


L 3M 
---+ - 
R® 2R° 





V=- 


) 416.0. (11) 


In calculating the constants of Eq. (10) for the 
water molecule we have assumed the electrical 
structure proposed by Bernal and Fowler," who 
place a charge +4e at each of the H nuclei, and 
a charge —e on the bisector of the H—H dis- 
tance, 0.20A away from the 0 nucleus. The mass 
structure, needed to locate the center of mass 
about which quadrupole and octupole moments 
must be computed, is that determined by Mecke 
and Baumann.” Details are depicted in Fig. 1, 
where C represents the center of mass. They give 
the known value p=1.87X10-"* e.s.u. for the 
dipole moment. Furthermore, we find 


L=35.2K10-* es.u., M=56.6X10-* e.s.u., 
N=99.6X10-® e.s.u. 


With these values, Eq. (9) becomes 





8.27X10-72 39.9«K10-* 
(V2?) = ( + 
R* R’ 


183 X 10-™ 
+———) ergs’. (12) 
4 Rw 


At R=3A, which is near the distance of closest 
approach of two molecules, the three terms here 
retained are, respectively, 11.3, 6.1, and 3.1 


uJ. D. Bernal and R. H. Fowler, J. Chem. Phys. 1, 515 
(1933). 

12 R, Mecke and W. Baumann, Physik. Zeits. Sowjet- 
union 33, 833 (1932). K. Freudenberg and R. Mecke, Zeits. 
f. Physik 81, 465 (1933). 
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X10-*? erg? which shows that omission of the 
quadrupole moment entails a very serious error. 
On the other hand, the approximate V of Eq. 
(11) becomes 





3.521078 8.5010-* 
-( +) (0.40). (13) 


R 


It gives a (V*)y which agrees with (12) in the 
first two terms, of course, but makes the coef- 
ficient of the last term much too small (48.2 
instead of 183). This, however, will be regarded 
as acceptable in view of the inaccuracies of the 
model. 


INDUCTION EFFECT 


In describing the induction effect, the analysis 
of reference 4 will be followed. If the potential 
energy produced by one molecule at the position 
of the other is gy, the energy induced in the 
latter is —}a(V¢)*, where a is the molecular po- 
larizability. To be sure, this energy depends on 
the orientation of the inducing molecule (and to 
the extent to which a@ is a tensor upon the 
orientation of the other as well), but this 
dependence is much less pronounced than in the 
preceding effect. In particular, the energy cannot 
change sign. Instead of introducing an angular 
function in the integral ( exp (— V/kT)dQ, we 
shall therefore average over directions at once. 
This leads to the induction energy 


L 3M 
—4a((V¢)*w=—a(—+-—+---), (14) 
R® 2R8 
wherein Z and M are the former invariants. For 
the two interacting molecules, this expression 
has to be multiplied by 2. In the next section, 
reasons will be given for choosing the value of a 
to be 1.4410-** cm*. We thus obtain as the 
contribution of the induction effect to the inter- 
molecular energy in water- 


10-*° 24x10-76 
--( a + = +--+) ergs (15) 





The leading term of this expression is only about 
10 percent of the first part of (13) at R=3A, 
which indicates the relative importance of the 
two effects. 
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DISPERSION EFFECT 


The dispersion effect produces an interaction 
energy also of the form —(C/R*)—(D/R®). The 
value of C calculated in reference 4 and used in 
previous work on the second virial coefficient®® 
was obtained with the use of London’s'* formula 
in which the “resonance’’ energy was set equal 
to the ionization energy of the molecule, and a 
was taken to be 1.5X10~-*4 cm*. This procedure 
can be improved by employing optical dispersion 
formulas, two of which are available.'* 

Dorsey,!® in his valuable compendium on 
water, assumes as the best value of the polariza- 
bility a=1.59X10-*4 cm’. This is derived from 
Stranathan’s'* and other measurements of the 
dielectric constants. As might be expected in 
that type of work, the measurements are some- 
what discordant, and Stranathan admits a pos- 
sible error of 10 percent in a. 

On the other hand, Dorsey lists, without com- 
ment on the implied discrepancy, two one-term 
dispersion formulas, one due to Cuthbertson, 


2.63 X 107? 


10700 X 1077 — vp? 





n—1= 


and one proposed by Hélemann and Goldschmidt, 
3.07 X 1077 


= . (17) 
12500 X 1027 — »? 





They permit computation of a@ by extrapolation 
to v=0. Both yield a=1.44X10-** cm*. We 
believe this extrapolation to be more accurate 
than the former value and shall, therefore, use it 
in the present work. 

In a similar way, the force parameters C and 
D will be computed from (16) and (17). Equation 
(16) involves a resonance energy of 13.55 ev 
and an oscillator strength of 2.41; the respective 
quantities in (17) being 14.7 and 2.81 ev. Both 
resonance energies are smaller than the ioniza- 
tion energy (18 ev). The constant D may be 
computed from a formula previously published.” 


13 F, London, Zeits. f. physik. Chemie B11, 222 (1930). 

4 C, and M. Cuthbertson, Phil. Trans. A213, 1 (1913). 
P. Hélemann and H. Goldschmidt, Zeits. f. physik. 
Chemie B24, 199 (1934). 

146 N. E. Dorsey, Properties of Ordinary Water-Substance 
(Reinhold Publishing Corporation, New York, 1940). 

16]. D. Stranathan, Phys. Rev. 48, 538 (1935). 

17H. Margenau, J. Chem. Phys. 6, 896 (1938). 
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This procedure leads to C=33.7X10-" erg 
cm®, D=95X10-"* erg cm® on the basis of (16); 
C=36.4X10-® erg cm*, D=95X10-"* erg cm’ 
from (17). The value of D is to be regarded as 
quite uncertain. We take as the contribution of 
the dispersion effect to the intermolecular po- 
tential 





(= x10-" 95x10-76 


- - + = ) ergs. (18) 


SECOND VIRIAL COEFFICIENT 


After determination of all constituents of the 
long range van der Waals energy, attention will 
now be directed to the short range repulsive 
forces. Concerning them there is but meager 
information. Attempts to derive specific details 
from a large variety of empirical data led us to 
abandon the search for a potential the uniqueness 
of which could plausibly be accepted before 
calculating the equation of state: 


in particular the second virial coefficient B. The 
procedure has therefore been inverted: The re- 
pulsive potential is to be determined by seeking 
agreement with known values of B at different 
temperatures. 

Two facts, however, may be assumed. A survey 
of the material on the crystal structure of ice,!® 
considerations by Bernal and Fowler, the density 
of water, etc., indicate that the minimum of the 
total potential must lie somewhere between 2.7 
and 2.9A. While the concept of a minimum has 
no precise meaning for forces which depend on 
angles, this realization nevertheless proves to 
be a guide. It is possible, for example, to calculate 
the average orientation of adjacent water mole- 
cules in ice and fix the minimum with some pre- 
cision. The value of R chosen for this was about 
2.80A. Such deliberations, however, serve only 
for preliminary orientation; it develops that 
choices widely different from this do not yield 
proper values of B. 

The second fact concerns the form of the 
repulsive force law. Calculations on simpler 


1 W.H. Bragg, Proc. Phys. Soc. London 34, 98 (1922). 
B. E. Warren, J. Chem. Phys. 6, 666 (1938). 








molecules and atoms of spherical symmetry show 
that the exchange energy, in a sufficiently large 
range of distances, may be represented in the 
form A exp (—R/p), which will here be chosen. 

At first thought, then, one would try a total 
potential pieced together from (13), (15), (18) 
and the latter expression, i.e., 


V =a(R)+0(R) f(A, 62, ¢) (19) 
in which’® 


a(R) =A exp (—R/p) 








45x10-” 11910-76 
sa sen eet a 
R’ R 
3.52X10-% 8.5xX10-* 
b(R) = + 
R R5 


The parameters A and p are not quite as dis- 
posable as they might seem, even if no condition 
is imposed on the position of the minimum. The 
condition that there be a minimum at all 
drastically limits the choice. Result (19) is to 
be inserted into the well-known formula 


“ dQ ° 
B=2nN? f R%dR (1—exp (— V/kT))—. (20) 
2 


T 


The second integration yields 


dQ 
l~oe [—a(R)/AT] f exp (—bf/kT)— 


=1—exp (—a/kT)F(b/kT). 
Let b/kT =x. Then 


—xf)\dQ2 «© f (xf)*™do 
F —-, 
wnE {bias 





the last because /f*d2=0 for odd \. Now 


fra. 1 2.4+++2n = (2) . 
8x Intl 1-3---(n+1)ioX\I). 


19 Because of an oversight, a different value of the 
coefficient in the last term of a(R) has been used in the 
subsequent numerical calculations. We have, unfortu- 
nately, neglected the quadrupole moment in computing 
the induction effect, so that the coefficient in question 
was 95 instead of 119. Revision was not made because 
this term makes only a minor contribution to B and is not 
known with sufficient accuracy. 
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TABLE I. Comparison of experimental and calculated 
joe | of B in cm*/mole on the basis of the hard sphere 
model. 











T Bealc Bexp Rel. error 
700°K — 65.9 —65.4 —0.8 percent 
600°K —112.4 —99.0 — 13.5 percent 
500°K —215.2 — 166.4 — 29.3 percent 
400°K — 542 — 346.9 — 56.2 percent 








TABLE II. Comparison of experimental and values of 
B in cm*/mole calculated on the assumption of a potential 
minimum. 




















T Bealc Bexp Rel. error 
700°K — 65.9 —65.4 —0.8 percent 
600°K — 109.9 —99.0 —11 percent 
500°K — 200.2 — 166.4 — 20 percent 

so that 
bd Gon 
F(x) =>, a, (21) 
n=O (2n) ! 


While it is quite possible to integrate e~*/ without 
expansion—the result being either a series in 
Bessel functions or a definite integral over Bessel 
functions—the present form of F(x) has been 
found most convenient for problems of this type. 
We have prepared a numerical table of F(x) with 
x ranging from 0 to 10. The remaining integration, 


B=20N { Re} 1—exp (a(R) ary] -— lar 


was carried out graphically. This is much less 
tedious than expansion of the integrand and sub- 
sequent term-by-term integration. 

If A+ while p—0, V reduces to the poten- 
tial characteristic of hard spheres. This was one 
of the first models tried. It involves only one 
disposable parameter, and this is adjusted by 
making B come out correctly at one temperature. 
Measurements with which comparison is here 
made are those of Keyes, Smith, and Gerry.?® 

In the calculation with the hard sphere model, 
only the attractive terms of the potential are 
used. If the infinite barrier is placed at R= 2.87A, 
sufficient agreement results for 7=700°K. Table 
I shows what happens at lower temperatures. 

This outcome, while at first unpleasant, is par- 


20F, G. Keyes, L. B. Smith, and H. T. Gerry, Proc. 
Am. Acad. Arts Sci. 70, 319 (1936). 


ticularly interesting in view of the fact that 
Stockmayer® was able to obtain good agreement 
at all temperatures with the use of the hard 
sphere model, when he neglected the quadrupole 
moment. The discrepancies incurred can, how- 
ever, be made to tell us more about the repulsive 
interaction. They are due, it seems, to the large 
coefficient of the angular part of the potential, 
which makes a proportionately larger contribu- 
tion to B at lower temperatures. We therefore 
attempt to neutralize this effect by using finite 
values for A and p. 

A typical calculation involved the use of 
p=0.28A (a value in this neighborhood seems to 
be required in order to produce a minimum at all; 
cf. also the findings of R. Heller**) and A =3.155 
X10-* erg. These correspond to a potential 
minimum, computed with the proper value of f 
as explained, at 2.80A. The results are given in 


Table II. While the correct B is obtained at, 


700°K, the same sort of discrepancy appears at 
lower temperatures, though it has diminished 
somewhat in magnitude. Numerous trials were 
made with other sets of parameters, notably 
those corresponding to a minimum at 2.76A 
(O—O distance in ice), at 2.90A (intermolecular 
distance in liquid water), and at 3.10A. These 
yielded no better success. 

The conclusion is therefore inescapable that a 
potential of the form (19) cannot be used for the 
purpose at hand. This might indeed have been 
expected, for it implies that the function f, 
which is characteristic of polar behavior, is 
effective all the way into the origin R=0. On 
the other hand, the repulsive potential which is 
taken to be independent of angles provides no 
mechanism for counteracting this fictitious 
polarity at close distances of separation. A more 
reasonable representation of the potential ap- 
pears to be this. Assume the polar forces to be 
active for R>Rmin, the range in which the 
molecules may indeed be regarded as permanent 
multipoles. But for R< Rmin, let all interactions 
be of spherical symmetry. Though the latter 
hypothesis is but an approximation to the truth, 
it is a better approximation than the very ex- 
treme and specific one of dipolar symmetry. 

Such a representation does not yield a con- 


*tR. Heller, J. Chem. Phys. 9, 154 (1941). 
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tinuous curve for all orientations of the mole- 
cules. This analytic fault will have to be accepted 
in order to avoid the artificial difficulties which 
the construction of a continuous function would 
entail. Thus we take 


Vi =a’'(R), R< Ruin 
Ve =a(R) +b(R) f(0:0e¢)R > Ruin, 


where a’ differs from a in having a different set 
of constants A’, p’. For A and p, the preceding 
calculations strongly suggest 3.25X10-° erg, 
p=0.28A as most likely values. This leaves A’ 
and p’ to be determined. There appears to be 
some latitude in their choice, and we cannot be 
sure that the adjustment here proposed is unique. 
But the effect of reasonable alterations on B is 
not large. 

One relation between them will be obtained 
by requiring continuity between V; and V2 at 
Rmin, that is, at 2.80A for some definite tem- 
perature, say 700°K. Thus 


dQ 
8r 
=exp [—a(R)/kT ]- FLb(R)/kT | 


exp (— Vi/kT) = fe (— V2/kT) 


at T=700°K and R=2.80A. (If this fit were 
tried at R=2.90A, it would require the use of a 
negative A’. This again shows one of the self- 
stabilizing aspects of the procedure.) The second 
relation is provided by the requirement that B 
shall have its experimental value at some tem- 
perature at the middle of the range. The work 
involved here is obviously tedious, and no claim 
will be made that the present investigation is 
exhaustive; slightly different values of the con- 
stant may produce an even better fit. On taking 
A’=2.4X10-* erg, p’=0.15A—a model which 








WATER MOLECULES 315 


TABLE III. Experimental and calculated values of B in 
cm*/mole by use of adjusted constants. 











T Beale Bexp Rel. error 
700°K — 68.9 —65.4 — 5.4 percent 
600°K —101.1 —99.0 — 2.1 percent 
500°K — 164.4 — 166.4 +1.2 percent 
400°K — 329.8 — 346.9 +4.9 percent 








approaches the hard sphere case—the tempera- 
ture dependence of B (see Table III) is obtained. 
While the agreement-is not perfect, the improve- 
ment is very considerable and indicates that the 
repulsive potential adopted is not far from correct. 

The small quantum correction, computed by 
Stockmayer® for dipoles, does not improve the 
calculated values materially. An estimate for the 
present case shows that it reduces the error at 
700°K to about —5 percent while increasing that 
at 400°K to about 5.5 percent. 

Summarizing the work, we find that 


2.4X10-* exp (— R/0.15A—c(R) 
if R<2.80A 
V =< 3.25X10~* exp (—R/0.28A) 
—c(R) —b(R)(2 cos 6; cos 62 


(22) 








—sin 6, sin 6,cos ¢) if R>2.80A 
with 
45x10-" 95x10-78 
c(R) = ( + ) ergs, 
R* R’ 
3.52K10-** 8.5010-* 
b(R) -( + ) ergs 
R: R 


represents the best intermolecular potential for 
water molecules in view of available facts. 
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Measurements were made in a Townsend gap to determine the average number of photons 
liberated by an electron for each secondary electron it liberates. In He, this number was found 
to decrease rapidly with increasing ratio of field strength to gas pressure in the range from 50 
to 150 volts/cm/mm, and at pressures of the order of 1 mm. Determination of the absolute 
number of photons per electron depends on the unknown efficiency of the photoelectric cell used 
to count the photons. However, this number is estimated as of order unity. The effective radia- 
tion has an absorptign coefficient of 0.55 cm at a pressure of 1 mm of Hg. This large coefficient 
indicates that only those photons of rather high energy were counted. There is corroborative 
evidence that the brass surface of the photoelectric cell had a high work function. Thus only 
high energy photons would be expected to register. 





INTRODUCTION 


HOTO-IONIZATION is increasingly ac- 
cepted as an active mechanism responsible, 
at least in part, for the appearance of the second 
Townsend coefficient. In 1934, Cravath,! investi- 
gating the corona discharge, found that photo- 
electric effects were produced both in the gas and 
at the cathode. Dechene,? in the following year, 
showed the same phenomena to be present in co- 
rona and sparks. Kenty,’ working with the posi- 
tive column, and Greiner,‘ with Geiger counters, 
also contributed important evidence of photo- 
electric action. Photoelectrons produced in the gas 
by a spark were found by Raether® who detected 
them with a cloud chamber. The studies of Bowls® 
and Hale’ have indicated that in certain cases the 
second Townsend coefficient can be attributed 
almost entirely to photoelectric effects at the 
cathode. Recent work on corona discharge*"! 
by investigators at the University of California 
has shown the dependence of this phenomenon 
on photoelectric effects in the gas. Raether’s” 
* Now at the Applied Physics Laboratory, University of 
Washington, Seattle, Washington. 
1A. M. Cravath, Phys. Rev. 47, 254 (1935). 
2G. Dechene, J. de phys. et rad. 7, 533 (1936). 
3 Carl Kenty, Phys. Rev. 44, 891 (1933). 
4 E. Greiner, Zeits. f. Physik 81, 543 (1933). 
5H. Raether, Zeits. f. Physik 110, 611 (1938). 
*W. E. Bowls, Phys. Rev. 53, 293 (1938). 
7 D.H. Hale, Phys. Rev. 54, 241 (1938) ; ibid. 55, 815 (1939). 
8 A. F. Kip, Phys. Rev. 54, 139 (1938) ; ibid. 55, 549 (1939). 
* G. W. Trichel, Phys. Rev. 54, 1078 (1938) ; ibid. 55, 382 
(1939). 
10S, Gorril, Thesis, (University of California, 1939). 
11 Loeb, Hudson, Kip, and Bennett, Phys. Rev. 60, 714 
(1941). 
12H. Raether, Zeits. f. Physik 107, 91 (1937); ibid. 112, 
464 (1939). 
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investigations, and Loeb and Meek® in their 
book, have demonstrated that the streamer 
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Fic. 1. Arrangement of chamber. A, cathode; B, anode; 
C, “doughnut;”’ D,, De, collecting screens; E,, E2, measuring 
electrodes; F, grounded screen; G, quartz graded seal and 
window; H, Pyrex envelope. 


3... B. Loeb and J. M. Meek, The Mechanism of the 
Electric Spark (Stanford University Press, 1941). 
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“F, M. Penning, Physica 5, 286 (1938). 
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Fic. 2. Photoelectric currents at pressures of 0.70 mm, 
0.81 mm, 0.935 mm, 1.52 mm, and 1.82 mm of Hg as 
functions of X/p. 


theory of spark discharge requires the assump- 
tion of an adequate supply of photoelectrons. 
Penning" has made calculations of the efficiency 
of electrons in ionizing and exciting gas molecules, 
these efficiencies being given as functions of 
field strength divided by pressure (X/p in 
volts/cm/mm). However, he was forced to use 
data of a doubtful nature and to extrapolate the 
data to extreme values. 

In view of the incomplete nature of the above 
data, it was felt desirable that experiments be 
performed which would give an indication of the 
quantity of photons produced in a Townsend dis- 
charge. This paper sets forth the results of such 
an investigation, the gas used being He. Measure- 
ments made by Dr. L. H. Fisher in a number of 
other gases will be reported in the near future. 

Consider a function @ defined as the number of 
photons produced by an electron in traveling one 
cm in the field direction of a discharge. In analogy 
with Townsend’s a it is to be expected that 


6/p=F(X/p). (1) 
From @ and @ we may construct a function 
n(X/p) =0/a. (2) 
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Fic. 3. Photoelectric currents at pressures of 1.16 mm and 
1.34 mm of Hg as functions of X/p. 


This represents the number of photons produced 
by an electron for each secondary electron it 
produces. The function n(X /p) lends itself readily 
to the calculations to follow. 6, a, and m are 
defined for a uniform field, in which case X/p 
determines the mean energy of the electrons. 
Let zo be the initial electron current liberated 
from the cathode of a plane-parallel gap, which 
cathode is illuminated with ultraviolet light from 
a mercury arc. Then in the interval dx, at a dis- 
tance x from the cathode, ipae**dx secondary 
electrons will start. Hence, in this interval, 
nigae“*dx photons will be produced in the gas. 
The measuring apparatus will count only those 
photons which pass through a slit system and 
liberate electrons from a brass electrode. Thus 
the solid angle subtended by the slit at each point 
of the discharge must be calculated. It is assumed 
that the photons are directed isotropically. At the 
low pressures (of the order of 1 mm of Hg) 
necessary in this experiment, diffusion of the 
electrons becomes appreciable and must be in- 
cluded in the calculation. It is also necessary to 
correct for the absorption of the photons by gas 
molecules. This correction, to a sufficient accu- 
racy, is made by a factor exp [ —uoPp | where yo 
is the average absorption coefficient at 1 mm of 
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Fic. 4. Photoelectric current at 0.70 mm of Hg for X/p 
up to 280. 


Hg and p is the distance from the center of the 
discharge to the collecting electrode. 

Performance of the necessary integration leads 
to the expression 


34.0 4 
n(- -)=— ~~ oe [uopp lh(a) 
=n, exp [uopp]. (3) 


Here 7 is the photo-efficiency of the brasscounting 
surface, to be discussed later; 7 is the measured 
photoelectric current; 7,=ie” is the total elec- 
tron current to the anode; A(a) is a function 
calculated from the geometry; and m, is the 
quantity actually found by measurement, as the 
absorption coefficient is not known beforehand. 
uo can be determined from values of m, at two 
different pressures as well as from two electrodes 
at different distances from the discharge. 


DESCRIPTION OF APPARATUS AND PROCEDURE 


The chamber, originally designed by J. M. 
Meek, is shown diagrammatically in Fig. 1. A is 
the cathode, upon which ultraviolet light falls 
through the quartz window G. The anode B and 
measuring electrode system C, EZ, and D can be 
raised and lowered together by means of a screw- 
thread, ball-bearing race, and an iron armature 
turned by an external electromagnet. D is a 


screen kept at +22.5 volts to remove electrons 
liberated from the brass electrode. The “‘dough- 
nut,” containing the slit system, was kept at a 
constant potential of —67.5 volts. E, D, and C 
are insulated from each other by quartz supports. 
F is a nickel screen kept at ground potential. 77 is 
the glass envelope. The electrodes A and B are 
4.0 cm in diameter, and are parallel to better than 
0.5°. The inside diameter of the ‘“‘doughnut”’ is 
6.0 cm. Plate separations were read by means of 
a cathetometer. 

The sets of collecting electrodes E,, D,; and 
E>, Dz are at 4.5 and 6.0 cm from the center, re- 
spectively, and were used as a pair to make 
measurements of the absorption coefficient. FE, 
and Es could be connected to a Dolezalek 
electrometer through sulphur-insulated platinum 
switches. The sensitivity of the electrometer was 
1500 mm/volt, and the capacity of the system 
including E,, the cable, and the electrometer 
was 96.5 cm. 

The anode B was kept at ground potential, and 
a negative voltage applied to A. This potential 
was supplied either by batteries or by a voltage 
stabilizer. An Heraeus mercury arc was used asa 
source of ultraviolet light. The image of a slit, 
illuminated by the arc, was focused on the 
cathode. The arc was operated on storage bat- 
teries and gave a quite constant intensity. The 
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Fic. 5. np as a function of X/p for various pressures. 
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Fic. 6. The number of photons per secondary electron, 
n, as a function of X/p. Ordinates must be divided by the 

hotoelectric efficiency to obtain absolute values. Vertical 
fines indicate spread. 


visible image was rectangular, about 0.5 cm 
by 2 cm. 

The photoelectric current was measured by the 
rate of charging of the electrometer. Backgrounds 
were determined with voltages applied but ultra- 
violet light cut off. They were occasionally as low 
as 10-'* ampere, and in most cases less than 0.1 
of the photoelectric current. To keep the field 
distortion as constant and small’as possible, all 
measurements were made at the smallest plate 
distance, 1.39 cm. 

The total electron current to the anode was 
determined for each photoelectric measurement. 
This was necessary since the “‘doughnut”’ pro- 
duced considerable field distortion. The plate 
separation could not be made less than 1.39 cm 
without interfering with the beam of ultraviolet 
light. Thus, measurements on the linear portion 
of the In i,/i9 vs. plate distance curve were im- 
possible for X/p in excess of 110. Values of a 
under these conditions were determined by 
measuring 7%» at lower X /p and using this quantity 
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and the total electron current i, at the X/p 
desired. It is interesting to note that in spite of 
the large and variable field distortion, the graphs 
of Ini,/t9 vs. plate distance are unmistakably 
linear at low X/p. The current 7, was determined 
by reading the voltage drop across a 10*-ohm 
S. S. White resistor, or if too large for this, a 10’- 
ohm I. R. C. calibrated resistor. 

The chamber was baked out at 200°C and 
10-5 to 10-* mm of Hg for several hours. Meas- 
urements were made in two fillings of the cham- 
ber. Tank He was passed over Cu at 350°C and 
then through a spiral liquid air trap. The cham- 
ber, of 6 liters volume, was filled to pressures of 
1 to 2 mm of Hg in a period of approximately one- 
half hour. A McLeod gauge was used to read 
pressures. At no time was the chamber opened to 
the rest of the vacuum system without the 
presence of liquid air around the trap. 


EXPERIMENTAL RESULTS 


Measurements were made from X/p 40 to 
X/p 150. In one case the current was followed up 
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Fic. 7. The number ofphotons produced per cm per mm 
of Hg, @/p, as a function of X/p. Ordinates must be divided 
by the photoelectric’ efficiency,to obtain absolute values. 
Vertical lines indicate,spread, 
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to an X/p of 280. The observed variation of 
photoelectric current at various pressures is 
shown in Figs. 2 to 5. 

Values of , were calculated from the observed 
data and the function h(a). These values are 
given in Fig. 6 in arbitrary units. To obtain abso- 
lute quantities they must be divided by 7, the 
photoelectric efficiency of the brass counting 
surface, which is unknown. 

Independent determinations of the absorption 
coefficient, one at constant plate separation and 
the other at constant pressure, give 0.55 cm™ and 
0.42 cm—, respectively, at one mm of Hg. Correc- 
tion of , for absorption brings the curves to- 
gether, as shown in Fig. 7. The spread in values 
is indicated by the vertical lines. This spread is 
comparable with that in the measured a. 


DISCUSSION AND CONCLUSION 


Penning" has made calculations of the per- 
centages of electron energy picked up in the field 
which appear as excitation,-ionization, and in 
other forms. He has plotted these percentages as 
functions of X/p. The result of dividing the 
excitation curve by that for ionization is a varia- 
tion similar to that,in Fig. 7. His curve would 
remain at high values for low X /p. The design of 
the present chamber precludes knowledge of con- 
ditions at low a, hence, low X/p, and so this 
aspect cannot be verified. A quantitative com- 
parison with Penning’s calculations would require 
a knowledge of two additional quantities: (1) the 
photo-efficiency of the counting electrode ; (2) the 
fraction of the total excitation energy represented 
by those photons capable of releaSing electrons 
from the electrode. These quantities are not 
known for the experimental arrangement under 
consideration. 

Since the counting surface was not degassed, 
but was in equilibrium with H, at a pressure of 
one or two mm, a value of 7 approaching 0.5 is 
not unlikely. If we assume this as a rough esti- 
mate, then at an X/p of 80 and at 1 mm of Hg, 
we find @ is about one photon/cm. From ex- 
trapolated data of Townsend and Bailey,’ 
Bradbury and Nielsen,!* and the kinetic free 
"4 F, M. Penning, Physica 5, 286 (1938). 

18 Townsend and Bailey, Phil. Mag. 42, 873 (1921). 


16N. E. Bradbury and R. A. Nielsen, Phys. Rev. 49, 388 
(1936). 


path, it is found that the number of collisions 
made by an electron per cm of path in the 
direction of the field under the above conditions 
is of the order of 100. Thus the inelastic collisions 
resulting in photons of the energies here con- 
sidered represent a percent or so of the total 
number. 

Hale’s’ curves for y in He, with Pt and NaH 
cathodes, show a photoelectric peak in the 
neighborhood of X /p= 140. The absence of such 
a phenomenon in the present work may be ex- 
plained by one or both of the following: (1) Such 
a peak may not appear with a brass electrode; 
(2) the brass electrode may have such a high work 
function that the photons responsible for the 
peaks found by Hale are not sufficiently energetic 
to register. Calculations of y for the anode show 
it to be of the order of 0.001, a tenth or less of 
Hale’s values. Thus the brass surfaces present in 
the chamber are not photoelectrically sensitive. 

Raether® has reported an absorption coefficient 
in Hz of 0.9 cm~ at 760 mm of Hg. His measure- 
ments were made above 235 mm of Hg. Greiner‘ 
gives a coefficient equivalent to 1.4 cm at 760 
mm of Hg measured at pressures of the order of 
100 mm of Hg. At such high pressures any com- 
ponent of radiation with an absorption coefficient 
as large as that found in the present work would 
not be measured. It thus appears that the brass 
electrodes have surfaces of a particularly insensi- 
tive nature and are affected by only the most 
energetic photons. 

Weissler'’ has shown that in a positive point-to- 
plane corona in pure H2 no pre-onset streamers 
occur. Since the present work indicates that in 
such a discharge photons are produced in an 
appreciable quantity, it can only be concluded 
that either they are not sufficiently energetic for 
photo-ionization in the gas, or are absorbed too 
slowly to contribute photoelectrons to the forma- 
tion of a streamer. 

In conclusion, the author wishes to express 
gratitude to Professor L. B. Loeb, without whose 
guidance this undertaking would have been im- 
possible. He also thanks Dr. Leon H. Fisher for 
many helpful discussions, and Mr. E. H. Guyon, 
glassblower, for his patience and skill. 


17G. L. Weissler, Phys. Rev. 63, 96 (1943). 
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The Mesotron Momentum Spectrum at 4.35-km Altitude 


Davip B. HALL 
Ryerson Physical Laboratory, The University of Chicago, Chicago, Illinois 


(Received August 19, 1942) 


The mesotron momentum spectrum has been determined at Mount Evans, Colorado, from 
the shape of the absorption curve in lead as measured by Geiger-Miiller counter telescopes. 
Special attention was paid to eliminating the contribution of single electrons to the measured 
intensity. A large number of low energy mesotrons were found, resulting in a sharp maximum 
of the momentum spectrum lying between 1 and 2X10* ev/c. A simple power law is not 
adequate to describe the spectrum in the entire region investigated. 





INTRODUCTION the usual relation to energy of E=(p?+y?)!—y?, 
where the mass of the mesotron has been taken 
as 10® ev/c?, corresponding to 200 electron 


masses.’ 


HE mesotron spectrum at intermediate 

altitudes has usually been approximated 
in calculations by assuming a distribution in 
energy similar to that found at sea level.'~* The 
function generally used has been in the form of 
a power law with an exponent of about 3, ex- 
tending from about 10° ev to very high energies. 
In the present investigation, the distribution was 
determined by absorption measurements in lead, 
with thicknesses extending from 0.5 cm to 91.5 
cm of lead. Because the mesotrons lose energy in 
solid material by ionization processes only, and 
therefore have a definite range depending only 
on the thickness of material traversed and the 
energy of the particle, a lower energy limit can 
be assigned to those mesotrons which are able to 
penetrate any given thickness of material. The 


EXPERIMENTAL PROCEDURE 


The experiments were carried out at the 
Mount Evans Laboratory, Colorado (4.35-km 
altitude, 45.7-cm Hg pressure). Two independent 
Geiger-Miiller counter-tube telescopes shown in 
Fig. 1 were used. The individual counting rates 
were adjusted for best agreement at two points. 
In the arrangement of Fig. 1a, counter tubes 
were connected in pairs and arranged to give 
fourfold coincidence between counter pairs 1, 2, 
3, and 4, measured with a recording set having a 
resolving power of about 10~‘ second. Each set 


oS 


NOON 1 


expressions for energy loss as a function of ma- ANYON 

a YOO 2 
terial traversed, such as the Bethe-Bloch formula, : 
have been well tested and found to be reliable at . 1 O 
least in the region of 2 to 7108 ev.® The for- ° OC 
mulae given by B. Rossi and K. Greisen® were os A 
found convenient to use here. Since the absorp- 
tion curve of mesotrons represents an integral anaes 
spectrum of these particles, the usual form of MW 


differential energy spectrum can be directly 
obtained from the experimental points. For con- 








venience in comparison with cloud-chamber SSN 3 
measurements, the results to be given are NOK 4 
expressed as momenta, in units of ev/c, bearing 
1P, M. S. Blackett and R. B. Brode, Proc. Roy. Soc. “=~ a 
A154, 573 (1936). 
?P. M.S. Blackett, Proc. Roy. Soc. A159, 1 (1937). . b 


* Hayden Jones, Rev. Mod. Phys. 11, 235 (1939). 

‘ Te Hughes, Phys. Rev. 57, 592 (1940). 

y . Wilson, Proc. Roy. Soc. A172, 517 (1939). cmnemsamensmemee 

*B. Rossi and K. Greisen, Rev. Mod, Phys. 13, 240 7L. Le Prince Ringet, E. Nageotte, S. Gorodetzky, and 
(1941). R. Richard Fory, Zeits. f. Physik. 120, 588 (1943). 
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Fic. 1. Arrangement of counters and absorbers. 
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of four tubes was housed in 5 cm of lead to reduce 
the effects of soft side showers. In order to correct 
for penetrating side showers, the two boxes were 
placed at the same height so that one particle 
could not cause a coincidence. This correction, 
which was of the order of 4 percent, was sub- 
tracted from the measured intensity for each 
thickness of absorber used. No scattering cor- 
rection was made since, with the counter arrange- 
ment used, the chance for scattering into coin- 
cidence is almost the same as for scattering out 
of coincidence. Absorber thicknesses ranging from 
15 cm to 91.5 cm of lead were used in this 
arrangement. The lower limit of 15 cm of lead 
was sufficient to assure that only mesotrons were 
being detected. 

The arrangement shown in Fig. 1b was used 
with absorber thicknesses from 0 to 20 cm of lead. 
Coincidences of counters 1, 2, 3 which defined the 
solid angle subtended by the vertical telescope 
were recorded with a circuit of higher resolving 
power (r=10-*° sec.). Corrections for side 
showers were determined by measuring out of 
line coincidences with the top counter displaced. 
With less than 10 cm of lead in the telescope it 
became necessary to distinguish between meso- 
trons and electrons. This was done by placing a 
1.5-cm lead block in position A and recording 
the resultant showers of two or more particles 
emerging from the lead, coincident with a dis- 
charge of counters 1, 2,3. These fivefold dis- 
charges of the combined shower and vertical set 
could be caused by (a) knock-on electrons from 
the lead block A accompanying mesotrons, (b) 
electrons generating a cascade shower in lead A, 
(c) air showers. 

Only those events of type (b) were subtracted 
from the measured vertical intensity (coin- 
cidence 1, 2, 3). Events (c) have no relation to 
the absorber B and only provide a constant 
background counting rate for the shower set 
(4,5). These events (c) were determined by 
measuring the counting rate with the lead A 
removed. Type (a) was separated from (b) in 
order that the correction contains effects owing 
only to incident electrons. This separation was 
accomplished by assuming that 10 cm of lead in 
position B would be sufficient to stop all electrons 
which entered the apparatus as single particles. 
Those electrons of high enough energy to pene- 
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trate 11.5 cm would have high probability of 
being accompanied by additional particles when 
entering the apparatus, and hence would be 
eliminated by the displaced counter measure- 
ments. After the background was eliminated, the 


frequency of fivefold coincidences with 10 cm of ¢ 


lead in position B was taken to be owing to knock- 
on electrons accompanying mesotrons. To a 
first-order approximation at least, this type of 
event is a constant percentage of the incident 
mesotrons. The increase, then, in the rate of 
shower-producing rays with decreasing thickness 
of absorber in position B was attributed to the 
electron component of the total cosmic-ray 
intensity. Corrections were made for the change 
in efficiency of the shower tray for different 
energies of incident electrons. 


CALCULATIONS 


The efficiency of the shower-detecting tray as 
a function of the number of particles emerging 
from the lead A can be calculated from the 
geometrical arrangement of the counter tubes. 
The number of particles to be expected from 
electrons of various energies was estimated from 
Rossi and Greisen’s discussion of the cascade 
process.* In calculating: the efficiency of the 
shower tray, it was assumed that in order for 
the bottom counter 3 to be discharged, at least 
one particle must have gone through the area 
covered by the tray and discharged one of the 
counters belonging to either 4 or 5. The problem 
then becomes the determination of the compound 
probability that (a) if N particles emerge from 
the lead block A, what is the chance that they 
will emerge in such a direction as to fall within 
the area of the tray, and (b) what is the chance 
that the tray will respond to these particles since 
only half of the area is available for registering 
a coincidence. Let P, be defined as the prob- 
ability that if N particles emerge from the lead 
exactly n of them will strike the tray; let p, be 
the probability that if m particles fall in the 
sensitive area, the tray will register a coin- 
cidence. The probability of detecting an N-fold 
shower may then be written as 

N-1 
Ivy = >> Papa. 


n=l 


The value of P,, was estimated from the data on 
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cascade showers reported by Froman and 
Stearns.’ For the experimental arrangement of 
Fig. 1b, at least 80 percent of the shower par- 
ticles should emerge from the lead in such a 
direction that they will fall within the active 
area of the tray. This value was assumed to be 
independent of the total amount of absorber in 
the telescope. The error introduced by this 
assumption will tend to discriminate against the 
showers produced by high energy electrons where 
the efficiency of detection is already high and the 
frequency is expected to be low. The value of p, 
was taken as 0.5 since the six.tubes comprising 
the tray were connected so that alternate ones 
were in parallel. Any inefficiency of the counter 
tubes themselves which was of the order of 1 
percent could be neglected as a second-order 
effect, 

The initial energy of an electron required to 
penetrate a given thickness of absorber and 
emerge with one descendant was estimated from 
calculations of the cascade theory as developed 
by Arley® and also by Rossi and Greisen.!° From 
the same calculations, the number of particles to 
be expected after passing through 1.5 cm of lead 
was determined for each incident energy. The 
change in shower-producing rays observed with 
changing thickness of absorber could then be 
ascribed to certain energy ranges of incident elec- 
trons and an estimate made of the probable size 
of showers accompanying these rays. This per- 
mitted assigning an average efficiency of detec- 
tion to each range of absorbing thickness. Table 
I lists the result of these calculations; their 
application to the data is given in Table II. 


EXPERIMENTAL RESULTS 


The experimental results are listed in Table III 
and shown graphically in Fig. 2. The values of 
N(e) cannot be interpreted accurately as an 
energy spectrum of the electron component, 
since (a) air showers, and consequently high 
energy electrons, have been eliminated, and (b) 
the penetration of cosmic-ray electrons is not 
an exact measure of their energy. The value 
given for the electron intensity at 0.5 cm is an 
extrapolated figure. The mesotron intensity at 

*D. Froman and J. C. Stearns, Rev. Mod. Phys. 10, 
157 (1938). 


*N. Arley, Proc. Roy. Soc. 168, 519 (1938). 
1° B. Rossi and K. Greisen, reference 6, p. 309. 
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TABLE I. Results of calculation of efficiency of detection 
of electrons. 








No. of assoc. 





rticles of 
Total Min. energy >10 Mev 
absorber of incident below 1.5 Efficiency 
cm of Pb electron cm of*Pb assigned 
1.5 0.5 108 ev 1 
3.5 2.0 108 ev 4 ye 
5.5 5.5 X 108 ev 10 1.00 
11.5 20. K 108 ev 100 . 








TABLE II. The electron intensity corrected for efficiency 
of detection. The measured intensity was adjusted by 
subtracting the contribution owing to air showers (0.24 
c.p.m.) and knock-on electrons. This latter component 
was estimated to be one percent of the vertical penetrating 
particles from measurements with 11.5 cm of lead. 











5-fold shower Electron 

Total Vertical and vertical Electron intensity 

absorber intensity combination intensity (corrected) 
cm of Pb c.p.m, c.p.m. c.p.m. c.p.m. 
1.5 7.00 0.98 0.67 1.06 
3.5 6.10 0.73 0.43 0.46 
50 5.43 0.43 0.14 0.14 
7.5 5.02 0.32 0.03 0.03 
11.5 4.57 0.28 0.00 0.00 








TABLE III. Compilation of data. Ni is the measured 
vertical intensity together with the standard deviation of 
the measurements; N;.) is the electron correction taken 
from Table II. Ni) is the resultant vertical mesotron 
component; J(m) gives the mesotron intensity in percentages 
of the total. 








Cm Minimum 
of momentum WN (é) Std. N(m) I(n) 
Pb X10%ev/e — c.p.m. dev. N(e) corrected % 





0 0 9.50 .09 6.25 100 
0.5 0.65 8.20 .09 2.00 6.20 99.4 
1.5 0.95 7.00 .08 1.06* 5.94 95.0 
3.5 1.35 6.10 11 46 5.64 90.2 
5.5 1.70 5.43 .07 14 5.29 84.6 
7.5 2.00 5.02 .08 .03 4.99 79.9 
11.5 2.50 4.57 .08 .00 4.57 73.0 
15.2 3.00 4.33 .04 4.33 69.2 
20.3 3.70 4.15 .07 4.15 66.5 
30.5 5.00 3.49 .07 3.49 55.8 
50.8 7.80 2.83 .07 2.83 45.3 
91.5 13.75 2.08 OS 2.08 33.3 








* Note that the electron component at 1.5-cm lead corresponds to 
about 18 percent of the total mesotron intensity. Considering a mini- 
mum energy of 1 X10* ev for these electrons, this is consistent with 
results recently reported by W. E. Hazen, Phys. Rev. 65, 67 (1944). 


0 cm was taken so that the number of very slow 
mesotrons, which are heavily ionizing, was 
limited to less than one percent of the total 
intensity, in agreement with cloud-chamber 
studies at this altitude." The resulting spectrum, 
however, is not sensitive to either of these extra- 
polations. 


uC, A, Nielsen and W. M. Powell, Phys. Rev. 63, 384 
(1943). 
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Fic. 2. Absorp- 
tion of cosmic radi- 
ation at Mount 
Evans, Colorado (el. 
4.35 km). It repre- 
sents the total meas- 
ured intensity J,, 
the shower-produc- 
ing component after 
correcting for effi- 
ciency. The differ- 
ence curve J», gives 
the mesotron ab- 
sorptionadjusted to 
be 100 at 0 cm. 
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Cm Pb 


5 10 1S 20 2 30 


Figure 3 is obtained from the smooth’curve of 
Fig. 2. The units are percent of total intensity 
per hundred million electron volts. The shape of 
the entire spectrum cannot be represented‘by a 
simple mathematical function. Between 210° 
ev/c and 10° ev/c, however, a rough fit can be 
obtained for the momentum spectrum of meso- 


o &¢& ©» 8&8 8© 6&6 H 


trons by dN/dP =Ap~-*, where A represents an 
arbitrary constant. In view of the very different 
spectrum observed at sea level, there is no jus- 
tification in extending this analytical expression 
to predict the shape of the curve beyond the 
region covered by the experimental data. 

The presence of the prominent maximum in 


Fic. 3. Differential momentum 
spectrum at Mount Evans (4.35 
km). Abscissa is momenta values 
obtained by range-momenta rela- 
tionships. Ordinate is the percent 
of total intensity per interval of 
108 ev/c. 
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the region between 1 and 2108 ev/c cannot 
easily be accounted for on the usual basis of a 
mesotron production layer at 16 km and sub- 
sequent spontaneous decay in the atmosphere. 
The shift of the maximum to lower energies and 
the large number of mesotrons between 10* and 
10° ev/c would require special assumptions in 
regard to the distribution with which they are 
created. It is highly probable, therefore, that 
many of these mesotrons are created relatively 
close to the point of observation (Mount Evans). 
Since the number of protons observed near this 
altitude" does not seem large enough to be the 
principal agent of production, neutral radiation, 
such as photons or other neutral particles, may 
be responsible for most of these low energy 
mesotrons. Although these production processes 
at low altitudes have been discussed before,!? !* 
the importance of their contribution in the de- 
scription of the spectrum has not been em- 
phasized. 

If the observed momentum distribution of 
Fig. 3 is used to calculate a mesotron spectrum 
for the elevation of Echo Lake, Colorado (3.24 
km), using the usual formulae for decay and 
ionization loss, one obtains a curve shown in 
Fig. 4 that is in reasonably good agreement with 
the data obtained by anticoincidence measure- 
ments of Rossi, Greisen, Stearns, Froman, and 
Koontz'* taken at this station. This would 
indicate that the mesotron production processes 
do not contribute an appreciable quantity of 
mesotrons in the energy range between 3108 
to 12 10* ev for altitudes intermediate between 
that of Mount Evans and Echo Lake. 

Since only mesotrons of momenta greater than 


2 M. Schein and V. C. Wilson, Rev. Mod. Phys. 11, 292 
(1939); M. Schein, W. P. Jesse, and E. O. Wollan, Phys. 
Rev. 57, 874 (1940); M. Schein, E. O. Wollan, and G. 
Groetzinger, Phys. Rev. 58, 1027 (1940), 

’V.H. Regener, Phys. Rev. 61, 105 (1942). 

“B. Rossi, K. Greisen, J. C. Stearns, D. K. Froman, 
and P. G. Koontz, Phys. Rev, 61, 678 (1942). 
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Fic. 4. Mesotron momentum spectra at several alti- 
tudes. Curve (a) is the measured curve for Mount Evans 
(4.35 km), curves (b), (c), and (d) are calculated from 
curve (a) for altitudes of 3.24, 1.6, and 0 km considering 
ionization losses in air and mesotron decay with ¢9 = 3X 10-* 
sec. The rectangles shown in broken lines are the data of 
Rossi et al. taken at Echo Lake (3.24 km). 


1X 10* ev/c at Mount Evans can reach sea level, 
the spectrum as measured cannot be satis- 
factorily used as a basis for calculating a sea- 
level distribution of mesotrons. The position of 
the maximum, however, would be expected to 
occur at about 510* ev/c if mesotron produc- 
tion processes below 4-km altitude can be dis- 
regarded. The existing cloud-chamber experi- 
ments at sea level have tended to discriminate 
against such low momentum particles, producing 
a false maximum in the region of 10° ev/c. 
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The existence of an empty (conduction) band 1.5 or 2.0 
volts above the highest filled band in metallic-looking 
semi-conductors such as silicon and certain sulfides is 
inferred from the optical properties of these substances. 
Reflectivity data of silicon are examined with the aid of 
the expression for the complex index of refraction as given 
by classical electromagnetic theory. Through use of this 
expression, which is checked for internal consistency, an os- 
cillator strength of 1.6electrons is calculated for the absorp- 
tion, and a dielectric constant at low frequencies of 12.5 
is computed. The Wigner-Seitz-Slater method of computing 
electronic energy bands in crystals is used to determine the 
band structure of silicon. A combination of one s, three , 
and three d functions for each of the two atoms in the unit 
cell is made, and through use of boundary conditions of 
continuity of value, of normal, and of tangential deriva- 


tives, a solution is obtained for the plane x= y. The result- 
ing band structure closely resembles that obtained by 
Kimball for the diamond, except that silicon is more nearly 
metallic than is the diamond. The 3s level of the silicon 
atom splits into two bands, the 3p level into six bands, and 
the 3d level into six bands; the two 3s bands and the 
lowest two 3p bands are completely filled by electrons. At 
the observed half-distance between nearest neighbors, the 
gap between the uppermost filled band and the lowest 
empty band is much greater than that expected from con- 
sideration of the optical data. The form of the optical ab- 
sorption as expected from band structure considerations is 
proposed. The width of the filled bands as observed in soft 
x-ray emission spectra is about equal to the width of the 
computed bands. 





I 


HERE is in addition to the metals another 

group of substances which possesses the 
so-called characteristic metallic luster. In this 
group are the elements silicon and germanium, 
and many minerals, mostly sulfides, such as 
pyrite (fool’s gold, FeS,), molybdenite (MoS,), 
galena (PbS), and stibnite (Sb2S;). All of these 
substances conduct electricity feebly at room 
temperature and may be classed electrically as 
semi-conductors. To our knowledge there are no 
substances which look metallic and which are not 
metals or semi-conductors. 

The question at once arises what connection 
exists between metallic appearance and con- 
ductivity. It is first necessary to determine what 
is meant by metallic appearance. Inspection of 
reflectivity data! reveals that only metallic- 
looking substances such as the previously-men- 
tioned semi-conductors have reflection coef- 
ficients above 35 percent throughout most of the 
visible spectrum. There are some substances such 
as the aniline dyes and KMnQ, which reflect 
metallically in narrow wave-length regions; in 
these, the strong reflection is brought about by 

* A dissertation submitted to the Faculty of the Grad- 
uate School of Arts and Sciences of the Catholic University 
of America in partial fulfillment of the requirements for 
the degree of Doctor of Philosophy. 


1 International Critical Tables (McGraw-Hill Book Com- 
pany, Inc., New York, 1929), Vol. V, p. 254. 


a definite, allowed electronic transition within 
the molecule or ion. Our interest, however, is in 
metallic reflection extending over a wide range 
of frequencies, over most of the visible and per- 
haps into the ultraviolet. 

Careful examination of the equations relating 
absorption and reflection shows that one cannot 
get a reflection coefficient of 35 percent or more 
over a wide range except by transitions corre- 
sponding to a total electron strength of the order 
of unity. In metals there are enough electrons in 
the conduction bands and these overlapping 
bands are sufficiently broad to permit within the 
bands transitions which give sufficient reflection. 
It is obvious, however, that in semi-conductors 
the conduction electrons, which are few in 
number and cannot contribute appreciably to 
the absorption, do not play an important part in 
the high reflectivity of many members of this 
class of substances. 

In semi-conductors, the metallic appearance 
must be due in general to the existence of a strong 
(i.e., f~1), wide absorption band covering the 
visible. This can, in the simplest manner, be 
ascribed to the existence of a broad, empty (con- 
duction) band with a lower limit 12,000—15,000 
wave numbers (1.5-2.0 volts) above the upper- 
most filled band.? 


2Such a view has been proposed by F. Méglich and 
R. Rompe, Zeits. f. Physik 119, 472 (1942). 
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From this brief consideration of the conduc- 
tivity and reflectivity of metallic-looking semi- 
conductors, certain conclusions seem to be jus- 
tified. There exists in these semi-conductors an 
empty (conduction) band, with a lower limit 
1.5-2.0 volts above the highest filled band, and 
a width of several volts. Transitions into this 
band from the uppermost filled band are allowed, 
and give correspondingly strong absorption and 
reflection in the visible. At high temperatures, 
some electrons possess enough energy to go 
across the 1.5—2.0 volt forbidden region between 
the highest filled band and the empty band; 
these high energy electrons in the conduction 
band and the holes produced in the filled band 
by their leaving are responsible for the intrinsic 
conductivity of these substances. Thus the 
nearness of this empty band to the filled band 
both makes these solids intrinsic semi-conductors*® 
at high temperatures, and gives them metallic 
appearance. 


II 


It is the purpose of this section to examine the 
refractive indices , and the coefficients of ab- 
sorption mx, of silicon in the visible region as 
determined by experiment to provide quanti- 
tative evidence of the existence of a band 
structure as described in the first section. The 
experimental values of m and « are to be used to 
find the oscillator strength of the band. Then, as 
a check, theoretical values of n?—n*x? are to be 
calculated by use of an expression developed for 
band absorption through the classical theory of 


dispersion. The optical properties of silicon are. 


also to be used to calculate the dielectric con- 
stant of silicon at low frequencies. 

For a single absorption line, it can be shown‘ 
from the classical theory of dispersion that at 
the frequency v 

Ne? f 
n?(1 —1x)?=1+— ; (1) 


° ; 
rm ve—v?+ivy 





where N is the number of atoms per cubic cen- 
timeter, vo the frequency of the absorption line, 
vy’ a quantity measuring its width, and f the 


An intrinsic semi-conductor is a substance which in the 
pure state has a small electronic conductivity that increases 
greatly with rising temperature. 

‘See, for example, Slater and Frank, Introduction to 
Theoretical Physics (McGraw-Hill Book Company, Inc., 
New York, 1933), p. 280. 
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“electron number” or oscillator strength to which 
the strength of the line is proportional. Classi- 
cally, f should be 1. 

The absorption coefficient mx is defined so that 
the amplitude of light of frequency », after 
penetrating a thickness x of the material, is 
diminished by the factor 


2rvnxx 
gies |- | 
Cc 





If, as in silicon, there is a broad absorption band 
instead of a single line, the relationship (1) can 
be used if the oscillator strength is distributed 
over the band; f is replaced by (df/dvo)dvo, the 
oscillator strength in the interval between vp» and 
vot+dyo, and the second term of the right side 
of (1) is integrated over the entire band. Ex- 
pression (1) then becomes 





Ne 1 df 
n?(1—ix)? =1+— f —dv. (2) 
Vv 


™m o—vt+ivy’ dvo 


If the real and imaginary parts of (2) are 
separated, the following equation is obtained for 
the imaginary part: 

Ne vy df 
12n?x =i—— —dy. (3) 


™m (ve? — v?)?+ vv”? dyo 





The quantity »’, which measures the width of 
the absorption due to an oscillator element of 
strength df at a given yp and arises from radiation 
damping, can be assumed to be vanishingly 
small, since the radiation damping can be shown 
to be negligible. The limit, then, of (3) can be 
taken as v’ approaches zero, and the result can 
be solved for df/dvo. In this way it is found that 
at the frequency v» the value of df/dyvo is given by 

df 4m 

— =—-n' xv. (4) 

dy, Ne 


To evaluate df/dvo, it is necessary to have 
measured values of m and mx throughout the 
band. Using a method developed by Drude,‘ 
Pfestorf? has determined values of m and nx at 
intervals throughout the band. From Pfestorf’s 
data, which are presented in Fig. 1, the values 

5H. M. O'Bryan, J. Opt. Soc. Am. 26, 125 (1936). 


6 P, Drude, Wied. Ann. 64, 162 (1898). 
7G. Pfestorf, Ann. d. Physik 81, 925 (1926). 
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of nx listed in Table I and plotted in Fig. 2 were 
determined. From the plotted values of m?x in 
(4), df/dvo was calculated at various frequencies 
as presented in Table II and the second curve of 
Fig. 2. As shown in Fig. 2, the broad absorption 
band lies well within the given frequency limits. 
In the following, the letter * indicates wave 
numbers, and df/di=(1/c)(df/dyv). 

By integrating (df/dv)d> over the entire band 
using numerical methods, one obtains® an oscil- 
lator strength of 1.6. This absorption is much 
too strong to come from transitions to or from 
an impurity band. The simplest assumption is 
that the absorption is brought about by transi- 
tions from a full band to an empty (conduction) 
band of the lattice. 

Attention can now be turned to the real part 


of (2), 





—d yo. (5) 


(vo? — v?)?+ vv’? do 


Ne? f (v2?—»*) df 


mm 











l l 
° sats 2 4,000 A SPOOA 








WAVELENGTH 


Fic. 1. Refractive index and coefficient of absorption 
of silicon; data of G. Pfestorf, Ann. d. Physik 81, 925 
(1926). N in the figure should have been n. 


If df/dvo were known for all frequencies, the 
integral in (5) could be evaluated with the 
limits zero and infinity. As it is, the values of 
df/dvo have been calculated through use of 
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optical data, so that integration can be done only 
over the frequencies of the visible range. How- 
ever, the assumption can be made that all other 
absorption bands lie so far in the ultraviolet that 
they do not contribute to the absorption in the 
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Fic. 2. Oscillator strength df/d> and nx of silicon. N in 
the figure should have been n. 


visible range, and add only a frequency inde- 
pendent quantity m?—1 to the expression for 
n?—n’x* for wave-lengths longer than 2500A. 
As the limit v’ approaches zero, the expression 
for the refractive index m and the coefficient of 
absorption nx at the wave number 7, is given by 


Ne 1 af 
nt—ntt=ni+— [ —_—as, 6 


rmc* 





in which 7 is the wave number at which the oscil- 
lator strength (df/dv)d7 is located, and 


Ne? 1 








d 
Ns (7) 


amc ¥*? — 7," dv 


is the contribution this oscillator strength makes 
to n?—n*x* at %,; the integration is performed 
over the visible region using graphical methods. 
The values of this integral for various wave 
numbers %, are listed in the fourth line of Table 
III; in the third line are given the values of 
n* —n*x* determined by experiment. The last line 
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is the difference between the third and fourth 
lines, and should be equal to m,?. The data are 
plotted in Fig. 3. 

The values of the differences are all positive, 
as they must be if they are to represent m,*; they 
show no reasonable trend. It can be assumed 
that the variations are due to experimental 
uncertainties, and as the value of m,? the arith- 
metical average may be taken: mo? =2.6. 

From these data it is possible to calculate the 
dielectric constant. The electrical properties of 
a material like silicon depend on the bound 
electrons and on the free or conduction electrons. 
For the dielectric constant as used in space 
charge formulae, the contribution of the con- 
duction electrons is to be omitted, as the con- 
tribution of these electrons to the electric field 
has already been taken into account in the space 
charge, and would be counted twice if it were also 
included in the dielectric constant. 
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Fic. 3. Observed n?—n*x* and computed value of integral 
(7). NV in the figure should have been n. 


If the refractive index m is known at zero fre- 
y 


quency, the dielectric constant e{can be found 
from 


lim n? =e. 
v0 


Placing » in (6) equal to zero, one finds that the 
expression for the dielectric constant at low fre- 


quency becomes 
id Y.. 
—foa (8) 
med 7 rr 


Graphical integration of the second term leads 
to a value of 9.9, and therefore to 


e=12.5. 


Apart from the uncertainty of the optical data, 
which is difficult to judge, the accuracy of e¢ is 
estimated to be +0.5. 

In the preceding, any influence of the con- 
duction electrons on the optical properties has 
been completely omitted. If these are completely 
free, their contribution to m*? would be — (Ne?) / 
(rmc*)1/7?. With Ni=6X10", *=20,000, the 
value of the above is —0.024, which is below 
experimental accuracy. 


Ill 


prem... 





Quantitative investigation of the electronic 
structure of solids has been confined largely to 
studies of the structure of metallic and ionic 
crystals. The atoms which form valence crystals 
have so many outer electrons which are affected 
appreciably by entrance into the crystalline 
state that computation involved in determining 
approximately the electronic structure becomes 
difficult. As a consequence of this difficulty, the 
work on valence crystals has not been extensive, 
and no band structure studies of the silicon 
crystal have been made. In this section, a method 
of quantum mechanical approximation is to be 
used to see if the origin of the observed absorp- 
tion of silicon can be found in the band structure 
as given by this approximation. 

A method of obtaining approximate solutions 
of the Schrédinger equation which has been used 
successfully for valence crystals is the well-known 
method developed by Wigner and Seitz* and 
extended by Slater.* Wigner and Seitz observed 
that in a crystal having a high degree of rota- 
tional symmetry relative to the nucleus, the 
potential field acting on an electron in the 
neighborhood of the nucleus is spherically sym- 

SE. Wigner and F. Seitz, Phys. Rev. 43, 804 (1933); 
46, 509 (1934). 


*J. C. Slater, Phys. Rev. 45, 794 (1934); Rev. Mod. 
Phys. 6, 209 (1934). 
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metrical to a good approximation. As in the case 
of isolated atoms, the Schrédinger equation can 
be solved by separation of variables and by 
numerical integration of the radial equation for 
the radial function. Each nucleus in the lattice 
can be considered to be the center of a roughly 
spherical cell which is bounded by planes so that 
together with the cells about other nuclei it will 
fill all space. The wave function for the entire 
lattice can be obtained by solving the Schrédinger 
equation in one of the cells, subject to certain 
conditions at the cell boundary. This method 
was used in the determination of the electronic 
energy bands of the diamond by Kimball!® and 
by Hufd" and Mrowka.” 

The structure of the silicon crystal like that 
of the diamond is face-centered cubic with two 
atoms in the unit cell; the edge of the unit cell 
is 5.42A in length,!* and the distance between 
nearest neighbors is 2.35A or 4.44 Bohr radii. As 
in the work of Kimball,!° one atom of the unit 
cell is considered to be at the origin of a Cartesian 
coordinate system and the other at the point 
(a, a, a); the entire lattice may be formed by 
translating these atoms of the unit cell by vectors 
of the form ¥,;-+-2V2+3;v; where the com- 
ponents of v; are 0, 2a, 2a; those of v2 are 2a, 
0, 2a; and those of v; are 2a, 2a, 0; and where 
m1, M2, and m3 are arbitrary integers. 

The Wigner-Seitz polyhedron for silicon is 
made up of a tetrahedron formed by the planes 
bisecting perpendicularly the lines joining the 
atom to its four nearest neighbors. The corners 
of this tetrahedron are cut off by the planes 
bisecting perpendicularly the twelve lines to the 
next nearest neighbors. The result is a sixteen- 
sided polyhedron having four large hexagonal 
faces.14 

As there are two atoms in the unit cell of 
silicon, it is necessary to consider the be- 
havior of wave functions in two Wigner-Seitz 
polyhedra simultaneously. The polyhedra sur- 
rounding the two atoms of the unit cell are 
identical, but are oppositely oriented. The wave 


10 G. E. Kimball, J. Chem. Phys. 3, 560 (1935). We often 
use Kimball’s notation. 

1 F, Hund, Physik. Zeits. 36, 888 (1935). 

12 F, Hund and B. Mrowka, Ber. d. Sachs. Akad. d. Wiss. 
math-phys. KI. 87, 185 (1935). 

18H. Kiistner and H. Remy, Physik. Zeits. 24, 25 (1923). 
a A drawing of this polyhedron is presented in reference 
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TABLE I. Experimental data. 




















A in A 5780 5460 4360 4060 3660 3350 2810 2540 

v 17300 18300 22930 24630 27320 29850 35600 39370 

nx 3.5 3.9 94 11.5 12.5 11.0 7.9 3.7 
TABLE II. Oscillator strength. 

d 16000 18000 20000 22000 24000 26000 28000 30000 32000 

55 X108 is 20 30 SO 75 0 9% £93 ~~ 89 

5 36000 38000 40000 4”700 


7p X10 78 59 37 15 








TABLE III. Comparison of experimental and theoretical 
data for n?—n?x?, 








r 5451 4286 3623 3331 2998 
"1 18300 23300 27600 30000 33300 
n?— nx? 17.5 17.6 10.0 —1.5 — 6.8 
Integral (7) 16.8 13.3 4.7 —3.1 —8.1 
Difference 0.7 4.3 5.3 1.6 iz 








function can be considered to be made up of an 
even part U* and an odd part U*. The even 
part U* remains unchanged in sign in going from 
the midpoint of a hexagonal face of one poly- 
hedron to the midpoint of the diametrically 
oppositely oriented face of an adjacent poly- 
hedron; the sign of the function and the orien- 
tation of the face in each case are determined 
using a coordinate system having its origin at 
the center of the polyhedron being used. The 
odd part U“ changes sign in such translations. 
For convenience, the wave function U is written 
in the form 
U=U9+1U". 

In the study of the electronic structure of the 
diamond, it was necessary to consider only the 
L-shell electrons in which only s and ? states 
exist. A combination of eight wave functions and 
the eight simple boundary conditions of con- 
tinuity of the wave function and its normal 
derivative at the four pairs of midpoints suf- 
ficed. In silicon, the electrons of interest are in 
the M shell, and it can be expected that d states 
will deserve consideration. 

The geometry of the Wigner-Seitz polyhedron 
for the diamond lattice is such that only the 
centers of the four hexagonal faces can be used 
for application of the boundary conditions. Any 
combination of wave functions which is to take 
into consideration the s, p, and d levels must 
contain more than eight functions, the maximum 
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that can be determined by use of the usual 
boundary conditions. It is necessary, therefore, 
to introduce additional boundary conditions. 

It was first attempted to include all five 3d 
levels in the combination of wave function. Ac- 
cordingly, the odd and the even part of the total 
wave function were made up of the 3s function, 
the three 3p functions, the five 3d functions, and 
one 4f function. As there were four points in each 
polyhedron at which the boundary conditions 
could be applied, it was thought that the number 
of terms in the total wave function must be a 
multiple of four; it was for this reason that an 
f function was added to the odd and even parts 
of the total wave function to give twenty coef- 
ficients. Five boundary conditions at each of the 
four joining points were sought. The boundary 
conditions selected were the usual conditions of 
continuity of the value and of the first normal 
derivative of the functions; the new conditions 
were the continuity of the second normal de- 
rivative, and the continuity of the two com- 
ponents of the tangential derivative.’® Solutions 
of the twenty simultaneous equations were 
sought for the condition’”® M=L=0. It was found 
that the equations were not all independent and 
would not give twenty solutions. No way of avoid- 
ing the difficulty was found, and it was decided to 
abandon the attempt to use twenty functions. 

A revision of the boundary conditions was then 
made in an attempt to get as many independent 
equations as there were functions. The revised 
boundary conditions were the continuity of the 
function and of first, second, and third normal 
derivatives. The total wave function was as 
before, except that only the three 3d functions 
having the surface harmonics xy/r’, yz/r?, xz/r’, 
were used. Independent solutions were found for 
the condition M=L=0, but when numerical 
values were calculated, the solutions appeared 
completely unreasonable. It appears that the con- 
tinuity of normal derivatives of order higher than 
the first are not satisfactory boundary conditions. 

The conditions of continuity of the second and 
third normal derivatives were then replaced by 
the continuity of the two components of the 
tangential derivative. The same combination of 


18 Continuity of the tangential derivatives was suggested 
as a boundary condition by W. Shockley, Phys. Rev. 52, 
866 (1937). 
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wave functions was used as in the preceding case. 
The first attempt at solution was made for the 
condition L=M=0, but it again appeared that 
there were too few independent equations to 
determine the s, ~, and d function coefficients. 
However, when the equations were examined for 
the condition which Kimball*® uses, K = L, it was 
found that fourteen unknowns could be deter- 
mined from the sixteen equations. Accordingly, 
the 4f functions in the odd and in the even part 
of the total function were dropped, and the 
values of the coefficients of the remaining four- 
teen s, p, and d functions were determined from 
the sixteen equations. As this procedure was used 
as the basis of the numerical work which follows, 
the solution is given in detail. 

The boundary conditions of continuity of the 
value of the wave function, of the normal de- 
rivative, and of the two components of the 
tangential derivative may be stated more ex- 
plicitly as :1* 

16 The form of the differential operator in the expressions 
giving the components of the tangential derivative may 
be developed as follows: at the midpoint x;, y:, 2: of a 
hexagonal face of the polyhedron let there be a normal 
vector with direction cosines x;/f1, yi/r1, and 2:/r1, where 
r, is the distance to the origin, which is at the center of 
the polyhedron. The normal vector to the midpoint of the 
face will pass through the origin. In the plane of the face 
let there be a second vector through the midpoint with 
components dx, dy, and dz and direction cosines dx/ds, 


dy/ds, and dz/ds, where ds = (dx*+dy*+dz*). Then from 
the familiar relationship 


xi/r, dx/ds+yi/r1 dy/ds+2:/r; dz/ds=0 (A) 


it follows that 
x,dx+y,dy 
21 7 


The total differential of F(x, y, z) at x1, 91, 21 is 


dz=— (B) 


OF oF OF 
CPo et eta (C) 


Using (B), we may then express the tangential differ- 
ential as 


OF _ #1 9F oF *\d 
6F= G Z oa + dy 2,02)" (D) 
The wave functions listed in Table IV involve x, y, z, 


and r 
F(x, y, 2) =(x, ¥, 2, 7) 


where r = (x*+-y*+-2*)!. The partial derivatives of the wave 
function then are: 


()- [=] + =| = (E) 
(=) -[3)... +2171... ©) 
a Es ee Eo (G) 


The differential quotients ‘enclosed in square brackets 
[ ] refer to differentiation of the variable as it appears 
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Ue+iUe =e'®-)(Ue +i), (9a) 
U'e+iU =e V(U'2 i), (9b) 


({;-|- “5 [eto 
--eon((S}z)pervon. 
(Fie 
--een([=]- *[--]))eti09. a 


In the boundary conditions (9c) and (9d), the 
differentiation set in square brackets [ ] involves 
only the variable indicated as it appears ex- 
plicitly in the function operated upon; r is con- 
sidered as though it were independent of x, y, 
and z. 

In the above, k is the wave number vector, and 
V is the translational vector from the midpoint 
of one hexagonal face to the midpoint of the 
diametrically oppositely oriented hexagonal face 
of an adjacent polyhedron. By combining terms, 
Eqs. (9) reduce to: 


explicitly in the function operated upon, with no regard 


to the dependence of r on x, y, z. 
On substitution of (E), (F), and (G) in (D), the tan- 
gential differential at the point x1, y:, 21, takes the form 


= ( [| 4 22/94] 21 84] _ 21 xf ov 
ar=([%]+ =(*]- =| Ti = | ae 
ay] , fay] _ nm 21 vil dy 
+( a+ {FI- te |-# ak *|)2 ” 


When the expressions by which dx and dy are multiplied 
fulfill the condition of continuity at a face of the poly- 
hedron, the tangential derivative is continuous in any 
direction and the function is continuous not only at the 
midpoint but in its neighborhood. For want of a better 
name, the coefficients of dx and dy are referred to in the 
text as the components of the tangential derivative. 

As can be seen from (H), the operators involved in 
these tangential boundary conditions are 


([z-]-SLe) @ 
([3]-Zle)) @ 


_ Thus the so-called components of the tangential deriva- 
tive do not involve explicitly differentiation of ¥(x, y, z, r) 
with respect to r. 


and 
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07] x7 a7 07 07 
LOxt 2,Ldz ay] z,Laz! 
U's U* 
=— =—=tan—. (10) 
U' Us 2 


The joining points of the polyhedron about the 
atom at 0, 0, 0 are: Pi = (a/2, a/2,a/2), P2=(a/2, 
—a/2, —a/2), P3=(—a/2, a/2, —a/2), and 
P,=(-—a/2, —a/2, a/2). For the polyhedron 
about the atom at (a, a, a) the joining points are 
P,=(-—a’/2, —a’/2, —a’/2) or (a/2,a/2,a/2), 
P,=(—a’/2,a’/2,a'/2) or (a/2, 3a/2, 3a/2), 

=(a’/2, —a’/2,a’/2) or (3a/2,a/2, 3a/2), 


and 
Ps=(a'/2,a'/2, —a’/2) or (3a/2, 3a/2, a/2). 


TABLE IV. Values of independent wave functions. 











Wave Value in polyhedron Value in polyhedron 
function about (0, 0,0) about (a, a, a) 

v1 S(r) — S(r’) 

v2 V3 (x/r) P(r) V3 (x’/r’) P(r’) 

V3 v3 (y/r) P(r) V3 (y'/r’) P(r’) 

v4 v3(2/r) P(r) v3 (2'/r’) P(r’) 

vs 3(xy/r?) D(r) — 3(x’y’/r) D(r’) 
vs 3(yz/r*)D(r) — 3(y's’/r") D(r’) 
v7 3(xz/r?)D(r) — 3(x’'s'/r'*)D(r’) 
Vs S(r) S(r’) 

ve V3 (x/r) P(r) — v3 (x'/r') P(r’) 
v0 V3(y/r) P(r) — v3(y'/r’) P(r’) 
vir V3 (2/r) P(r) — V3(2'/r’) P(r’) 
Vie 3(xy/r?) D(r) 3(x’y'/r)D(r’) 
V3 3(yz/r*) D(r) 3(y’2'/r'?) D(r’) 
vis 3(xz/r?)D(r) 3(x’s'/r’2) D(r’) 








In the above, the primed quantities refer to a 
coordinate system with origin at a, a, a and the 
unprimed quantities to a system with origin at 
0, 0,0. If the boundary conditions are applied 
to these two adjacent cells, the midpoints of the 
oppositely oriented faces and the translational 
vectors V involved are: P; and P; with V=0 (the 
faces coincide); P, and Ps with V=v,; (com- 
ponents 0, 2a, 2a); Ps and P; with V=v, 
(2a, 0, 2a); and P, and Ps with V =v; (2a, 2a, 0). 

To satisfy the boundary conditions at the mid- 
points, the wave equation 


U = 1a, +taw. + 1a3/3 + tas + lass 
H+1aetiaw:t+bist bot bio 
tbautbspiet bois tou 
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was built up of the fourteen independent func- 
tions whose values in the two polyhedra are 
listed in Table IV. 

In Table IV, S(r), P(r), and D(r) are the 
radial parts of the wave functions; x, y, and z are 
Cartesian coordinates with origin at (0, 0, 0) and 
r is the distance from (0, 0, 0); x’, y’, 2’, and 7’ 
are similarly related to the point (a, a, a). 

The equations obtained from the boundary 
conditions are: 


$(a1) + p(a2+a3+as)+d(as+ag+a7)=0, (11a) 
$(a1) + p(a2—a3 — a4) +d( —a5+a¢4—a7) 
+K[s(bi) + )(b2— bs — by) 

+d(—bs+b.—b7)]=0, (11b) 
$(a1) + p( -a2+a3— a4) +d(—as—ag+a7) 
+L[s(bi1) + p(—b2+b3— da) 

+d(—bs;—be+b7)]=0, (11c) 
$(a;) + p( —a2—a3+ 4,4) +d(as5—a5—a7) 

+ M[s(bi)+p( —b2—b3+0,) 

+d(bs—be—bz) ]=0, (11d) 
s'(b;) +p’ (be +b3+b,4) +d’ (bs +b6+b7) =0, (11e) 
—K[s’(a1)+p'(d2—a3—a4) . 

+d’ (—as+d5-—-a7) ]+s' (1) 

+ p’ (be—b3— by) +d’(—bs+b.—b;)=0, (11f) 
—L[s'(a:) +p’ (—a2+a3—a4) | 
+d’ (—as—adg+az) ]+s’(d;) 
+p’ (—be+bs—bs) +d’ (—bs—be +7) =0, (11g) 
— M[s’ (a1) +p’ (—a2—a3+<4) 
+d’ (as—as—az) ]+5'(b1) 

+p’ (—b2—bs+bs) +d’ (bs—be—b7)=0, (11h) 

p(b2—bs) ++d(bs—be)=0, (111) 

— K[p(a2+as) +d(—as—ag) ] 
+ p(be+b,)+d(—bs—be) =0, (113) 

—L[{p(a2—a4) +d (a5 —a¢) ]+p(b2—d,) 

+d(bs—be)=0, (11k) 

— M[p(a2+a,) +d(—as5—a¢) ] 
+ p(be+b,) +d(—bs—bs)=0, (111) 
b(bs—b,) +d(b;—b;) =0, (11m) 
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—K[p(a3—as) +d(as—az) J 
+ p(bs—bs)+d(bs—b;)=0, (11n) 
—L{p(as+a,)+d(—as—a) ] . 
+p(bs+hs)+d(—bs;—b;)=0, (110) 
— M[_p(as+a,) +d(—as—a7) ] 
+p(bst+h.)+d(—bs—b;)=0. (11p) 


In Eqs. (11) and (18), the @a’s and b’s are fac- 
tors, not arguments. 
In the above equations, 


s=S(R), s’=(dS/dr),-r, 
p=P(R), pb’ = (dP /dr),.r, 
d=D(R), d’=(dD/dr),.r, 


where R is the distance from the center of the 
polyhedron to the midpoint of the hexagonal 
face. The definitions of K, L, and M are 


K=tan (k-v,/2), L=tan (k-v,/2), 


and M=tan (k-v;/2). 

Equations (11a) to (11d) express the condition 
of continuity of the wave function itself at the 
four selected pairs of points, (1le) to (11h) that 
of the continuity of the normal derivative, (11i) 
to (111) that of the “x component” of the 
tangential derivative, and (11m) to (11p) that 
of the ““y component.” 

To make solution of these equations feasible, 
propagation of electron waves is restricted to the 
plane x=y, i.e., k,=0. Then, as in Kimball's 
solution,'“° K=L. The sixteen simultaneous 
equations can be broken into two sets of eight 
equations. The first set, obtained from the 
boundary conditions of continuity of the value 
and of the normal derivative, contains s, p, d, 
s’, p’, and d’, and the fourteen constants. The 
second set of eight equations, obtained from the 
boundary conditions of continuity of the com- 
ponents of the tangential derivative, contains 
only p and d and the twelve constants associated 
with p and d. The equations of this latter set are 
used to obtain the coefficients of d in terms of 
p/d and the coefficients of p; then the coefficients 
of d are eliminated from the first set. 

As the first step in this procedure,” Eq. (111) 
was subtracted from (11j) with the result 


(M—K)[p(a2+a,) —d(as+4¢) | =0 
pa2—da,= — (pa,—das). 


or 
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By subtracting (11i) from (11k) we find that 
pa — dag = pas, — das. 
Obviously, then, 


paz — das = pas — da, =0 
from which 


and 


paz = dag (12) 
pas = das. (13) 


By use of the above relationships in (11j) and 
(11k) it is found through similar reasoning that 
pb. =dbg (14) 


and 
pb, = dbs. (15) 


In the same way, (11m) and (11n) give the 


relationships 


and 
pas = daz. (17) 


Equations (110) and (11p) are identically ful- 
filled. . 

The first set of eight equations was then 
simplified by substituting for the coefficients of d 
their values in terms of the coefficients of ». These 
equations become 


$(a1) +2p(a2+a3+44) (18a) 


$(a;) +2p(a2—a3—44) 
+K[s(b:)+2p(b2.—b3—by) ]=0, (18b) 


$(a;) +2p( —a2+a3— 44) 
+L[s(b1) +2p(—b2+b3—b4) ]=0, (18c) 


$(a1) +2p( —a2—a3+44) 
+ M[s(b:)+2p(—b2—bs3 +4) ]=0, (18d) 


d’ 
(b+ (9+) otbr+b)=0, (18e) 
pd’ 
-K[ 9a) +(e'+—) (2-01-00) 
pd’ 
+3/(b)+( 0+) o—-b.—2)=0, (18f ) 
pd’ 
-1 sa) +(9+=-)(-ar+a.-a)| 


pd’ 
+3/(0)+(2'+—)(-bi+b—b)=0, (18g) 


pd’ 
=| s"(a)+ (2+) (-a-a+09 
pd’ 
+3'(b)+(2'+—) (bade) = (), (18h) 


Up to here, no use has been made of the restric- 
tion K=L. These equations are identical with 
Eqs. (1a) to (ih) in Kimball’s paper'® on the 
diamond, with p and p’ of Kimball’s equations 
replaced, respectively, with 2p and p’+(pd’) /d. 
Kimball has obtained a solution of these equa- 
tions 
—(10K?M?+8K?+8KM+10M?+ 16) 


w= » (19) 
3K?M?+4K?—-4KM+3M 





where w = (u?+1)/u and p=sp’/s’p. This solution 
can be used for Eqs. (18) if u is redefined as 


_s(b"/p+d'/d) 
7 2s’ 


(20) 





yu 


Through use of Eq. (19), w can be evaluated 
directly from k, and since w is known through 
numerical computation as a function of E, the 
relationship of E and k can be determined. 
Determination of E from k involves the evalu- 
ation of s, s’, p, p’, d, and d’. These values were 
obtained through numerical integration of the 
Schrédinger equation; the details of this com- 
putation are given in Appendix I. A plot of w 
against uw gives a double-branch curve only the 
negative branch of which, as inspection of Eq. 
(19) will show, need be considered. As y ap- 
proaches 0 and —, w approaches —@ as a 
limit; this corresponds to K = M=0. The values 
of s, s’, p, p’, d, and d’ which will make |w| = 
form the limits of the bands at which k =0; these 
values are s=0, s’=0, p’/p+d’'/d=0, p=0, and 
d=0. The negative branch of w has a maximum 
of —2 when n»=-—1, which corresponds to the 
place in the band at which & has the highest 
value. At p= —1, s/s’+2/(p'/p+d’/d) =0. 
When Eq. (19) is used in the process of finding 
E asa function of k, it is found that to each value 
of k there are six values of Z. As shown in Fig. 4, 
at the observed half-internuclear distance the 
lowest of these values lies in the region between 
curve (I) s’=0, (k=0), and (II) the lowest 
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branch of s/s’+2/(p’/p+d’/d) =0, (k= ~); the 
second value is in the region bounded by this 
latter curve and (IV) the lower branch of p’/p 
+d’/d=0, (k=0); the region of the third value 
is bounded by (III) s=0, and by (V) the middle 
branch of s/s’+2/(p’/p+d’/d) =0, (k= @); the 
fourth value is in the region bounded by (V), 
and by (VI) p=0, (k=0); the regions of the 
fifth and sixth values are not shown in Fig. 4 
since they are of such high energy that they are 
of little interest. The region of the fifth value 
lies between the upper branch of p’/p+d’/d=0, 
(k=0), and-the highest of the three branches of 
s/s'+2/(p'/p+d’/d)=0 (k=); the region of 





ENERGY (RYDBERG UNITS) 





| i. 4 i 
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Fic. 4. Energy bands of silicon as a function of R 
(one-half the distance between nearest neighbors). The 
curves are: I, s’=0; II, s/s’+2/(p'/p+d’/d) =0; III, s=0; 
IV, (doubly degenerate) p’/p+d’/d=0; V, s/s’+-2/(p'/p 
+d’/d) =0; VI, (doubly degenerate) p=0. 


the sixth value is between this latter branch and 
the curve d=0, (k=0). 

In addition to the above six band regions, 
degenerate bands arise from Eqs. (3a) and (3b) 
in Kimball’s paper.!° These equations may be 
rewritten in the form 


2p(a2—d3) +K2p(b2—bs) =0, 
—K(p'+pd'/d) (a2—a3) + (p’ + pd’/d) (bz — bs) =0. 


They are satisfied by any value of k if p=0, or 
if p’+(pd’)/d=0. 

The values of the coefficients were determined 
for the case p=0, and are found to be 
a,=[—M/(2K—M) }(a2+as), (21a) 


b=[K/(2K—M)][(K—M)a:—Kas], (2b) 


bs=(—K/(2K—M)][Ka.—(K—M)as], (21c) 
bs =[(KM/(2K —M) ][a2+as]. (21d) 


For the condition p’/p+d’/d=0, the coefficients 
are 
a4= — (d2+4s3), (22a) 


b,=[1/(MK) ][ —Ka.+(M—K)as], (22b) 
bs=[1/(MK)][(M—K)a:—Kas],  (22c) 
bs=[1/K ][a2+<as]. (22d) 


Thus two bands of zero width follow p=0, and 
two follow each of the two branches of p’/p 
+d’ /d=0. 

Two additional bands of zero width can be 
shown to follow the curve d=0. To determine 
this it is necessary to go back to the original con- 
ditions, and to examine them for this special 
case. If d2=0, p0, it follows from relationships 
(11j) to (11p) derived from the continuity of the 
components of the tangential derivative that 
dz = 3 = a,=b,=b; =b,=0. By use of these values 
in Eqs. (11a) and (11b), we find that sa,;=0, and 
sb,=0; it follows that Eqs. (11c) to (11h) may 


be written as 
d'(bs+be+67)=0, (23a) 


— Kd'( —ds+d5—47) +d'( —bs+b.—b7) =0, (23b) 
— Kd'(—as—ag+a;)+d'(—bs;—bg+57) =0, (23c) 
— Md'(a,;—a¢—a7) +d’ (bs — bg —b;) =0. (23d) 


Through simple algebra it may be shown from 
the above that 


bs = —(bg+52), (24a) 
ds= —(1/K)(bs+57), (24b) 
dg= —(1/2)[(1/K —2/M)(be+57) 
—(1/K)(bs—57)], (24c) 
a, = — (1/2)[(1/K —2/M)(b6+7) 
+(1/K)(bs—b;)J. (24d) 


Thus it is possible to express all but two of the 
unknowns in terms of the others. From this it is 
seen that two bands of zero width are contained 
in the solution d=0. 

Certain essential features of silicon can be 
understood qualitatively through study of Fig. 4. 
One would expect that at absolute zero the eight 
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electrons of the unit cell, four electrons to the 
atom, at the observed half-internuclear distance 
would fill completely the two lower bands and 
the two degenerate bands following curve IV. 
With no partially empty bands, silicon would not 
conduct electricity at lower temperatures. 

That. silicon just misses being a metal can 
also be seen in Fig. 4. At _ half-internuclear 
distances greater than 3.5, the two degenerate 
bands following curve IV touch the empty 
bands bounded by curves IV and VI. If silicon 
could exist at these increased half-internuclear 
distances, it would be a metal of good conduc- 
tivity. 

To determine the extent by which the band 
structure obtained by including consideration of 
three 3d functions, besides giving a d band of 
high energy, differed from the band structure 
obtained by using Kimball's solutions!® in which 
no d functions were considered, the values of s, 
s’, p, and p’ calculated for silicon were used in 
Kimball’s solutions. The results are shown in 
Fig. 5 as a plot of energy against half-distance 
between nearest neighbors. Kimball’s solution 
p’ =0 (curve IV of Fig. 5) replaces p’/p+d’/d=0 
(curve IV of Fig. 4). The curve p’ =0 is signifi- 
cantly lower than the curve ~’/p+d’/d=0; at 
the observed half-internuclear distance, the 
energy difference of curves III and IV is 0.44 
Rydberg units more in Fig. 5 than it is in Fig. 4. 

Experimental information about the width of 
the M band of silicon is furnished by the soft 
x-ray emission spectra obtained for silicon by 
O’Bryan and Skinner.” The emission bands 
represent transitions from the broad M band toa 
relatively sharp L level, so that the width of the 
observed x-ray band is approximately equal to 
the energy width of the M band. The width of the 
x-ray band of silicon was found to be 19.2+2.5 
volts. The width of the band of Fig. 4 at the 
observed half-internuclear distance is about 16 
volts; that of Fig. 5 is about 9.5 volts. The 
photometer curve of the x-ray band is smooth 
with a gradual rise and a more abrupt fall, 
indicating more transitions of high energy elec- 
trons of the M shell than of low energy electrons. 
One might expect that there would be a well- 
defined large peak toward the high energy end 


17H. M. O'Bryan and H. W. B. Skinner, Phys. Rev. 45, 
370 (1934). 
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Fic. 5. Energy bands of silicon as a function of R 
(one-half the distance between nearest neighbors), using 
Kimball's solution for the diamond lattice. The curves are: 
I, s’=0; II, s/s’+p/p’ =0; III, s=0; IV (doubly degen- 
erate), p’=0; V, s/s'+p/p’=0; VI (doubly degenerate) 
p=0. 


of the x-ray band corresponding to transitions 
from the narrow bands (of zero width in our 
approximation) following curve IV of Fig. 4. It 
has been shown,'* however, that the optical transi- 
tion probability is not simply related to the level 
density in a given energy region. 

Unfortunately the results expressed in Fig. 4 
do not agree quantitatively with the conclusions 
drawn in the first part from the optical properties 
and the intrinsic conductivity. The optical 
properties lead us to expect an energy gap 
between the full band and the lower conduction 
band of not more than 1.5 to 2.0 volts, or about 
0.15 Rydberg unit; the conductivity measure- 
ments point in the same direction. The energy 
gap in Fig. 4, as estimated by extrapolation, is at 
least 0.8 Rydberg unit. 

It is interesting to compare this energy gap 
with that for the diamond and that obtained for 
silicon using Kimball’s solutions (Fig. 5). For 
diamond (2s and 2p orbits) Kimball'® has found 
a gap of 4 or 5 Rydberg units, which is also much 
too high, since diamond absorbs at about 1800A, 
and should have a gap of about 0.5 Rydberg 
unit. Figure 5, which takes only 3s and 3p orbits 
into account, gives an estimated energy gap for 
silicon of about 1.3 Rydberg units. The intro- 
duction of the 3d orbits brings this down to about 
0.8 unit. 


18H. Jones, N. F. Mott, and H. W. B. Skinner, Phys, 
Rev. 45, 379 (1934). 
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This disagreement of the calculated energy 
gaps with experiment is probably explained by 
Shockley’s statement!® that the method of Wigner 
and Seitz is fairly good for the ground state, i.e., 
the full bands—in Fig. 4 the region between 
curves I, II, and IV—but leads to completely 
erroneous results for excited states, i.e., the 
region between curves III, V, and VI. It must 
be expected that in reality curve III lies much 
lower on the left and is much flatter. 

The average electron energy for a given half- 
internuclear distance may be computed from 
Figs. 4 and 5 by means of?° 


Ew =0.625Er+0.250En+0.125E;. (25) 


The lowest average electron energy of the bands 
of Fig. 4 for the region between 2 and 3 units of 
length is found to be —0.55 for r=2.35; the 
observed half-distance between nearest neighbors 
is y= 2.22 atomic units. In Fig. 5, obtained from 
Kimball’s solutions, the point of lowest average 
electron energy as given by (25) is about —0.85 
Rydberg unit, occurring almost exactly at the 
observed half-internuclear distance of 2.22 atomic 
units. 

If we trust the general character of our results 
apart from the quantitative side, we can say 
that the absorption band consists of several 
overlapping bands. There is first the possibility 
of an absorption starting in the doubly degener- 
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Fic. 6. Schematic plot of absorption bands of silicon based 
on electronic structure as shown in Fig. 4. 


1*W. Shockley, Phys. Rev. 52, 866 (1937). 

*® This expression is the weighted average of the energies 
of the occupied bands and is obtained by taking one-fourth 
the sum of: the energy of the two degenerate bands 
following curve IV, 2E:y; the average energy of the band 
between curves II and IV, approximately }(E1v+£n); 
and the average ene of the band between curves II 
and I, aggreutamately tUlies-+- ik). 





ate level IV; since level IV contains all values 
of k, transitions can occur from IV to any level 
in the range III to VI. Since the electron level 
density is usually highest near the middle of an 
allowed band, and since VI is doubly degenerate, 
the absorption band starting in IV should have 
the following character: The long wave limit is 
due to a transition from IV to III; then, to 
shorter wave-lengths, follow transitions from IV 
to the interior of the band between III and V, 
with first increasing, then decreasing intensities, 
until V is reached ; the end state moves then into 
the band between V and VI, with a maximum 
intensity for an end state near the middle be- 
tween V and VI, whereupon the intensity de- 
creases again. At the short wave-length limit, 
there is almost a strong line, due to transitions 
between the degenerate levels IV and VI. This 
band, therefore, should have its long wave limit 
for transitions between IV and III, an intense 
short wave limit due to transitions between IV 
and VI, and two maxima in between. 

Four other bands in the same region, probably 
of much smaller absorption coefficient because of 
the spread of the filled levels, are brought about 
by electrons coming from the two bands between 
I and IV and going to the region between II] 
and VI. 

There should be one band with a long wave 
limit due to transitions between IV and III, 
which coincides with the long wave limit of the 
band described before, and a short wave limit 
due to transitions between II and V; imme- 
diately adjacent is a band with a long wave 
limit due to transitions between II and V and a 
short wave limit due to transitions between I and 
VI. This extends much farther toward the ultra- 
violet than does the band originating in the 
degenerate level IV. 

In addition to these, there should be a band 
with a short wave limit due to transitions be- 
tween II and V and with a long wave limit due 
to transitions between IV and VI. Finally, there 
should be a band where the long wave limit is 
due to a transition between II and V and the 
short wave limit to a transition between I and 
Ill. 

A schematic drawing of the absorption as 
brought about by these transitions is presented 
in Fig. 6. It is, of course, impossible to give 
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TABLE V. Values of Z, and of unnormalized wave functions P at various energies. 











Wave 

function 3s 3s 3s 3s 3s 3s 3p 3p 3p 3p 3d 3d 3d 3d 4f 4f 4f 
Negative 

energy 

(Rydberg 

—_, 0.0 0.2 0.3 0.6 0.9 1.2 0.0 0.2 0.3 0.6 | 0.0 0.2 0.3 0.6 0.0 0.3 0.6 

r p 

0.0 14.00 0.00 0 0.00 0.00 0.00 
0.01 13.52 0.87 0.01 0.00 0.00 
0.02 13.06 1.50 0.03 0.00 0.00 
0.03 12.62 1.94 0.07 0.00 0.00 
0.04 12.21 2.21 0.12 0.00 0.00 
0.05 11.83 2.35 0.18 0.00 0.00 
0.06 11.48 2.37 0.24 0.00 0.00 
0.07 11.15 2.31 0.31 0.00 0.00 
0.08 10.84 2.17 0.38 0.00 0.00 
0.09 10.56 1.98 0.45 0.01 0.00 
0.10 10.29 1.73 0.52 0.01 0.00 
0.12 9.79 1.15 0.66 0.01 0.00 
0.14 9.34 0.49 0.79 0.02 0.00 
0.16 8.92 -—0.20 0.92 0.02 0.00 
0.18 852 —0.88 1.03 0.03 0.00 
0.20 8.15 —1.52 1.13 0.04 0.00 
0.24 746 -2.66 -—2.66 -—2.66 -2.66 -—2.66 —2.67 1.29 0.06 0.00 
0.28 683 -3.53 -—3.54 -3.54 -3.54 -3.54 -—3.55 1.39 0.08 0.00 
0.32 6.26 -—4.14 -—4.14 —4.14 -4.15 -4.16 -—4.17 1.44 1.44 1.44 1.45 0.11 0.01 
0.36 5.76 -—448 -—4.49 —449 -—4.50 -—4.51 —4.53 1.44 1.44 1.45 145 0.14 0.01 
0.40 532 -460 -—461 -4.61 -—4.63 -4.65 -—4.67 1.41 1.41 1.42 1.43 0.17 0.01 
0.44 4.94 -4.52 -—4.53 -4.54 -—4.57 -4.59 —4.62 1.34 1.35 1.35 1.37 0.21 0.01 
0.48 4.60 -—4.27 —430 -4.31 -—434 -4.38 -—4.42 1.25 1.26 1.26 1.28 0.24 0.02 
0.52 431 -3.90 -—3.93 -3.94 -—3.99 -4.04 -—4.08 1.14 1.15 1.15 1.18 0.29 0.02 
0.56 405 -342 -346 -—3.47 -—3.53 -3.59 -—3.65 1.01 1.02 1.03 1.06 0.33 0.03 
0.60 3.83 -—2.86 -—2.90 —2.93 -—3.00 -—3.07 -—3.14 0.86 0.88 0.89 0.92 0.38 0.04 
0.68 3.47 —-1.59 —1.64 -—1.68 -—1.77 -—1.87 —1.95 0.55 0.57 0.59 0.63 048 048 048 049 0.06 0.06 0.06 
0.76 3.119 -—0.21 -—0.27 -0.32 -—043 -0.55 —0.66 0.21 0.24 0.26 0.31 0.59 0.59 0.60 0.61 0.08 0.08 [0.09 
0.84 2.97 1.19 1.11 1.05 0.94 0.80 0.67 -—0.113 -0.09 -—0.08 -—0.01 0.70 0.71 0.72 0.73 0.12 0.12 0.12 
0.92 2.79 2.53 2.45 2.38 2.26 2.11 1.98 -047 -043 -041 —0.34 083 0.84 0.85 0.87 0.16 0.16 0.16 
1.00 2.63 3.78 3.70 3.62 3.51 3.36 3.23 -0.80 -0.75 -0.73 -0.66 0.95 0.97 0.98 1.01 0.20 0.21 0.21 
1.08 2.49 4.90 4.83 4.75 4.66 4.51 4.39 —1.64 —0.9 0.26 0.27 ‘0.27 
1.16 2.36 5.88 5.82 5.75 5.68 5.56 546 -—-140 —1.36 -1.33 -—1.26 1.22 1.26 1.27 1.33 0.33 0.33 0.34 
1.32 2.13 7.36 7.37 7.33 7.36 7.32 7.30 -1.91 -—1.88 -1.86 -—1.79 1.50 1.56 1.59 1.68 049 O51 0.52 
148 1.93 8.24 8.36 8.36 8.54 8.64 8.77 —2.32 -—2.31 -2.30 -—2.26 1.79 1.88 1.92 2.06 0.69 0.73 0.76 
1.64 1.75 8.57 8.83 8.90 9.28 9.59 9.92 -—2.63 -2.65 -—2.66 -—2.67 2.09 2.21 2.28 248 0.95 1.01 1.06 
1.80 1.60 8.44 8.86 9.02 9.66 10.23 10.83 -—2.85 -—2.91 -—2.95 —3.02 2.39 2.57 2.66 2.95 1.27 1.36 1.44 
1.96 1.46 7.92 8.54 8.80 9.75 1064 11.58 -—2.98 -—3.11 -3.17 -—3.33 2.69 2.94 3.06 1.65 1.79 1.92 
2.12 1.34 7.09 7.93 8.31 9.62 10.89 12.25 -3.04 -—3.24 -3.34 -3.61 3.00 3.33 349 404 2.09 2.30 2.51 
2.28 1.24 6.04 7.11 7.62 9.32 11.04 1289 -3.04 -3.31 -346 -—3.86 3.32 3.74 3.96 2.62 2.92 3.22 
2.44 1.17 4.81 6.12 6.76 8.90 11.14 13.57 -2.98 -—3.35 -—3.54 -—4.11 3.64 4.18 446 542 3.22 3.65 4.09 
260 1.10 3.47 5.01 5.80 840 11.22 14.34 -—2.87 —3.35 —4.35 3.97 465 5.01 3.91 4.50 5.13 
2.76 1.05 2.05 3.82 4.75 7.84 11.31 15.23 -2.73 -—332 -3.64 -4.60 4.30 5.14 560 7.15 4.70 5.50 6.37 
2.92 1.02 0.59 2.58 3.64 11.45 16.30 -—2.55 -—3.26 —4.87 4.63 5.67 6.23 5.58 6.12 7.84 
3.08 100 -—0.88 1.30 2.50 6.65 11.64 —2.33 -—3.18 -—3.65 -—5.14 4.96 6.22 6.91 9.35 6.56 7.98 9.58 
3.24 100 —2.33 0.00 1.33 11.91 —2.10 —3.09 —5.44 5.29 6.80 7.65 9.50 11.61 
340 100 -—3.75 -—1.29. 0.15 5.41 12.26 —1.84 -2.97 -—3.61 —5.77 562 741 8.43 12.13 11.21 13.98 
3.56 100 -—5.11 -—2.57 —1.03 12.71 —-1.56 —2.84 —6.12 5.93 805 9.26 13.14 16.73 
3.72 100 -640 -—3.83 4.18 13.27 —-1.26 -—2.69 -—3.53 -—6.50 6.23 8.70 10.14 15.58° 
4.04 100 -—8.70 -—6.24 2.97 14.77 —0.63 -2.36 -—3.42 -—7.37 6.78 10.08 12.06 19.82 
4.36 1.00 —8.46 1.80 16.89 0.02 —-1.99 -—3.28 -—8.41 11.53 14.19 25.03 
4.68 1.00 —10.46 0.65 19.78 0.67 -—1.58 -—3.12 -—9.67 13.06 16.55 31.41 
5.00 1,00 —12.21 —0.49 23.63 1.30 -—1.15 -—2.95 —11.19 14:64 19.16 39.22 








estimates of the relative strengths of these bands 
without extensive calculations.'* 

The experimental data (Fig. 1) seem to show 
that the absorption band is complex, but there 
is no sign of an absorption line due to transitions 
between degenerate states IV and VI and no sign 
of the long ultraviolet “wing” that should come 
from transitions originating in the band between 
I and II and ending in the band between V and 
VI. 

The writer wishes to acknowledge his in- 
debtedness to Professor K. F. Herzfeld who sug- 
gested the problem and the method of solution 
and who gave much valuable aid in the work. 





He also wishes to thank Dr. G. E. Kimball of 
Columbia University and Assistant Professors 
C. A. Beck and F. E. Fox of Catholic University 
for suggestions and discussion. 


APPENDIX 


The field used in the numerical integration of the 
Schrédinger equation is a self-consistent field in which no 
account is taken of the correlations giving rise to exchange 
energy. Since a modification of Hartree’s method of ob- 
taining solutions of the Schrédinger equation for the free 
atom was used in calculating the solutions for the solid™ 


21 An account of this method together with a discussion 
of the error introduced by neglect of exchange may be 
found in Chapter IX, F. Seitz, The Modern Theory of 
Solids (McGraw-Hill Book Company, Inc., New York, 1940). 
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the field of the core functions (with small changes) ob- 
tained by McDougall through the simple Hartree self- 
consistent field method was used in preference to the field 
calculated more recently, in which exchange terms were 
considered. The core functions were estimated for neutral 
silicon by comparing the core functions obtained by 
McDougall for Si V with those obtained by Donley™ for 
Si IV and Si III and then applying small corrections. The 
extrapolated values differed only slightly from the original 
values; the error introduced by extrapolation is believed 
to be far less than that brought in by other approximations. 
To this core field was added three-quarters the field of a 
combination of three 3p electrons and one 3s electron 
normalized to 3.18 atomic units, the radius of the sphere 


2 J. McDougall, Proc. Roy. Soc. A138, 550 (1932). 

*W. Hartree, D. R. Hartree, and M. F. Manning, 
Phys. Rev. 60, 857 (1941). It is noted here that the optical 
term values calculated for SilV and Si V from the field 
with exchange do not agree as well with experiment as 
those obtained by McDougall (reference 22) using the 
wave functions found by methods neglecting exchange in 
the wave functions but including exchange terms in the 
expression for energy. The authors conclude that neglect 
of exchange has two counteracting effects which may 
balance each other. 

*H.L. Donley, Phys. Rev. 50, 1012 (1936). 


having the same volume as the cellular polyhedron of 
silicon. In consideration of the symmetry of the silicon 
lattice, the use of this spherically symmetrical valence 
electron field (three-quarters the field of three 3p electrons 
and one 3s electron) was thought to be preferable to use of 
a non-symmetrical field of three valence electrons (the 
field of two 3s electrons and one 3 electron, or of one 3s 
and two 3p electrons), as would be used in Hartree’s 
procedure for obtaining wave functions for the free atom. 
Successive approximations were made with energy values 
for the 3s and 3p electrons approximately equal to those of 
the free atom until the expected difference between initial 
and final fields was not greater than about 0.01 at any 
point. It was considered needless to seek greater accuracy 
in the field than this since the other approximations 
required in building up the wave function of the solid 
introduce unavoidable error far greater than that brought 
in by a field of this accuracy. 

In Table V is presented, for various values of r, a tabu- 
lation of Z,, the effective nuclear charge for potential, i.e., 
the point charge which would give, if it were placed at the 
nucleus, the same potential at radius r as that of the actual 
field. Also presented in Table V are the unnormalized 
functions P for various energies; P=r times the radial 
part of the solution of the Schrédinger equation. 














PHYSICAL REVIEW VOLUME 66, 


NUMBERS 11 


AND 12 DECEMBER 1 AND 15, 1944 


On Milne’s Theory of Gravitation 


A. SCHILD 
Department of Mathematics, University of Toronto, Toronto, Canada 


(Received September 29, 1944) 


The present paper is a criticism of Milne’s theory of the gravitational two-body problem. In 
particular, it is shown that Milne’s theory is not invariant under a Lorentz transformation and 
therefore, since it is based on “kinematical relativity,” is inconsistent. 





1 


N a series of papers,' Milne claims to have 
obtained an ‘‘expression of the inverse square 
law, or its potential 1/r, in a form relativistically 
invariant under Lorentz-transformation from any 
one fundamental observer to any equivalent 
observer.”? The purpose of this paper is to show 
that Milne’s theory of gravitation, contrary to 
this and other similar assertions, is not invariant 
under a Lorentz-transformation. 

Milne, in his writing, deliberately foregoes 
“any mention of the construct ‘space-time.””® 
However, the four-dimensional tensor-notation 
not only enables us to rewrite formulae in a 
form more compact than is possible in the 
inappropriate three-dimensional vector calculus; 
it also contributes to greater clarity. The co- 
variance of a relation under Lorentz transforma- 
tions is easily tested if the four-dimensional 
tensor calculus is used consistently. 

We denote the coordinates of an event by x‘ 
(¢=0, 1, 2, 3) or alternatively, when the occasion 
arises, by (¢, x, y, z). The usual summation-con- 
vention applies to double indices. Since we are 
concerned only with special relativity, our metric 
form is that of Minkowski space: 


ds? = nidx'dx!,  nig= 05 
= dt? —dx* —dy* —dz’. 


(The velocity of light ¢ is taken equal to 1.) 


(1) 


In order to avoid confusion, it is well to dis- 
tinguish clearly between a particle and an event. 
Much of the criticism in this paper is caused by 
Milne’s failure to differentiate sufficiently be- 

1E. A. Milne, I, II, III, Proc. Roy. Soc. A156, 62-85 
(1936) ; 160, 1-23, 24-36 (1937) ; and IV, Phil. Mag. [VII] 
34, 82-101 (1943) ; hereafter referred to as I, Il, II, IV, 
respectively. 

2 IV, Sec. 35, p. 96. 

SIV, Sec. 13, p. 84. 


tween these concepts. For our purpose, we may 
define a particle simply as a world line,‘ i.e., a 
curve x‘=x‘(@) in space-time, @ being some 
arbitrary parameter along the curve (usually the 
“propertime” s or the coordinate ¢). With a 
particle or world line may be associated physical 
scalars such as the rest mass m, the electric 
charge e, etc. An event, of course, is a point com- 
pletely determined by its four coordinates x‘. 


2 


The fundamental problem which Newton's 
theory of gravitation solves in principle is the 
following: Given initial positions and velocities 
of an isolated system of gravitating particles, 
i.e., given the points of intersection of their 
world lines with the plane ¢=const. and the 
directions of the world lines at these points, to 
find the world lines. Newton’s theory enables us 
to write down a second-order differential equation 
for each world line; the unique solution of the 
system of equations with the given initial condi- 
tions can then be shown to exist. 

As is well known, the classical theory of gravi- 
tation, like all theories involving action-at-a- 
distance force laws, is not “Lorentzinvariant.” 
This is understood to mean that the differential 
equation of the world line of a particle is changed 
in form under Lorentz transformations; alter- 
natively, if observers in different Lorentz frames 
use such a.theory to calculate the world line of a 
particle, they will in general obtain different 
curves. This deficiency of the Newtonian law of 
attraction is immediately traced to the fact that 
in the theory of relativity simultaneity of distant 
events is not an absolute concept, but is meaning- 
ful only in a fixed inertial frame. Precisely the 

4In order to adhere to common usage, the expression 
“world line of a particle’”’ will be utilized in the following, 


although our definition permits us to replace it by the one 
word “particle.” 
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MILNE’S THEORY OF GRAVITATION 


same argument also applies to Milne’s formula- 
tion of the inverse square law. 


3 


First, however, mention will be made of a 
minor difficulty in Milne’s analysis. The coordi- 
nates x‘of an event are in general not components 
of a vector ; only the coordinate-distance between 
two events x‘ and x)‘, i.e., x'—x@', is a vector. 
Milne considers only Lorentz frames (his ‘‘funda- 
mental observers’) with a common origin O, 
chosen at the point in space-time where his 
model of the expanding universe is initially con- 
centrated, i.e., the common point of intersection 
of all world lines. Thus, limiting himself to 
Lorentz transformations which correspond to 
rotations only of the coordinate frame, he can 
consistently put x~@)'=0=(0, 0, 0,0) and hence 
treat x‘ as a contravariant vector. In the follow- 
ing we shall retain this particular aspect of 
Milne’s cosmological model and limit ourselves to 
Lorentz frames with the common “natural 
origin.” But it is well to emphasize that this 
aspect of Milne’s cosmology (which is common 
also to some other cosmological theories) is a 
prerequisite to his formulation of the inverse 
square law;® if the former is repudiated, the 
latter loses all significance. 


4 


A system of two particles under mutual gravi- 
tational attraction is considered, the rest masses 
being mq) and m.). Let us assume that the world 
lines of these particles have been obtained and 
investigate whether Milne’s differential equations 
of motion can be satisfied by the world lines 
irrespective of the Lorentz frame. 

Let Pa) =x‘, Pe) =X‘ be events on the two 
world lines. We construct a “potential” x : 


M (1) (2) 


Mo 





x=—- 


X(1y jX (2)? 


_ (2) 


. ’ 
. 
{ (2¢¢1y aX 2")? — X24 c2y* + Xr XK ary"} # 





where M, is a constant which plays the role of 
the constant of gravitation and x;=7,,x*, as is 
suggested by the notation. Milne’s equations of 

5Also to Milne’s theory of rational electrodynamics, 


Proc. Roy. Soc. A165, 313-332, 333-357 (1938); arid Phil. 
Mag. [VII] 34, 73-101, 197-211, 235-258 (1943). 
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motion,® written in tensor notation, are 


d dx1)* 
——{ maka! ) 















































ds) $(1) 
Mayéqay! dx" dx‘ 
5 rca” X1y*—Xcye° . 
X(1ykX (1) dsqay dsqy 
dx¢1)* d , 
- -——(m &y*) = n# (3) 
dsqa) dsq) - Xa)? 
d dx 2)* 
——{_ m2) & 24 
ds 2) aS (2) 
M2) &(2)* dx 2)" dx,2)* 
a ener X(2)*— X(2)r° : 
X(2)k%X (2) Say ds(2) 
dx,2)* d re 
- ; (m2) Ec2)*) = n*# (4) 
ds (2) dS (2) OX (2)? 
where’ . 
( =) ( dx (2)/\? 
X(1)j %(2) 5 
ds. ds (2) 
Ey = : and &)= — 
X1yeX 1) X (2) kX (2) 


If Eqs. (3) and (4) together with initial 
conditions (see Section 2) are to yield solutions 
Xa) *=Xy (Sy), Xa) *=X (Sa), a relationship be- 
tween the two parameters sq) and sg must be 
postulated, to be applied to the partial deriva- 
tives of x in the equations of motion. To put it 
differently, given a value of sq, i.e., given an 
event Py on the first world line, the differential 
equation (3) tells us about the behavior of the 
world line at Pq only if we know how to deter- 
mine 0x/0%X)* which depends not only on Pa)(sq) 
but also on an event P,2)(sq)) on the second world 
line. We must have a rule for associating an 
event Py(sq)) with an event Py)(sq), ie., a 
relation between sq) and sq. A corresponding 
argument applies to Eq. (4). It should be men- 
tioned, however, that the relations between 
Sa) and Sq) need not in general be symmetric 
and that we obtain the relation which is to be 
applied to Eq. (4) by interchanging suffixes ‘*(1)"’ 
and ‘‘(2)”’ in that applied to Eq. (3). 

Equations (3) and (4) are apparently in 
tensor form, but they are in general true tensorial 
equations if and only if x behaves like a scalar 
under Lorentz transformations. This is the case 
~ STI, Sec. 11, p. 12, Eqs. (45), (45’), (46), (46’). 


7 ds 1) and ds 2) are given by (1), i.e., 45 (1)? = 95 dX (1)*dx,1)! ; 
ds 9)? = 9; jdX(2)*dx(2)/, 
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if in Eq. (3) [or (4)] an event x)‘ (or xq‘) is 
associated with each event xq‘ (or x@*) in a 
manner independent of the frame of reference. 
Otherwise Milne’s equations are not invariant 
. under Lorentz transformations. Milne uses simul- 
taneity*—t) =¢)—to set up a symmetric corre- 
spondence between events on the two world 
lines. This correspondence does, of course, depend 
on the particular frame of reference. If one ob- 
server finds that the two world lines satisfy 
Eqs. (3) and (4) with ¢) =¢,), a second observer 
may fail to do so. Thus the theory is not invariant 
under a Lorentz transformation. 

Repeating, we may say that Eqs. (3) and (4) 
alone are invariant but meaningless as equations 
of motion. By adding the relation of simultaneity, 
the equations of motion become meaningful but 
lose their Lorentz invariance. Milne anticipates 
this in part when he says, ‘‘Of course the interpre- 
tative relation fq) =?) =¢ is not invariant,’’® but 
he fails to draw the conclusion that the invariant 
character of his whole gravitational theory is 
affected. 

It is not difficult to construct a correspondence 
which is invariant. As an example’® we consider 
the light cone with xq) ‘=P as vertex and with 
its axis pointing to the past. Its equation in 
any Lorentz frame is (x‘— xq) *)(x;— xi) =0. We 
now define Py) corresponding to Pa) as the 
point of intersection of the cone and the second 
world line. Interchanging suffixes ‘‘(1)”’and“*(2),” 
we obtain the relation to be applied to 0x/dx,)‘ 
in Eq. (4). This is the well-known method of the 
retarded potential. However, the advanced po- 
tential or any linear combination of the two 
could be used equally well. Thus, to make the 
theory relativistically invariant, new and to 
some extent arbitrary assumptions have to be 
added which destroy the epistomological char- 
acter of Milne’s theory. 


5 


There is an exceptional case in which one of 
the two equations of motion is invariant under 
Lorentz transformations. If the world line of 
the attracting particle, say the second, is a 


8 II, Sec. 11, p. 13. 
* II, Sec. 9, (1), p. 7. 
1°Compare Milne, Proc. Roy. Soc. A165, Sec. 12, 340 


(1938). 
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straight line through the common origin O, i.e., 
if its equation is x~@)'=A*‘s),!! where A‘ is an 
arbitrary unit vector, then 


X19? 
Mo =} (XqA*)?—X~paxy*}? 


and x is independent of s,2) and hence of xq)‘. We 
then have a one-body problem with the equation 
of motion (3), and the arguments of the previous 
paragraph do not apply. In fact, we have the 
case of one particle moving in a given external 
field of force. In all applications of his theory, 
this is the case which Milne assumes, although 
no explicit mention is made of it. A coordinate 
system is chosen in which the attracting particle 
is at rest, i.e., such that its world line coincides 
with the ¢ axis; this of course implies that the 
world line is a straight line through the origin." 

In Section 15 of II, Milne generalizes his theory 
to a system of n particles. Again, the theory is 
relativistically invariant only if the world lines 
of m—1 particles are straight lines through the 
origin, and then only the equation of motion of 
the last particle is given. 

That Milne does not himself wish to limit his 
theory to the one-body problem discussed in 
this section can be gathered from the following 
quotation : 

“The core of the present paper is . . . to ob- 
tain a general formula for the potential energy x 
of any two massive particles mq) and mq)... . 
It is applied to obtain the relativistic equations 
of motion ... of any two particles in one 
another’s presence.’’?* 


mM (1)M (2) rr 1 (5) 





© 


1 It should be noticed that Eqs. (3) and (4), which are 
Milne’s Eqs. (45), (45’), (46), (46’) of II, Sec. 11, p. 12 are 
not both satisfied by x;2)'=A*sc2)(m 1), m2), finite) ; since in 


this case £2) =1, a Xi, and Eq. (4) (on multiplication 





(2) 


by nj) yields 








9x _ 9 _ My M2) X¢1yaX1y* 
Oxy9 Myf (capi cay*)?—XcpiXcy* -Xcmexcy*}! 
» ((1)1%(2)*% (2) j — X(2yiXc2y*Xy4) 
_ MaMa) | Xcyaxcry* 





Mo = {(Xcpid*)?— xy iXcy*}! 
1 : 
*— + (XAG —X(y))- 
5(2) 


Hence x,1)iA*A; = X(1)j, Xj iS parallel to A;, and thus finally 
%(1)j =AjS). This shows that the two particles coincide, x 
becomes infinite and the problem of gravitation loses all 
significance. 
12 [], Sec. 17, p. 15; and IV, Secs. 33, 34, p. 95; the argu- 
ment of Sec. 34 only holds in the case we are considering. 
13], Sec. 2, p. 2; italics are my own. 
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Sec. I. This paper inverts certain quantum mechanical 
ideas by constructing a “probability of presence’’ from first 
principles. If x;=/;(s,a:,@2) represents a 2-parameter 
family of curves in 3-dimensional space, the density of the 
functions f; with respect to the parameters @:@2 is defined 
by means of an expression for the normal cross section of a 
sheaf of curves. This density plays the role of a probability. 
Its derivative along a curve of the family is independent of 
the parameters @,@2. This circumstance permits a special 
density to be constructed, called the ‘“‘characteristic” or 
“smoothest” density, which is independent of the parame- 
ters and which can be found by solving a linear second- 
order partial differential equation. The extension to 
n-dimensional space and (m — 1) parameters is given. Sec. IT. 
The concept of orbital density is the dynamical equivalent 
of the density of Sec. I. Dividing the orbital density by the 
velocity of the representative point gives the quantity we call 
kinematical probability P, which satisfies the fundamental 
equation (2.8), and which has the same force as any other 
probability of presence. Equation (2.8) isused toderive from 
the Hamilton-Jacobi equationa “wave” equation for »/P for 
a system of m planets moving around a sun under the condi- 
tion that the space average of the mutual interactions of 
the planets is a minimum. This provides a new approach 
to planetary dynamics and yields an equation that is very 


much like Schroedinger's equation for the stationary states 
of nm particles, and has the same basic interpretation. 
Sec. III. A non-conservative system of n particles is as- 
sumed to be in a state such that the action function can be 
written W=—Et+R(x;i, t) in which R still contains the 
time explicitly. Starting with the Hamilton-Jacobi equa- 
tion, using Eq. (2.8), and imposing the condition that the 
space average of 9R/dt is a minimum, an equation 
{viz., (3.10)] is derived which differs from Schroedinger's 
equation for stationary states only in certain respects that 
cannot be regarded as fundamental. Despite this basic 
equivalence of Schroedinger’s equation and our equation 
(3.10), conflicts in interpretation and formalism arise which 
are briefly discussed. Sec. 1V. Complex U’s are introduced 
into Eq. (3.10) in a simple way which fails to bridge the 
gap between our formalism and that of quantum mechanics. 
The kinematical probability defined by Eq. (3.10) is shown 
to have the “smoothest” property. To account physically 
for the stationary states an electromagnetic wave system 
is proposed which has the property that a particle moving 
in it may be in “phase-mechanical” resonance with it. 
Macroscopic examples of this kind of resonance are given. 
This concept represents an inversion of physical ideas of 
quantum mechanics corresponding to the mathematical 
inversion with which we began. 





I. INTRODUCTION 


HE mathematical theory presented here is 
an inversion of quantum-mechanical ideas 
in the sense that instead of determining a proba- 
bility density or probability ‘tof presence” by 
solving a prescribed differential equation, we 
construct such a probability ab initio from ele- 
mentary geometric or mechanical concepts and 
determine some of its properties. The first results 
evolving from this construct were presented at a 
recent meeting of the American Physical Society ;! 
here they will be generalized and amplified. 
Let the solutions of the system of differential 
equations 


dx dx dxs 
anu@bemetiae (1.1) 
X,; Xe Xs 
be of the form 
xi=fi(S, di, G2). (1.2) 


Here the X,, as functions of the x;, are single- 
1B. Liebowitz, Phys. Rev. 65, 255 (1944). 
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valued or else confined to one branch; if the X; 
contain any parameters };, they will be supposed 
fixed. Then Eq. (1.2) defines a 2-parameter sys- 
tem of curves in a 3-dimensional space through 
each point of which passes one and only one 
curve (except at singular points). The inde- 
pendent variable s in Eq. (1.2) measures distance 
along any such curve. These equations determine 
a system of unit tangent vectors t given by 


dx; 
t=—=X(LX/) (1.3) 
ds 


which may be expressed as functions of the x, or 
of s, a1, Gs. 
By means of the displacement 


Ox, OXe OX; 
és= (45+ Jao 
Oa, Oa; da; 


os Os 
=—46a; (definition of —), 
0a, 0a 
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in which a; alone is varied, we pass from a given 
point on a given curve to a nearby point on a 
neighboring curve. The normal distance between 
the two curves is 


6ni= (=x t)éar, 


where the components of t are given by Eq. (1.3). 
Likewise if a2 alone is varied, the normal distance 
to a different neighboring curve, assuming a; and 
a2 are independent, is given by 


sns= (—xt) ans 


The vector 6n;X dn, is tangential to the curve; 
its magnitude represents the normal area pierced 
by the sheaf of curves comprised between 
x;(S, @1, @2) and x;(s, @1+6a1, d2+5a2). The den- 
sity of the curves is obviously inversely pro- 
portional to this area. Hence we define the density 
Q of the functions f; with respect to the parame- 
ters @:@2 by the equation 


t _ 6m Xéns 
-(—x xt) x (=xt) (1.4) 
QO. “Qbdee 0a, 
ds OS 
= —-—xt)t 
0a, Oa2 


or 
ty te ts 


Ox; OX, OX; 
—_ —_—_ |. (1.5) 


1 
Q 0a, 0a; 0a; 


OX, OX2 OX; 


Oa Ode Ode 








The determinant is to be taken positive. It 
depends obviously on the parameters @;d2, but 
the derivative of log Q in the direction of the 
curve does not. For it is not difficult to prove that 

om -t, (1.6) 
Q ds 


and since the components of t can be obtained by 
Eq. (1.3) directly from the differential equation 
(1.1), it follows that the derivative of the loga- 
rithm of the density Q taken along a curve (1.2) 
is independent of the choice of the parameters 
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@1a2. The derivative in any other direction ts de- 
pendent on the a’s. 

The density Q has the force of a probability of 
presence. It will be clear from the definition that 
if a point 0 is constrained to lie on a curve of the 
set defined by Eq. (1.2) and if the a’s are chosen 
so that all values of a; and a2 are equally probable, 
then the probability of finding the point 0 in the 
volume element dxdydz is proportional to Qdxdydz. 
If it is desired to know the average value of some 
property of the functions (1.2), e.g., the square 
of the gradient of the density Q, it is only neces- 
sary to form the integral /(VQ)*Qdr/ fQdr over 
the domain of interest. In general, the density Q 
will depend on the choice of a; and az, but it is 
possible to define a density Q. which will be called 
the ‘‘smoothest” density or ‘“‘characteristic”’ 
density, which has special properties and which 
does not depend on the a’s, but does depend on 
certain boundary conditions, as will be shown. 

If the derivative of Q in all directions perpen- 
dicular to t should vanish, then (1/Q)(dQ/ds) 
would represent the gradient of log Q. Obviously 
this would be a very special case. More generally, 
[(1/Q)(dQ/ds) ? represents the smallest value 
which the square of the gradient of log Q can 
have. Let G? represent the difference between the 
square of the actual gradient of log Q and the 
square of this minimum value. Then the charac- 
teristic density is defined as that density which, 
on the average, makes this difference a minimum. 
In symbols 


d 2 
G*=(V log Q)?— (— log ) = (Vlog Q)?—(V-t)?, 


5 f G*0.dr=0. 


Putting Q.= ¢*, we find the variational equa- 


tion is 
if |x(— log °) -1(0-t)!| pdr=0, 


which leads, by the usual rules, to the necessary 


condition 

09 V-t\? 

Bry eno, an 
Ox,? 2 


subject to appropriate boundary conditions. 
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This result depends, of course, on the existence 
of the integral to be minimized, but not on the 
prenormalization of ¢* because g can be multi- 
plied by an arbitrary constant before the varia- 
tion and/or after without altering anything. The 
equivalent of Eq. (1.6) isgiven for twodimensions 
in reference 1; Eq. (1.7) was presented at the 
meeting cited but is not printed in the reference. 

The generalization to an (m—1)-parameter 
family in n-dimensional space is immediate. For 
example, the generalization of the definition 
Eq. (1.5) for the density Q is 




















ty te cee f. 
Ox, OXe OXn 
0a, 0a, 0a, 
1 Ox, OX OXn 
—=| — eee ss (1.8) 
Q 0a2 das 0a2 
Ox) OXs OXn 
OGn-1 OAn-1 OAn-1 


Equations (1.6) and (1.7) remain formally un- 
changed and so does the meaning of Q as a 
probability of presence. The corresponding sys- 
tem of differential equations is now 


dx,/X,=dx2/Xo= +++ =dXn/Xn. (1.9) 


These concepts, definitions, and equations 
represent a definite step toward the construction 
of a “calculus of finite precision,” which, it is 
believed, is of intrinsic interest for both mathe- 
matics and physics. We shall not undertake here 
to illustrate the application of this calculus to 
purely mathematical matters; in the next section 
we shall present an application to dynamics. 


Il. APPLICATION TO DYNAMICS 


We consider a dynamical system of m particles 
moving in a field of force determined by a 
potential V which may contain the time ex- 
plicitly. The action function W, which is the 
solution of the Hamilton-Jacobi equation 


OW 1a 1 sdWy? 
—+-$—(—) +7=0, (24) 
Ot = 2 pms MeN OX 


will in general contain, beside the 3m coordinates 





x;, 3n parameters c;, and the time ¢. We suppose 
that the parameters c; are fixed, and that the 
time is frozen at some instant t. Then the 
3n-dimensional configuration space will be filled 
with a frozen system of action surfaces 


W(x, Cj, to) =C (2.2) 


and of orbits normal thereto. 
The momenta 
pr=dW/dx, (2.3) 


determine a system of unit tangents to the 
orbits, viz., 


Pi aw aw\?\-4 
i=—=)D; 2)\-} — —_ —— , 2.4 
m= (EBM —{7—) } (2.4) 


so that the equations of the orbits may be 
written 


dx /p1=dx2/p2= oo? =dX3n/Psn- (2.5) 


The only difference between Eqs. (1.9) and 
(2.5) lies in the circumstance that the X; in 
(1.9) are given functions of the x; and the 
parameters b, (supposed fixed), whereas the /; of 
(2.5) are unknown but nevertheless definite 
functions of the x; and the parameters ¢ (also 
supposed fixed). The apparatus of the calculus 
of finite precision described in the first section is 
therefore applicable; the only change we make is 
one of nomenclature—we call the density of 
solutions of Eq. (2.5) the orbital density, but 
retain the same symbol Q. Using the values of ¢; 
given by Eq. (2.4) we have 


Ot; 0 pi 
v-t=D—=L—— 
Ox; OX; p 
(2.6) 
-(v dp 
=_— ‘ees > 
p ds 


so that Eq. (1.6) for the orbital density be- 
comes 


V-p 
P 
Now the kinematical probability P, or the 

probability of presence of the representative 

point in the volume element dx;---dxs,, is di- 


rectly proportional to the orbital density Q and 
inversely proportional to the velocity (>> (#;,)*)! 


(2.7) 


=f 


qd @Q 
ds p 
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= (>> (p;/mz;)*)*. For the present, let it be assumed 
that all the masses are equal; then P=const. 
XQ/p and Eq. (2.7) becomes 


d V-p 
—— log P=—. (2.8) 
ds p 


This equation is given in reference 1 for one 
particle. It can be derived for n particles by the 
method there used, but the present derivation is 
more general and reveals the dependence of 
grad P on the choice of parameters like the a; 
of Eq. (1.2). The restriction to equal masses is 
not essential; it can be removed at any time by 
going over to a space £;=x,/mj. 

Equation (2.8) is the basic theorem of our 
dynamics. For subsequent applications it is 
useful to transform it by first using p>=dW/ds, so 
that it becomes 


clog P V-p 


dW op? 





and then transforming W by writing? 
W=+h log U. (2.9) 
Equation (2.8) thus becomes 
d log P V-p 


o%)—. (2.10) 
d log U p? 





Our first application is to the following classical 
problem : Given m planets moving about a sun, 
under what condition will the average of the 
mutual interdction potential of the planets be a 
minimum? Let the masses be all equal (a non- 
essential restriction, as mentioned above) ; using 
an obvious notation, write V; for —G>> Mm,/r; 
and V2; for —G>-mm;,/r;;; then, since the system 
is conservative, the action function W may be 
written as —Et+S(x;) and the H.J. equation 
becomes 


~ § (— ) -e- Vi)=-V 
2m i= OX, 


Now introduce S=+K log U, multiply by U’, 
and integrate over the 3n-dimensional space; 


2Compare Schroedinger Abhandlungen zur Wellenme- 
chanik (J. A. Barth, Leipzig, 1928), p. 1 


the result is 


Sle) —pe-ror|e 


~ Re V2U"dr. 


If U? were identified with the kinematical proba- 
bility P, the right side of this equation would be 
equal to the average value of the mutual inter- 
action potential of the planets, and equating its 
variation to zero would furnish the answer 
sought. Assume this to be the case; then, subject 
to appropriate boundary conditions, the varia- 
tional equation leads, by the usual rules, to the 
necessary condition 


T+ U=0, (2.11) 


where ~;?=2m(E—V;). Using U=exp (+5S/K) 
and p;=0S/dx;, we find this condition becomes 


p?+p;? 





V-p=F (2.12) 


But if, as in the case of the solar system, the total 
momentum of the mutual perturbations of the 
planets is relatively very small, p? and p;? are 
nearly equal, and Eq. (2.12) becomes V-p 
= +2p?/K. When this is inserted in Eq. (2.10) 
(with K for k), the result is P=const.X U?. 
Thus the assumption about U is verified to the 
extent that the momentum owing to mutual per- 
turbations of the planets is small compared with 
the total momentum of the system. 

Equation (2.8) and its equivalent (2.10) are, 
I believe, entirely rigorous. The only lack of 
rigor in our derivation of Eq. (2.11) lies in the 
sufficiency hiatus, but it is unlikely that physi- 
cists will cavil about that. (The momentum 
approximation does not constitute a Jogical weak- 
ness.) There are two noteworthy features of these 
results: (1) a new approach has been found to 
the basic classical problem of planetary motion 
(cf. Bode’s law), and (2) an equation which is 
the same in form and meaning as Schroedinger's 
has been derived by substantially rigorous means, 
purely classically for a purely classical problem. 
The constant K can be determined, as soon as U 
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and a single component of momentum are known, 
from pj=+(K/U)(dU/dx;), subject, however, to 
certain characteristics of momenta whose diver- 
gence is zero, as will be brought out in the next 
section. 


Ill. APPLICATION TO ATOMIC DYNAMICS 


It is surprising how quickly the construct we 
have called kinematical probability leads to 
formalisms like those of quantum mechanics. 
Starting with the basic theorem (Eq. 2.8) 


d 


V-p 
—— log P=—-,, (3.1) 
ds p 


and transforming it by means of Eq. (2.9), 
W= xk log U, (3.2) 
the result as found in Section II (Eq. 2.10) is 
d log P V-p 


= -—Rk. (3.3) 
d log U p? 





Now if either of the two following equations is 
assumed, 


P=const. X U?, (3.4) 
op? 
V-p= t (3.5) 


the other follows. Also if Eq. (3.5) is entered with 
pi=0W/dx;= +(k/U)(AU/dx,), it becomes 


aU p? 
+—U=0. (3.6) 





p & 


Ox,;* 


Furthermore, the formula for the average of a 
momentum component becomes 


dU 
pix { piU%dr= 4k U—dr; 
Ox; 


likewise corresponding formulae for 7, etc. are 
obtained, provided that after partial integration, 
the integrated part UdU/dx; can be assumed to 
vanish. 

Equation (3.6) is just the same in form and 
meaning as Schroedinger’s and the formulae for 
the momentum averages are very similar to 
those of quantum theory. There are differences, 
to be sure, which will be discussed below. 
And there is conflict! Thus, where we write 


pi= +(k/U)(0U/dx;), quantum formalism writes 
p;U=(h/t)(8U/dx;) which has only operational 
significance and in which +/ —1 enters in place of 
our +1. Before discussing this conflicting paral- 
lellism we shall show that Eq. (3.6) can be derived 
straight from the H. J. equation by means of our 
basic theorem without assuming either (3.4) 
or (3.5). 

Except in the planetary problem, all our equa- 
tions are applicable to non-conservative systems, 
i.e., systems in which the potential energy con- 
tains the time explicitly. Let us assume, however, 
that the system under consideration, though non- 
conservative, is nevertheless in a ‘‘stationary 
state,” by which we mean that the system is in 
a quasi-conservative state such that it hovers 
around a definite energy E. Corresponding to 
this we may write 


W=—Et+R(x;,, t), (3.7) 
which differs from the usual form in that R still 


contains the time. The H. J. equation now reads 
(with all masses equal) 


1 oR\? OR 
nt) —(E-V)=-—. 
2m Ox; at 


Now use Eq. (3.2) (with R for W), multiply 
through by U?, and integrate over the 3n-dimen- 
sional volume; the result is 


S\n=(G) ~2- 
=~ f —var. (3.8) 


For a conservative system, the right side of 
Eq. (3.8) vanishes identically. Our non-conserva- 
tive system will be in a state nearest to con- 
servative if the space average of 0R/dt is a mini- 
mum. We therefore impose the condition 


OR 
-of —Pdr=0. (3.9) 
ot 


Once more we tentatively identify P with UV’. 
The resulting variational problem leads, by the 
usual rules, to the necessary condition 


0?U 2m 





(3.10) 
ax,? 
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subject, of course, to appropriate boundary con- 
ditions. Enter this with U=exp (+R/k) and 
p:=9R/dx;; it becomes — 


V-p=+F2p?/k. (3.5) 


If we use this in Eq. (3.3) we verify the identifi- 
cation of P with U? and complete the derivation 
of the Schroedinger equation (3.10) classically. 
Furthermore, it has been demonstrated that Eq. 
(3.5) is an alternative way of writing Eq. (3.10). 

| Although Eq. (3.10) has exactly the same form 
as Schroedinger’s equation for stationary states, 
and although the dependent variable in both 
instances is a ‘‘probability amplitude,”’ neverthe- 
less there are differences with respect to k, U, and 
V which must be discussed. Concerning k, 
Schroedinger identified hisconstant with Planck’s 
on phenomenological grounds; we do the: same, 
although we hope ultimately to find a logical 
basis for writing k=h. (Cf. the determination of 
K in the astronomical problem in Section II.) 
The first difference then is nil. 

The second difference concerns the admission 
of complex U’s into Eq. (3.10). This will be taken 
up in Section IV below, but here we wish to point 
out that so far as stationary states are concerned, 
the admission of complex y’s into Schroedinger’s 
equation does not seem to be a matter sine qua 
non; and since our theory thus far deals only with 
stationary states, it does not appear that the 
second difference is vital. 

The third difference concerns the fact that our 
dynamical system is non-conservative (i.e., V 
contains £) ; if it were conservative, the right side 
of Eq. (3.8) would vanish identically and our 
procedure would become as meaningless as it was 
in 1926, when Schroedinger? went through a 
similar performance. Two circumstances indicate 
that this third difference is mot of primary signifi- 
cance: (a) Our system is assumed to be in a quasi- 
conservative state in which the energy hovers 
about a definite value and to which we might, 
experimentally, ascribe a conservative character ; 
(b) the Schroedinger equation is not exact. It 
requires correction at least for “‘spin.’’ 

Compared with the outstanding equivalence of 

3If V=V; (conservative)+ V2 (non-conservative), Vs 
can be formally eliminated by taking it over to the right 
of Eq. (3.8) and minimizing it with 0R/dt. This introduces 
a momentum error like that of Sec. II and increases the 
“sufficiency”’ perils. 
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the two equations in form and meaning, these 
differences cannot be considered basic. It is 
submitted that Schroedinger's equation for station- 
ary states ts a consequence of classical physics. 

This violent conflict with accepted theory is 
emphasized by other conflicts; for instance, 
Eq. (3.2) which (without the + sign) was written 
by Schroedinger in 1926, has long since been read 
out of the quantum theory; here it plays a 
fundamental role as the connection between the 
probability and action functions. Of greater sig- 
nificance is a conflict concerning the occupation 
of the discrete stationary states predicted by 
Eq. (3.10); quantum theory asserts that only 
these states are occupiable, to the total exclusion 
of the intermediate continuum of states; our 
classical derivation, although not yet providing 
an “occupation” formula, cannot accept this 
complete exclusivity. It is submitted that there 
is experimental evidence in favor of non-ex- 
clusivity, viz., the overlapping of discrete and 
continuous spectra. 

How, despite these and other conflicts, the two 
theories may nevertheless reach the same or 
similar results, can be glimpsed from the simplest 
example—that of a uniform stream of particles 
moving with constant momentum. Our theory 
simply asserts that the probability in this case is 
constant @ priori; Eq. (3.10) above is not appli- 
cable here because its derivation breaks down 
completely. Quantum theory, however, accepts 
Eq. (3.10) at its face value for this (as any other) 
case, takes a complex exponential solution for U, 
and then washes it out by complex multiplication 
to get a constant probability. 

This example serves also to illustrate the care 
necessary in handling the theory presented here. 
To take a more general case, suppose the mo- 
mentum of the system can be resolved into two 
components pi=grad W,; and p2=grad W2 such 
that div p:+0, div pp=0. Then Eq. (3.5) becomes 


V-pi= F2(pitp2)?/k. (3.11) 


Putting Wi=+k log Ui, the corresponding 
Schroedinger equation becomes 


3? 





= 
+—(pi?+4pi-pe+2p2?)Ui=0. (3.12) 
Ox? k* 


L 


Perhaps we may be able to gain from results of 
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this kind an understanding, in the classical sense, 
of the non-classical rules of quantum theory for 
angular momentum, etc. 


IV. ATOMIC DYNAMICS, CONTINUED 


The conflict mentioned in Section III between 
our theory and quantum mechanics concerning 
the occupancy of intermediate states seems to be 
irreconcilable—although quantum mechanics ap- 
pears to have made some preconcession in the 
form of “‘weak quantization.’”’ The conflict of 
formalisms, per se, is not necessarily irrecon- 
cilable, for it is conceivable that some ingenious 
“frame of reference’’ may be found which will 
make the two formalisms agree, insofar as agree- 
ment is inherently possible. In any case, the 
restriction to real U’s in Eq. (3.10) is incon- 
venient, so we present here the results of an 
attempt to introduce complex U’s in the naive 
way. 

If, instead of transforming the action function 
W by means of Eq. (3.2), complex U’s are intro- 
duced by writing 


W=tk log UU*, 
where 
U=pexp (10), U*=pexp (—10), (», © real), 


and if U and U* are required to satisfy the same 
equation for stationary states, then it is found 
that: (a) U and U* must be related by 


aU 0? U* 
U*y—=UL : (4.1) 
Ox ;? Ox? 





which is equivalent to the following relation be- 
tween the amplitude p and phase 9: 


div. (p? grad 0) =0. (4.2) 


(b) The phase © must appear explicitly in U 
as exp (10). (c) The equation for stationary 
states becomes, in place of Eq. (3.10), 


9 


vu sp? 
= -+(] +070") =o, (4.3) 


Ox ;* 





i.e., the term (grad ©)? is added. 
Equation (4.1) is also equivalent to 


LX [2(dp/dx;)(90/dx;) + pd?O/dx;?]=0, 


so that if }°d°@/dx7 vanishes, the relation (4.1) 
implies that grad p shall be perpendicular to 
grad @. If, for example, in a 3-dimensional prob- 
lem in spherical coordinates, U= u(r, 6) exp (img), 
the resulting ‘“‘wave equation”’ is 


1a0/ ou 1 a _ Ou p? 
—_-— r—)4 sin) += 0, (4.4) 
r? or or r? sin 0 a0 06 





in which the term m?u/r? sin? @ is cancelled, so 
that m does not enter the equation at all; this is 
not surprising since P= UU* is independent of 
m and ¢ in this case. 

The concept of ‘‘smoothest density’ intro- 
duced in Section I for purely mathematical cases 
can be applied equally well to the kinematical 
probability P. Following the previous pattern, 
we define the “‘characteristic’’ or ‘‘smoothest”’ P, 
(call it P.) by the equations 


d 2 
(grad log P)?— rs log P) =(G? 
$ 


3 f GPar=0. 


Using Eq. (2.8), putting P.= U2, and carrying 
out the variation, we find the result is 


ar 








rere 


Inserting the value of div p given by Eq. (3.5) 
this is reduced to Eq. (3.10). Thus the kinematical 
probability P=U? defined by the Schroedinger 
equation (3.10) has the ‘‘smoothest’’ property. 

We conclude with a brief discussion of our 
“quasi-conservative stationary states.’’ They im- 
ply some sort of stabilizing forces which have 
hitherto been lacking on any theory. It will now 
be shown by a simple one-dimensional example 
how this hiatus may be filled by means of an 
electromagnetic wave field. Consider a particle 
moving in a field of force determined partly by a 
conservative potential V;(x), and for the rest by 
the potential of a traveling wave field 


V2(x, t) = A(x) sin 2nr( t— fdx/o(2)+a), 


where v(x) is the phase velocity and a is constant. 
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The equation of motion is 
mé= —OV;/dx—8V2/dx. (4.5) 


This equation is, in general, completely intrac- 
table, but it ceases to be so if the particle travels 
at or nearly at the phase velocity. We may call 
this condition “‘phase-mechanical resonance.”’ 

Suppose first that the particle’s velocity is 
exactly equal to that of the phase, i.e., that 
=v. Then t— f{dx/v is a constant which may 
be put equal to zero. Equation (4.5) now be- 
comes 


mé=—V;'—A’ sina+(2mv/v)A cosa. (4.6) 


Denote the values of x which satisfy Eq. (4.6) 
by y, and put x=y+7 where 7 is supposed so 
small that only its first power need be retained 
in the Taylor expansion of all functions involved 
and in all products of such expansions; also, 
let v vary sufficiently slowly so that /v can be 
written for /dn/v. Forming the expression for 
m(j+7) obtained under these conditions and 
using Eq. (4.6), we find the result is 


myn= —(Vi"+A” sin 2rva—(4rv/v)A’ cos 2rva 
—(2mv/v)?A sin 2rva). (4.7) 


Equation (4.7) shows that if the parenthesis 
on the right is positive, the particle is stabilized 
against small displacements from exact ‘‘reso- 
nance” even if V; and A vary slowly or are 
constant; and Eq. (4.6) shows that under suit- 
able circumstances Iarge interchanges of energy 
can take place between a wave and a particle in 
“‘resonance”’ with it. 

A macroscopic example of this sort of thing is 
furnished by the cyclotron or Lawrence’s linear 
ion accelerator which preceded the cyclotron. 
Better examples, I believe, exist in some elec- 
tronic devices. These examples require more 
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complicated functions for their representation, 
but otherwise they illustrate the phenomenon of 
“rectification” which arises when a particle in a 
wave field “rides with the wave.”’ This concept 
of phase-mechanical resonance goes further than 
supplying stabilizing means for electrons in 
special orbits; it also provides means for under- 
standing such phenomena as electron diffraction 
in terms of classical physics. Indeed, phase 
mechanical resonance represents an inversion of 
physical ideas of quantum mechanics correspond- 
ing to the inversion of mathematical ideas on 
which this paper is based. Thus, quantum 
mechanics ascribes the non-classical behavior of 
electrons to ‘‘matter waves’’ associated with 
electrons; our theory would put the waves as 
an “atmosphere” around the nucleus and ascribe 
to them a classical electromagnetic character. 
In the quantum case, the speed of the particle 
is exactly equal to the group velocity of the 
matter waves; in our case, the speed of a 
particle in phase-mechanical resonance is equal 
or nearly equal to the phase velocity of the 
electromagnetic waves. 

To develop these concepts further, one would 
have to consider stationary waves also, but the 
situation then becomes obscure on account of 
trigonometric terms of double frequency which 
arise. More important, however, is the question 
of how such wave systems can exist at all, for, 
according to ordinary classical theory, they would 
radiate themselves away in short order. For a 
discussion of this question the reader is referred 
to my paper on the “Development of electro- 
magnetic theory for non-homogeneous spaces.’’‘ 
We may look to an extension of classical theory 
of the kind there described for the existence and 
properties of the wave systems here postulated. 


4 B. Liebowitz, Phys. Rev. 64, 294 (1943). 
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Letters to the Editor 








ROMPT publication of brief reports of important dis- 
coveries in physics may be secured by addressing them 
to this department. The closing date for this department is the 
third of the month. Because of the late closing date for the sec- 
tion no proof can be shown to authors. The Board of Editors 
does not hold itself responsible for the opinions expressed by 
the correspondents. Communications should not in general ex- 
ceed 600 words in length. 





Comments on Bethe’s Theory of Diffraction of 
Electromagnetic Waves by Small Holes* 


C. L. PeKeris** 
Columbia University, New York, New York 
November 3, 1944 


N this interesting paper Professor Bethe develops a 

vectorial theory of diffraction of plane electromagnetic 
waves through a circular aperture in an infinite plane con- 
ducting screen, for the limiting case when the aperture is 
small in comparison with the wave-length. In the intro- 
duction Professor Bethe points out that if the problem is 
attacked by Kirchhoff’s method one obtains a solution 
which fails to satisfy the boundary conditions, and this 
result he designates as a “failure of the Kirchhoff theory.”’ 
The method adopted by Bethe is to replace the aperture 
by an equivalent distribution of magnetic currents, for 
which he obtains an explicit solution satisfying the bound- 
ary conditions. 

It should be pointed out that Kirchhoff’s method con- 
sists in making a certain approximation when solving 
Helmholtz’s equation which, as Professor Bethe states, is 
justified in the limit of short wave-lengths. Kirchhoff was 
careful to stress the limitation of his method to wave- 
lengths which are small in comparison with the size of the 
aperture. When he took up the case of diffraction by a 
grating he made the following statement: “The equation 
derived above makes the essential assumption that the 
dimensions of the aperture are very large in comparison 
with the wave-length, and its application to diffraction 
spectra produced by gratings whose lines are only a few 
wave-lengths wide, is not justified.’’ From this statement 
it is clear that Kirchhoff considered his approximation to 
be unjustified in the case where the wave-length is a 
moderate fraction of the aperture, let alone in the other 
extreme case, considered by Bethe, where the aperture is 
very small compared with the wave-length. For this case 
there is no ‘Kirchhoff theory,” and therefore one cannot 
speak of a failure of the Kirchhoff theory. 

A method of solving Helmholtz’s equation in the case of 
long wave-lengths, when Kirchhoff’s approximation cannot 
be made, was given some fifty years ago by Lord Rayleigh.' 
The problem solved by Rayleigh is scalar, but his method 
is essentially the one used by Bethe in the more compli- 
cated vectorial problem. Rayleigh’s results, especially for 
the pertinent case where the wave function is made to 
vanish on the screen,? agree with Bethe’s as far as order of 
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magnitude is concerned. In a later paper* Rayleigh extends 
his method (in the two-dimensional case) to apertures 
which are of the dimension of a wave-length. 

A systematic exposition of the application of the Helm- 
holtz-Kirchhoff theory, including the case of “‘long’’ wave- 
lengths, to the scalar problem will be found in Lamb's 
Hydrodynamics‘ where will also be found a discussion of 
Rayleigh’s paper.’ The question of the retention of the 
various terms in Helmholtz’s equation which Professor 
Bethe raises in Section 2 of his paper is taken up on page 
500. 


*H. A. Bethe, Phys. Rev. 66, 163 (1944). 
** On leave of absence from the Massachusetts Institute of Tech- 


nology. 
1 Lord Rayleigh, Sci. Papers, Vol. IV, p. 283. 
2 See reference 1, pp. 286, 287. 
3 Lord Rayleigh, Sci. Papers, Vol. VI, p. 11. 
“Lamb, Hydrodynamics, pp. 492-520. 
5 See reference 4, pp. 517, 518. 





Double Bragg Reflections of X-Rays 
in a Single Crystal 


C. J. Davisson AND F. E. HAwortH 
Bell Telephone Laboratories, Inc., New York, New York 
October 6, 1944 


HEN a beam of monochromatic x-rays traverses a 
crystal in a direction to satisfy the condition for hkl 
reflection, the resultant Bragg beam traverses the same 
crystal in a direction to satisfy the condition for (hkl), 
reflection. The two beams constitute a symmetrical pair. 
Radiation is reflected out of the primary beam into the 





Fic. 1. Double x-ray reflection in quartz. 


Bragg beam, and out of the Bragg beam into the primary 
beam. There is a continuous interchange of radiation be- 
tween the two beams with a net flow to the less intense. 

In the case in which two Bragg beams are produced 
simultaneously the two beams interchange radiation, not 
only with the primary beam, but also, through associated 
reflections, with one another. If the two reflections are 
Aykyl; and hekele, radiation flows from beam 1 to beam 2 via 
the reflection (h2—/:)(k2—k:)(l2—1:) and from 2 to 1 via 
(hy— he) (ki — 2) (li — 12). These are the associated reflections. 

In the still more special case in which the structure factor 
of one of the reflections is zero, no radiation flows to this 
beam from the primary beam directly, but radiation does 
flow to it indirectly via the other reflection and one of the 
associated reflections. A doubly-reflected beam appears in 
the place of the non-existent once-reflected beam. The con- 
dition for this occurrence is, to repeat, that the primary 
beam shall satisfy the Bragg condition for the non-existent 
reflection, and that it shall, at the same time, satisfy the 
Bragg condition for some allowed reflection. There is the 
further requirement that the associated reflection also is 
allowed. 
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Evidence of what appears to be the occurrence of 
forbidden reflections is met with, not infrequently, in 
electron diffraction patterns, more often in patterns pro- 
duced by single crystals, but in powder patterns as well. 
Among the explanations of these anomalies which have 
been proposed is this one of double reflections.! The possi- 
bility of double reflections in single crystals has thus long 
been recognized. 

We have succeeded recently in demonstrating double 
x-ray reflections in quartz: A beam from a copper target 
tube operated at 24 kv and 25 ma is brought onto an etched 
z-cut (normal to principal axis) quartz crystal at the CuKa 
Bragg angle for the non-existent 00-2 reflection (16° 36’). 
A film in a semicircular cylindrical holder is set to receive a 
regularly reflected beam through a 140-degree paper- 
covered window in the holder’s inner wall. The crystal and 
holder rotate together at a uniform rate about an axis 
normal to the crystal surface (z axis). The axis of the holder 
coincides with the axis of rotation, and the primary beam 
meets the crystal on this axis. Half of the time the crystal 
is eclipsed by the holder which is receiving the primary 
beam on its outer wall; the other half, the primary beam, 
reaches the crystal, and any beam which emerges in the 
direction of regular reflection, while the window is in 
opposition, reaches the film. At certain angular positions of 
the crystal the primary beam satisfies not only the condi- 
tion for 00-2 reflection but for some other reflection (or 
reflections) as well. As the rotating crystal passes through 
these positions, the doubly-reflected beams flash out. Those 
which reach the film produce a characteristic pattern, an 
example of which is reproduced in Fig. 1—exposure time 
72 hours. 

The chief features of the pattern are pairs of strong 
spots. The members of the pair are separated by a little 
more than 8 degrees, and the pairs themselves by 60 
degrees. These are produced by double reflections of the 
type (01-1+01-1). The spot midway between the strong 
pairs is itself a close pair,each member of which is produced 
by two simultaneous double reflections of the type 
(01-3+01-1) and (01-1+01-3). The identities of the other 
spots have not yet been worked out. Most of them will be 
found to have been produced by the simultaneous double 
reflections. 

The complete circular pattern has 6-fold symmetry about 
its center, and six 2-fold diametric axes. The number of 
spots discernible in a 60-degree interval of the original of 
the pattern here shown is 18, making 108 in the whole 
pattern. The total number of flashes per complete rotation 
can be calculated from the density of points in the quartz 
reciprocal lattice and other considerations—the number is 
about 700. The intensities of most are below the threshold 
of the present observations. 

The masking effects of the generally scattered radiation 
from the crystal and of the 00-3 Laue beam from the con- 
tinuous spectrum are suppressed by shielding and filtering; 
that of the 00-6 and higher order Laue beams is avoided by 

limiting the tube voltage to exclude the 00-6 wave-length 
from the spectrum. 

A number of interesting questions present themselves. 
There is, to mention but one, the question of the phase 
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relation between the two beams produced by simultaneous 
double reflections. These beams are coherent and super- 
posed, and must therefore interfere. 


1H, Raether, Zeits. f. Physik 78, 527 (sen; J. A. Darbyshire and E, 
. Cooper, Proc. Roy. Soc. A152, 104 (1935) 





Enhanced Thermionic Emission 


J. B. JoHnson 
Bell Telephone Laboratories, Incorporated, New York, New York 
November 1, 1944' 


NEW type of electron emission from oxide-coated 

thermionic cathodes is disclosed by a method of 
measuring simultaneous thermionic and secondary emis- 
sion. In this method the heated oxide target is bombarded 
by the primary electrons from a gun in which the high 
voltage electron stream can be turned on by a small 
auxiliary voltage. The steady collector voltage is first 
applied so that before the bombardment starts the circuit 
is in electrical equilibrium. The thermionic current is then 
read by a d.c. meter. Short-time changes in the current to 
or from the target are observed on an oscilloscope with an 
intervening amplifier from which d.c. is excluded. Time 
as short as 2X 107? sec. can be resolved. 

When now the primary bombardment is started the 
oscilloscope shows at first an abruptly rising flow of 
electrons from the target, followed by a gradual rise to an 
equilibrium value many microseconds later. When the 
bombardment is stopped, the electron current drops 
abruptly by the same amount as the initial rise, and there 
then follows a more slowly decreasing current. The bom- 
barding primary current shows none of this behavior but 
begins and ends sharply. 

The initial rise or final drop of current is interpreted as 
the beginning or end of normal secondary emission of the 
target. It changes little with temperature and varies in 
the usual way with the primary voltage. At a primary 
voltage of 1200 to 1500 volts, the ratio of the maximum 
value of the secondary current to primary current varies 
from 3 to 5. The emission which persists after the end of 
the bombardment can hardly be secondary emission but 
must be of thermionic origin, and presumably the equal 
rise of emission during the bombardment is of the same 
kind. This emission varies with the temperature of the 
target in about the same way as the steady emission, thus 
following roughly the Richardson law. It increases with 
bombarding voltage and current density, and may exceed 
the steady thermionic current in value. It is undoubtedly 
an enhanced thermionic emission excited by the electron 
bombardment. 

A natural assumption is that the increased emission is 
caused by a rise in target temperature caused by the bom- 
bardment. The temperature rise of the surface of the target 
may be calculated, and is far too small to explain the effect. 
One must conclude instead that the bombardment tem- 
porarily changes the thermionic activity of the oxide target. 

This effect no doubt explains the exponential rise with 
temperature that has been reported for the secondary 
emission factor of oxide cathodes." 


P 1N. Morgulis and A. Nagorsky, J. Tech. Phys. U.S.S.R. 5, 848 
1938). 
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Proceedings of the American Physical Society 


MINUTES OF THE MEETING OF THE NEW ENGLAND SECTION AT PROVIDENCE, OCTOBER 28, 1944 


HE twenty-fourth regular meeting of the New England Section of the American Physical 
Society was held at Brown University in Providence, Rhode Island, on Saturday, October 28, 
1944. At least sixty members were in attendance. The invited papers were as follows: 


The British Society for Freedom in Science. P. W. BripGMan, Harvard University. 
The Applied Mathematics Program at Brown University. DEAN R. G. D. Ricnarpson, Brown 


University. 


The Electrical Properties of Polarized Metal-Solution Interfaces. Davip C. Graname, Depart- 


ment of Chemistry, Amherst College. 


Some Metallurgical Applications of Radioactive Tracers. Joun W. IrvINE, JR., Department of 
Chemistry, Massachusetts Institute of Technology. 
Continuity in Mathematics and Physics. Morton Masius, Worcester Polytechnic Institute. 


Abstracts of some of the contributed papers are appended; others will appear in the American 


Journal of Physics. 


At the business meeting the following officers were elected for the ensuing year: 


Chairman, Morton Masius, Worcester Polytechnic Institute. 

Vice Chairman, R. B. Linpsay, Brown University. 

Secretary-Treasurer, MILDRED ALLEN, Mount Holyoke College. 

Members of Program Committee, V. E. Eaton, Wesleyan University; Francis W. SEARs, 


Massachusetts Institute of Technology. 


Conversion of Non-Rationalized c.g.s. to Rationalized 
m.k.s. Units in Electromagnetism. HERBERT JEHLE, Har- 
vard University.—The standard method of converting the 
numerical values of a physical quantity from one unit to 
another is the substitution of units whereby a physical 
quantity is considered as an invariant with respect to 
change of units, e.g., energy=37 erg=37 (10-7 joule) 
= (37.10-") joule. The conversion of the numerical values 
of E, B, H, and of nearly all other electromagnetic quanti- 
ties may be performed by such a substitution. With the 
usual choice of units, however, this conversion procedure 
breaks down for electrostatic flux and for D, because the 
rationalization of Gauss’ law was performed by a change 
of the physical quantity of D by a factor 4x. This of 
course invalidates the substitution method. The most 
convenient way to solve this dilemma is to keep con- 
sistently to the substitution method for converting units; 
and to choose the units of electrostatic flux in c.g.s., non- 
rationalized m.k.s., and rationalized m.k.s., respectively, 
as: 1 statlorentz=1 lorentz/3-10®°=1 coulomb/4xr-3- 10°. 
This proposed introduction of a new c.g.s. unit, the 
“statlorentz,”" does not involve any change of numerical 
values; and therefore leaves all physical equations un- 
touched, as it is the numerical values only which enter 
into the equations. 


Apparatus for the Production and Measurements of 
Small Motions. Gorpon FERRIE HULL, JR., Bell Telephone 
Laboratories, Inc. (now at Dartmouth College).—A canti- 
lever bar is described for producing small motions of the 
order of 10-* cm and greater together with an inter- 
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MILDRED ALLEN 
Secretary-Treasurer 


ferometer for measuring these motions. Two methods of 
using the interferometer are discussed. The first uses an 
uninterrupted monochromatic light source in which the 
interference fringes are observed to be visible and sta- 
tionary for particular amplitudes of oscillation. The second 
method uses a stroboscopic light source synchronized with 
the oscillating bar. Motions of the cantilever bar measured 
by these two methods are compared with each other and 
with the theory of the cantilever bar. 


Acoustic Intensity Distribution from a Piston Source. 
A. O. WiLuiaMs, JR., AND L. W. LaBaw, Brown University. 
—The acoustic intensity distribution from a circular 
piston-like source, earlier treated by several investigators, 
is evaluated by an approximate method. The velocity 
potential, expressed in cylindrical coordinates about the 
beam (and piston) axis, is expanded in powers of r, with 
coefficients depending on the known intensity variation 
along the z axis. The first few terms of the series are com- 
puted under certain restrictions of wave-length and dimen- 
sions, resulting in a simple, reasonably accurate approxi- 
mation particularly applicable to narrow supersonic beams. 
The peak of the intensity curve is a parabola of calculable 
curvature. This expression permits determination of the 
cone on the surface of which the intensity falls to one-tenth 
the axial value. Certain data taken in this laboratory with 
radiometers and microphones (in water at frequencies 
near 1 megacycle) allow quantitative comparisons. The 
general agreement is satisfactory, excelling that resulting 
from two other approximations used earlier. Since the 
theory discusses point values, while microphone and 
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radiometer measure average values over an area, correc- 
tions to the latter are necessary. Simple corrections are 
derived and shown to improve the agreement. 


Wave Front Determination in a Supersonic Beam. L. W. 
LaBAW, Brown University——An experimental method has 
been developed by which the loci of equal phase, in a 
unidirectional high frequency compressional wave beam in 
liquids, may be accurately determined without recourse to 
a standing wave pattern. This permits the tracing of wave 
fronts in the beam at different distances from the supersonic 
source, furnishing a new tool for the study of elastic wave 
diffraction in liquids. The electric signal generated by the 
piezoelectric action of a Rochelle salt microphone which 
has the phase of the compressional wave in the liquid is 
modulated by a small electric signal which has the phase of 
the voltage impressed across the quartz supersonic gener- 
ator. The sum of these two signals with the proper ampli- 
tudes gives a sinusoidal space variation of the root-mean- 
square value of the high frequency signal generated by 
the microphone as it is moved in the direction of propaga- 
tion of the beam. The locus of points of equal phase may 
thus be found by following a maximum in this sinusoidal 
space variation as the microphone is moved at right angles 
to the direction of propagation. This has the advantage 
over any optical method since it permits the wave front 
determination in a beam of circular symmetry. 


Research on the Origin of Petroleum. CLARK GOODMAN, 
Department of Physics, Massachusetts Institute of Technol- 
ogy.—During the past two years the American Petroleum 
Institute has sponsored research on the transformation of 
organic materials into petroleum. Under the general 
supervision of technical representatives from all of the 
contributing oil companies, investigations are in progress 
at the Scripps Institution of Oceanography on the role of 
micro-organisms, at Pennsylvania State College on allied 
bacteriological research and on chemical aspects of the 
problem, and at the Massachusetts Institute of Technology 
on the effects of terrestrial radioactivity. The objective of 
this project is to determine the origin or origins of petro- 
leum in the earth. In addition to being of general scientific 
interest, it is anticipated that this study will improve the 
methods of exploration for petroleum. The increased de- 
mands for new petroleum reserves imposed by a war 
economy emphasize the need for such improvements. 


A Theory of Singing Flames. ARTHUR TABER JONES, 
Smith College-—A gas flame may sing when surrounded by 
a resonant air tube. Rayleigh’s theory assumes standing 


waves in the gas supply tube, with most rapid efflux of gas 
approximately a quarter of a period before maximum com- 
pression in the air tube. It is now found that flames may 
sing without standing waves in the gas tube. In such cases 
any vibration in the air tube causes progressive waves to 
run backward through the gas tube, so that the phase of 
maximum efflux of gas is approximately that of minimum 
pressure in the air tube. The present work makes it likely 
that when a flame sings there is often a delay of about 
half a period between greatest rate of efflux and most rapid 
supply of Reat, thus providing the proper phase relation- 
ship for maintenance. When the gas tube ends below in a 
vessel of some size, standing waves are superposed on the 
running waves, and if the flame is small and hot, like those 
employed by Sondhauss, the frequencies maintained show 
gaps that do not occur when the flame is larger. 


The Energy Spectrum of the Primary Cosmic Radiation. 
SHuICHI KusaKA, Smith College.-—The experimental results 
of Millikan, Neher, and Pickering on the energy incident 
vertically and the total incident energy at different mag- 
netic latitudes in India and in North America are com- 
pared with the theoretical calculations based on the geo- 
magnetic theory of Lemaitre and Vallarta. It is found that 
the agreement between theory and experiment is within 
the statistical error if the distribution in number of the 
primaries is taken to be a curve which follows the inverse 
cube law for high energy and which flattens out on the low 
energy side at about 6 Bev. Hence this energy distribution 
is offered as a possible alternative to the discontinuous one 
suggested by Millikan and his colleagues based on their 
atom annihilation hypothesis. 


The Deflection Sensitivity of Parallel-Wire Lines in 
Cathode-Ray Oscillographs. H. GUNTHER RUDENBERG, 
Cruft Laboratory, Harvard University—The deflection 
angle @ of a cathode-ray beam accelerated by a potential 
V is shown to be proportional to the linear charge density p 
transverse to the beam on the interior of the electrostatic 
deflection system, 0=p/2 Veo. The voltage deflection sensi- 
tivity depends on the plate-to-plate capacitance ¢ per 
unit width, 6/E=c/2Ve . This includes edge effects and 
permits a determination of the sensitivity from capacitance 
measurements or a simple graphical evaluation from the 
dimensions of the deflection system. The use of parallel- 
wire lines as a deflection ‘system is suggested, and its 
deflection sensitivity evaluated. Such lines have low 
transit time error and are either easily tuned, or can be 
operated non-resonant with resistive termination. 
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MINUTES OF THE MEETING OF THE NEW YORK STATE SECTION AT 
SCHENECTADY, NOVEMBER 11, 1944 


HE Fall Meeting of the New York State Section of the American Physical Society was held 

in the Physics Building, Union College, Schenectady, New York on Saturday, November 11, 

1944. This meeting was sponsored jointly by the Section and the New York State Science Teachers 

Association. The attendance at the sessions numbered over one hundred and fifty. In view of this 

large attendance the resolutions committee recommended that joint meetings of the section and the 

Science Teachers Association receive consideration in the planning of future meetings. The following 
program of invited papers was presented. 


Morning Session 


Address of Welcome. D1xon Ryan Fox, President, Union College. 

Physics in the Secondary Schools of New York State. E. L. MANNING, State Education De- 
partment. 

Presentation of A.C. Phenomena for High School Students, Demonstrations. P. I. WoLpD, 
Union College. 

Applications of Electrical Contact Resistance, Polarized Light, and Electron Diffraction 
Techniques to the Study of Aluminum Surfaces for Spot Welding. R. A. Wyant, G. H. 
CARRAGAN, AND K. H. Moore, Rensselaer Polytechnic Institute. 


Afternoon Session 


Business Meeting. 

Demonstrations for Physical Optics, Using Very Short Radio Waves. C. L. ANDREws, New 
York State College for Teachers, and General Electric Company. 

Night Vision. Brian O'BRIEN, University of Rochester. 

New York State Aurora and Geomagnetism. C. W. GARTLEIN, Cornell University. 


The biennial election of officers of the Section was reported at the afternoon session and is as 


Terms expire 


Chairman, P. R. GLEASON, Colgate University, 1946. 
Vice Chairman, J. R. CoLutns, Cornell University, 1946. 
Secretary, W. R. FREDRICKSON, Syracuse University, 1946. 
Treasurer, G. H. CAMERON, Hamilton College, 1946. 
Members of Executive Committee, 
C. K. BoswortH, Glens Falls High School, 1948. 
V. H. RECKMEYER, Agfa Ansco Corporation, 1948, 


The Executive Committee also includes three others whose terms expire in 1946: W. B. Rayton, 
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MINUTES OF THE MEETING OF THE OHIO SECTION IN CONJUNCTION WITH THE 
A.A.A.S. AT CLEVELAND, SEPTEMBER 11-12, 1944 


HE American Physical Society was represented at the Cleveland meeting of the American 
Association for the Advancement of Science by its Ohio Section, which held in conjunction 
with Sections B and E of the Association the meetings of which the programmes and some abstracts 
are hereto appended. 
LEON E. SMITH 
Secretary 


Contributed Papers 


Equipment and Techniques for Evaporation of Metals. Leonarp O. O_sen, CHARLEs S. 
SMITH, AND E. C. CRITTENDEN, JR., Case School of Applied Science. ] 

The Evaporation of Metals from Crucibles. CHARLEs S. Smita, LEONARD O. OLSEN, AND E. C. ( 
CRITTENDEN, JR., Case School of Applied Science. 

A Convenient Virtual Image Locator for Elementary Optics. Pau. K. TayLor, Department of 
Physics, Oberlin College. : 

The Energy Balance at the Positive Crater of the Carbon Arc. H. G. MAcPuHersoN, National : 
Carbon Company. 

Calculation of the Piezoelectric Effect in Ionic Lattices of the Zinc Blende Type. HANns JAFFE, : 
The Brush Development Company. 

Non-K-Electron Capture by Nuclei of Relatively Low Isotopic Number. ELtiot Q. ADAms, . 
General Electric Company. , 

fc 


Invited Papers 

Medical Physics. Orro GLAsseR, The Cleveland Clinic. ol 

Biophysics in Modern War and Technology. DetLEv W. Bronk, The Eldridge Reeves Pp 

Johnson Foundation, University of Pennsylvania. 

Mathematical Biophysics in Its Relation to Experimental Biology. N. RasHevsky, University 

of Chicago. S. 

C 

Invited Papers on Physics in the Metal Industry ue 

Selection and Training of College Students for Industrial Research. A. W. Hut, General ev 

Electric Company. , A 

The Strength Properties of Metals. GeorGce Sacus, Case School of Applied Science. ha 

Recovery Phenomena in Metals of Continuous Casting. N. P. Goss, S.O.M. Center. we 

Heat Transfer Problems in the Foundry. H. A. Scuwartz, National Malleable and Steel by 

Castings Company. = 

Metallurgical Developments in the Production of Steel Castings. CHarLeEs W. Briccs, Steel is 

Founders’ Society of America. up 

The Spectroscopic Laboratory on the Furnace Charging Floor. J. T. Rozsa, Republic Steel tui 

Cor paration. to 

at 

. ‘ to 

Quartz Symposium ot 

Joint Meeting with Section E (Geology) of the A.A.A.S. 

alu 

Quartz for Radio Use—1944. ALLYN C. SWINNERTON, Office of Chief Signal Officer. the 

Quartz as a Mineral Resource. Ropert B. McCormick, War Production Board. Su, 

Geology of Quartz Deposits Outside of the United States. Cart ToLMAN, Foreign Economic dif 
Administration; RoBert D. BuTLer, Foreign Economic Administration in Brazil; RICHARD 

E. Storper, Office of Chief Signal Officer. d 

Geological Investigations of United States Quartz Crystal Deposits with Special Reference to uP 

0c: 


Those of Arkansas. A. E. J. ENGEL AND OTHERS, Department of Interior. 
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Frequency Adjustment of Quartz Oscillator-Plates by X-Rays. CLirrorp FRONDEL, Reeves 


Sound Laboratory. 


Twinning in Quartz. CLIFFORD FRONDEL, Reeves Sound Laboratory. 
The Interlacing of Piezoelectricity with Other Branches of Science. W. G. Capy, Wesleyan 


University. 


The Quartz Crystal Unit. R. A. Hetsinc, Bell Telephone Laboratories. 
The Aging of Quartz Crystal Plates. Vinci. E. Borrom, Camp Coles Laboratory. 
X-Ray Studies of Surface Layers of Crystals. EvizabetH J. ArmstronG, Bell Telephone 


Laboratories. 


The Continuing Development of Quartz Crystal Units. Kari S. VAN Dyke, Office of Chief 


Signal Officer. 


“Crystals Go to War.” A thirty-minute sound and color film. This will be exhibited through the 
courtesy of the Reeves Sound Laboratory, Inc. The motion picture shows the steps in the 
manufacture of crystal units in the Reeves plant. It is expected this will be the first public 
showing of a complete revision of the 1943 film of the same title. 


Equipment and Techniques for Evaporation of Metals. 
LEONARD O. OLSEN, CHARLES S. SMITH, AND E. C. 
CRITTENDEN, JR., Case School of Applied Science.—An 18- 
inch evaporation system has been designed and constructed 
along fairly standard lines. An externally operated, four- 
inch vacuum valve between the diffusion pump and the 
bell jar allows air to be admitted to the bell jar without 
cooling the diffusion pump. A large number of evaporations 
have been carried out using close-wound conical baskets 
and/or helical coils of nickel, iron, chromel, and tungsten, 
as well as crucibles (described in the following abstract) as 
the heating elements. Forty metals have been observed 
for evaporation behavior with as many of the above- 
mentioned sources as was deemed feasible. The information 
obtained supplements the excellent work previously re- 
ported by Caldwell.! 

1W. C. Caldwell, J. App. Phys. 12, 779 (1941). 


The Evaporation of Metals from Crucibles. CHar.Es S. 
SmitH, LEONARD O. OLSEN, AND E. C. CRITTENDEN, JR., 
Case School of Applied Science.—Metals which alloy readily 
with those metals ordinarily used for filaments cannot be 
evaporated in large quantities by the conventional methods. 
A technique for preparing crucibles for use in such cases 
has been developed. Ten- to 25-mil tungsten wire is close 
wound into a conical basket and cleaned and stress relieved 
by flashing in vacuum. The basket is then painted with a 
suspension of 325 mesh BeO in HCI (Ph=1). The crucible 
is air dried and repainted as often as necessary to build 
up a uniform and continuous coat of BeO bridging the 
turns of wire. The product is then slowly heated in vacuum 
to 1000°C by means of the tungsten element. After holding 
at 1000°C for 15 minutes and cooling the crucible is found 
to be hard and quite strong. Such crucibles have been made 
with a 1-cm’ capacity and used at 1750°C for short periods, 
and at 1550°C for long periods. A water suspension of fine 
alumina used in the same way is somewhat more convenient 
than BeO and is excellent for temperatures below 1200°C. 
Such alumina crucibles have been given repeated use with 
different metals. 


A Convenient Virtual Image Locator for Elementary 
Optics. Pau. K. Taytor, Oberlin College.—The instrument 
locates a virtual optical image with reference to the position 





of a real image on a ground glass plate, thus avoiding the 
parallax method which, for the beginning student, is often 
confusing and inaccurate. It is adaptable to experiments 
either with lenses or with plane and spherical mirrors. 
Correct placement is easily and quickly made with re- 
settings within one millimeter a usual occurrence. 


The Energy Balance at the Positive Crater of the Carbon 
Arc. H. G. MacPuerson, National Carbon Company.— 
Measurements were made of the electrical energy input to 
the crater of a low current carbon arc between solid elec- 
trodes and of the losses by radiation from the crater surface, 
by conduction through the electrode, and by evaporation 
of carbon. For a }” graphite anode at 10.5 amperes, 24 
percent of the energy was radiated from the crater, 10 
percent was lost by evaporation of the carbon, and 66 per- 
cent was conducted back along the electrode. In the case 
of a #" carbon anode at 12 amperes, 36 percent of the 
energy was dissipated by radiation from the crater, 11 
percent by evaporation, and 53 percent by conduction 
away from the crater. The relationships between the anode 
voltage drop, the temperature and area of the anode spot, 
and the physical characteristics of the anode material are 
discussed in the light of these results. 


Calculation of the Piezoelectric Effect in Ionic Lattices 
of the Zinc Blende Type. Hans Jarre, Brush Development 
Company.—Except for some early attempts around the 
year 1900, nd derivations of piezoelectric coefficients from 
atomic data alone have been published. As a first step, 
the piezoelectric coefficient e:4 (ratio of piezoelectric polar- 
ization to applied shearing strain) is calculated for an 
ionic crystal of the zinc blende structure type. The assump- 
tion is made that the applied shearing strain will leave 
the distance between nearest neighbor ions unchanged. 
The relative displacement between positive and negative 
ions is thereby completely determined and the resultant 
polarization is P,=¢,4-x,=ne/a*-x,, where ne is the charge 
per ion and a the lattice constant. The assumption is 
justified as a starting point by the successful use of fixed 
ionic radii in estimating stability ranges of structure types. 
For the crystal zinc blende, ZnS (a@=5.43A), we find 
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é14=1.08 coulomb/m* for »=2. From the best experi- 
mental value! of the piezoelectric compliance, di4=3.25 
-10-" coulomb/newton, and the elastic compliance s44, we 
find ¢:4=0.142 coulomb/meter*®. The polarity of the ob- 
served effect agrees with our derivation for the ionic 
lattice. As a completely covalent structure would lead to 
the opposite charge distribution, we conclude that the bond 
in ZnS is somewhat more ionic than covalent. 
1K. S. Knol, Proc. Amst. Akad. Sci. [1] 35, 99 (1932). 


Non-K-Electron Capture by Nuclei of Relatively Low 
Isotopic Number. ELtior Q. Apams, General Electric 
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Company.—Radioactive transformations in which nuclei 
of relatively low isotopic number capture an electron are 
customarily designated as “‘K capture.’’ Based on the 
table of Hartree and Hartree for the self-consistent field 
of mercury, it is estimated that only 87.6 percent of the 
electrons captured by such nuclei of mercury are from the K 
shell, the contribution from the outer shells being, respec- 
tively: L, 9.6 percent; M, 2.1 percent; N, 0.5 percent; 
O, 0.1 percent; and P, 0.006 percent. Since only S electrons 
have any appreciable probability of being found at the 
nucleus, the process is more properly called ‘‘S capture.” 
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